{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "480765be",
   "metadata": {},
   "outputs": [],
   "source": [
    "from astropy.io import fits\n",
    "import HAWCtoHerschel as HtoH"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "d5477674",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Source information\n",
    "SourceName = 'Fil5'\n",
    "FullName = 'Filament 5'\n",
    "# Loading the fits files\n",
    "SourceT = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Tdust_10_10.fits')\n",
    "SourceN = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_N_10_10.fits')\n",
    "SourcePol = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Gal.fits')\n",
    "SourcePlanck = fits.open('../FIELDMAPS_Plots/'+SourceName+'/Data/'+SourceName+'_Planck.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "8fbf97a2",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Plot parameters\n",
    "ScaleLength = 0.0139746 # 0.0168517 # 1 pc in degrees for the target\n",
    "\n",
    "# Region location, width, and height\n",
    "l_center = 18.565 # Position in galactic longitude (degrees)\n",
    "b_center = -0.11 # Position in galactic latitude (degrees)\n",
    "l_width = 0.415 # Width in galactic longitude (degrees)\n",
    "b_height = 0.225 # Height in galactic latitude (degrees)\n",
    "Region = [l_center,b_center,l_width,b_height]\n",
    "\n",
    "# Intensity contours for the plot\n",
    "Icontours = [10, 20, 30, 40, 50, 60] # in mJy/arcsec^2\n",
    "# Ranges for plots\n",
    "Nscale = [0,1.5e+22] # Range for column density in cm^-2\n",
    "Tscale = [15.0,25.0] # Range for temperature map in K\n",
    "Iscale = [0.0,20.0] # Range for intensity map in mJy/arcsec^2\n",
    "# Vector scaling fraction in polarization plots\n",
    "Pvscale = 0.6\n",
    "\n",
    "# Figure size\n",
    "MapFigSize = [10.67,6] # in inches\n",
    "\n",
    "# Effective HAWC+ beam size to use\n",
    "HAWCBeam = 0.00519 # 18.7'' in Band E"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "a2a163b6",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Beam Size of polarization data\n",
      "18.2\n",
      "Beam size of target pojection\n",
      "36.395999999999994\n",
      "Pixel size of the polarization data\n",
      "4.550000000000001\n",
      "Pixel size of the target projection\n",
      "14.0000004\n",
      "Target standard deviation for smoothing in arcseconds\n",
      "13.384761981071023\n",
      "Target standard deviation for smoothing in pixels\n",
      "2.941705929905719\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: nan_treatment='interpolate', however, NaN values detected post convolution. A contiguous region of NaN values, larger than the kernel size, are present in the input array. Increase the kernel size to avoid this. [astropy.convolution.convolve]\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:96: RuntimeWarning: invalid value encountered in divide\n",
      "  dPI = ((StokesQ*dQ)**2.0 + (StokesU*dU)**2.0)**0.5/PI_biased\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:106: RuntimeWarning: invalid value encountered in divide\n",
      "  P_biased = 100.0*PI_biased/StokesI\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:109: RuntimeWarning: invalid value encountered in divide\n",
      "  dP = np.absolute(P_biased*((dPI/PI_biased)**2.0 + (dI/StokesI)**2.0)**0.5)\n",
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:120: RuntimeWarning: invalid value encountered in divide\n",
      "  dO = (0.5*180.0/math.pi)*((StokesQ*dU)**2.0 +\n"
     ]
    }
   ],
   "source": [
    "# Reprojecting the HAWC+ Stokes data to the Herschel pixel scale\n",
    "NewPol = HtoH.ReprojectCube(SourcePol,SourceN, Smooth=True)\n",
    "#NewPol.writeto('test2.fits', overwrite=True)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "31eb9087",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "c:\\Users\\simon\\OneDrive\\Python\\FIELDMAPS_Testing_Herschel\\HAWCtoHerschel.py:228: RuntimeWarning: invalid value encountered in divide\n",
      "  imask_01 = np.where(Pdata[0].data/Pdata[1].data < IdI) # Total intensity SNR threshold\n"
     ]
    }
   ],
   "source": [
    "# Creating downsampled vector catalogs (independent vectors)\n",
    "# Creating the mask for the region\n",
    "DownsampleMask = HtoH.DownsampleVectors(NewPol,Step=2)\n",
    "# Creating the new catalogs for the region\n",
    "CatI, CatdI, CatQ, CatdQ, CatU, CatdU, CatP, CatdP, CatO, CatB,CatdO, CatPI, CatdPI, CatN, CatT = HtoH.MakePCats(SourceN, SourceT, SourcePol, NewPol, CatMask=DownsampleMask)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "2c4fb0cb",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "===\n",
      "Curve fit test\n",
      "Log scale fit: (-0.7384879360515205,33.47564164912725)\n",
      "Direct power law fit: (-0.8385745821176855,43.1184552246998)\n",
      "===\n",
      "\n",
      "Power law index: -0.7384879360515205 ± 0.043177236901130275\n",
      "Coefficient: 33.47564164912725 ± 3.044600805968337\n",
      "Chi-Squared: 259.1271200235861\n",
      "Number of elements: 134\n",
      "Reduced Chi-Squared: 1.963084242602925\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfUAAAF2CAYAAAB3bAFnAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB2E0lEQVR4nO3deVxU1fsH8M8M+zqCyDLIIqigooC4ZGqammbuKV+yVEDUNBQRzaXMBU3TX26IaRqCpiZhamVuRS5lriiuiaIoKJuIrAICc39/3GZkYIBZmYXn/XrdlzN37j33ucM4z5xzzz2HwzAMA0IIIYRoPa66AyCEEEKIclBSJ4QQQnQEJXVCCCFER1BSJ4QQQnQEJXVCCCFER1BSJ4QQQnQEJXVCCCFER1BSJ4QQQnQEJXVCCCFERzSbpF5RUYFly5ahoqJC3aEQQp9HonHoM6kbNCKpb926FV26dIGlpSUsLS3Rq1cvHDt2TPR6eXk5QkND0bJlS5ibm2Ps2LHIycmR6RgVFRVYvnw5fWCJRqDPI9E09JnUDRqR1Fu3bo2vvvoKSUlJuHLlCgYMGIBRo0bh9u3bAIA5c+bg119/RUJCAs6cOYPMzEy8//77ao5aPbZs2aITx1S0THn2l2UfabdtbDt1/L2aGn0m5d9f2Z9JZW2j7XThMyl3eYyGsrKyYr777jumoKCAMTAwYBISEkSv/fvvvwwA5vz581KXV1hYyABgCgsLVRFuk+nQoYNOHFPRMuXZX5Z9pN22se3qe11XPo8MQ59JRfZX9mdSkW3oM6lZx5S3PH3l/a5QjurqaiQkJKC0tBS9evVCUlISKisrMWjQINE2np6ecHZ2xvnz5/HGG29ILKeiokKsGamgoAAAUFhYqNL4Va26uhpFRUVaf0xFy5Rnf1n2kXbbxrar73Xhuqb+W6oCfSbl31/Zn0lFtqHPpOYck2EYVFRUoKCgAFzu6wZ1IyMjGBkZNbqzRrhx4wZjZmbG6OnpMTwej/ntt98YhmGYvXv3MoaGhnW27969OzN//vx6y1u6dCkDgBZaaKGFFlp0Ylm6dGmjuVRjauoeHh5ITk5GYWEhDhw4gMDAQJw5c0bu8hYtWoSIiAjR88LCQjg7OyMjIwOWlpb171hSAri7A+XlwNmzgLe33DHogqqqKjx9+hSOjo7Q19eYjwshhOisoqIiODk5IT09HTweT7S+0Vo6AI35ljY0NETbtm0BAH5+frh8+TI2bdqEgIAAvHr1CgUFBWjRooVo+5ycHNjb29dbXn3NFMIe9vWytATeew84eBA4fhzo21fuc9IFVVVVKCoqgqWlJSV1QghpQjwer+F8JYFG9H6XRCAQoKKiAn5+fjAwMEBiYqLotZSUFKSnp6NXr16qOfi4cey/Bw4ADKOaYxBCCCFKphFVr0WLFmHo0KFwdnZGcXEx9u3bh9OnT+PEiRPg8XgICQlBREQErK2tYWlpiVmzZqFXr171dpJT2LBhgKEhcO8ecPs24OWlmuMQQgghSqQRST03NxeTJk1CVlYWeDweunTpghMnTuCdd94BAGzYsAFcLhdjx45FRUUFhgwZgm+++UZ1AVlaAkOGAL/+ytbWKakTQgjRAhyGaR7ty0VFReDxeCgsLJTuGsXu3UBgINCxI1tbb6aqqqrw5MkTtG7dmq6pE0JIE5A5X9WgsdfU1W7kSMDAALhzB7h1S93REEIIIY2ipF6fFi2AoUPZxz/8oNZQCCGEEGlQUm/I+PHsv/v3Uy94QgghGo+SekNGjABMTYGHD4FLl9QdDSGEENIgSuoNMTMDRo1iH1MTPCGEEA1HSb0xwib4+Higulq9sRBCCCENoKTemCFDACsrIDsbUGAsekIIIUTVKKk3xtDw9bCxzawJXiAQoKCgAAKBQN2hEEIIkQIldWkIm+B/+gl49Uq9sTQhgUCAwsJCNJPxiQghROtRUpfGW28BDg7AixfAiRPqjoYQQgiRiJK6NPT0gIAA9vG+feqNRcNQEz0hhGgOSurSEjbB//ILUFqq3lg0CCV1QgjRHJTUpdW9O+DuDrx8ySZ2QgghRMNQUpcWhwN88AH7uJn1gieEEKIdKKnL4sMP2X+PHwfy89UbSxPJyclBVFQUMjMz1R0KIYSQRlBSl0XHjkCXLkBlJXDwoLqjaRLx8fFISUlBTEyMukMhhBDSCErqshJ2mNu7V71xNJGAgAB4eHggJCRE3aEQQghpBCV1WX34IXt9/fRp4PFjdUejcnZ2dggLCwOfz1d3KIQQNYiLiwOHwwGHwxGt69+/PzgcDoKCgtQXGJGIkrqsnJ2Bt99mH+/Zo95YCCFEAcLkLGk5fPgwAKBVq1bo2bMnevbsqd5glSgoKAgcDgf9+/eXavv63qPFixerNlA56Ks7AK00aRLw55/A7t3AZ5+xNXfSIIFAgKKiIlhaWoLLpd+ShGgSQ0ND+Pr6iq2ztrYGAAwbNgzDhg1TR1gax8fHB0ZGRqLnTk5OaoxGMvp2lcf77wOmpsC9e8ClS+qORivQIDWEaC4HBwdcuHBBbHnrrbcASG5+lyQwMBDt2rWDhYUFDA0N4eLigrCwMBQVFYm2EdaQXV1dsWvXLjg5OcHS0hLh4eEoKytDeHg4eDwenJ2dsW3bNrHyMzMzMXnyZPD5fBgaGsLNzQ0rVqxAVVWVaBthy8OkSZOwdOlSODg4wMrKChMmTEBxcTEAiI4NAGfOnBGd2+nTpxt9nw4dOiT2Hn388cdSvb9NiZK6PCws2MQOAP99OAghpDn7+eef8eLFC7i7u8PJyQnp6enYvHmzxE62WVlZ+OSTT2BkZITi4mJs2rQJfn5+iI2NBY/HQ0ZGBkJDQ/Hvv/8CAJ4/f4433ngDsbGxKCkpQYcOHZCRkYElS5Zg2rRpdcrfv38/NmzYABMTExQUFGDv3r346quvAAC+vr6wsbEBAFhYWIguLVhaWjZ6jt26dYOpqSk6deqEr776ChUVFYq8ZSpBSV1egYHsv/v3Axr4hyWEEGk9fvy4zvViWZ05cwZ5eXlITk7GgwcP8PnnnwMADh8+jPLycrFtX716hZMnT+LevXuiJuzU1FRcv34dd+7cgbGxMQQCAc6cOQMAiI6ORkZGBuzs7PDgwQNcv34dBw4cAMC2JKSmpoqVb2xsjH///Repqanw8/MDACQmJgJga9vCywldu3YV1bq7du3a4PlZWVmhdevWMDIywp07d7Bo0SJMmjRJ5vdJ1Sipy+vttwFHR3bmtt9+U3c0hBAiN0NDQ1GNVd5OcX/88Qe8vLxgYmICDoeDL7/8EgBQVVWFZ8+eiW1rZWWF3r17g8vlwtnZGQDg5eUFV1dXmJubw9bWFgA7+BUAXPrvMmdOTg5sbW3B4XAwevRoAADDMLh48aJY+QMGDICjoyO4XC48PT3FypLHhQsX8Pz5cyQnJ+Pp06cYMGAAAODHH39ERkaG3OWqAnWUk5eeHjBhArBmDdsEL2yOJ4QQLSO8pi6vvXv3Yt68eaKynJyckJeXh4cPHwIAqqurxbav2dStr69fZ52wpYBhGLH9LCws0LFjxzrHNzU1FXveokWLOuXXLksWNX/kmJqaYsyYMfjzzz8BABkZGRrVYY5q6ooQNr0cPQrU+iVKCCHNhfAHgYWFBdLS0nDx4kUMHjxYaeV3794dAJug9+/fL2oy//333/HJJ59gzJgxMpUn/BFQKsWMm2fPnsWBAwdEP0zKy8vx888/i153cXGR6diqRkldER07At26AVVV7LV1Qghphrp06QIAKC4uhpubG9zc3PDjjz8qrfzQ0FA4OjrixYsX8PDwgI+PD9zd3dGyZUsECvs3yUDYJH/lyhV07twZb7zxBsrKyiRu+/DhQ/j7+4PH46FLly7g8/n4448/AADBwcFwdHSU/8RUgJK6ooS1deoFLxe61Y0Q7RcSEoKIiAjY2NiguLgY/fv3R2RkpNLKb9WqFS5cuIDg4GC0bNkSt2/fRllZGfr27YsNGzbIXN7kyZMxduxY8Hg83Lp1CxcvXqxziUCoT58+mD59OpydnZGWlgaBQAA/Pz9s27YN27dvV/TUlI7DKHKhQYsUFRWBx+OhsLBQqlsXpPbsGcDns7X1W7eATp2UV7aaVVVV4fF/Q+G6uLiIrk3V3ubJkydo3bq1xNclbcflckUD0QgEAqn2J4SQ5kKRfEU1dUW1agW89x77+Pvv1RuLlhAm+JqDRhBCCFEcJXVlEF7T+f57oJ4mHPKaQCBASUkJNbkTQoiSUVJXhmHDgJYtgcxM4MQJdUdDCCGkmaKkrgxGRuw96wAQE6PeWDQU1c4JIUT1KKkri3B8419+ARQYuUhXUVInhBDVo6SuLJ07Az16sL3gqcMcIYQQNaCkrkzC2npMDNA87hQkhBCiQSipK9MHH7DzrN+9C/zzj7qjUToaKIYQQjQbJXVlsrQE/vc/9rGOdZgTCATIz89Hfn4+JXVCCNFQGpHUV69eje7du8PCwgK2trYYPXo0UlJSxLbp379/nfl+p0+frqaIGzBlCvtvfDxQVKTeWJRIIBCgsLBQrpmOhPvSjwFCtMOWLVvg6uoKY2Nj9OzZUzT1qSSurq51vps5HA5CQ0Mlbv/VV1+Bw+EgPDxc7bFLu8+yZcvqnJ9w/HhNoxFJ/cyZMwgNDRXNulNZWYnBgwfXmUFn6tSpyMrKEi1r165VU8QNePNNwMMDePmSTexEqh8E1LRPiGaIj49HREQEli5diqtXr8Lb2xtDhgxBbm6uxO0vX74s9r38+++/AwD8/f0lbvvtt9+KJoCRRv/+/REXF6eS2GXZp1OnTmLn+ffff0t9Dk1JI5L68ePHERQUhE6dOsHb2xtxcXFIT09HUlKS2Hampqawt7cXLUodw11ZOJzXtfXvvlNvLFqEkjrRdZcuXUL//v1hYmICT09PXLlyBdu3b8fIkSPVHZqY9evXY+rUqQgODkbHjh2xbds2mJqaYufOnRK3b9Wqldj38pEjR+Du7o5+/fqJbVdSUoKPPvoIO3bsgJWVlUbELss++vr6YudpY2OjknNQlEYk9doKCwsBANbW1mLr9+7dCxsbG3h5eWHRokV4+fJlvWVUVFSgqKhIbGkykyYB+vrApUvAzZtNd1xCiEa6cOEC+vXrh2HDhuHGjRvo0KEDIiMjsWbNGixfvlxs21WrVsHc3LzBJT09XSVxvnr1CklJSRg0aJBoHZfLxaBBg3D+/Hmp9t+zZw8mT54MDocj9lpoaCiGDRsmVrYyyRO7LPvcv38ffD4fbm5u+Oijj1T2N6ipdg6rqKhodB+NmxZLIBAgPDwcvXv3hpeXl2j9hx9+CBcXF/D5fNy4cQMLFixASkoKDh48KLGc1atX1/nP0mRsbYGRI4GDB9kOcxs3qicOQohGiIiIgL+/Pz799FMAwPjx4zF+/HiMGjUKvr6+YttOnz4d/xN2uK0Hn89XSZx5eXmorq6GnZ2d2Ho7OzvcvXu30f0PHz6MgoICBAUFia3fv38/rl69isuXLzdaxqpVq7Bq1SrR87KyMly4cAEzZ84Urbtz5w6cnZ0Vjl3afXr27Im4uDh4eHggKysLy5cvR9++fXHr1i1YWFg0ek7ycnJyEnu+dOlSLFu2rMF9NC6ph4aG4tatW3WuV0ybNk30uHPnznBwcMDAgQPx4MEDuLu71yln0aJFiIiIED0vKiqq8wapVEgIm9S//x746ivA2Ljpjk0I0RhPnjzB+fPn8fXXX4vW6evrg2EYiRUPa2vrOq2Uslq4cCHWrFnT4Db//vuv0jt7xcTEYOjQoWI/OjIyMjB79mz8/vvvMJbie7D2j5qPPvoIY8eOxfvvvy9ap6ofNfUZOnSo6HGXLl3Qs2dPuLi44Mcff0SIcHwSFcjIyBC7zGxkZNToPhqV1GfOnIkjR47g7NmzaN26dYPb9uzZEwCQmpoqMakbGRlJ9QaozJAhgJMTkJEBJCQAEyeqLxYtJBAIRHOuc7kaeZWIEKn8+++/AICuXbuK1qWkpKBHjx7o3Llzne1r11QlkVRTrWnu3Ll1asu1ubm51VlnY2MDPT095NQa6jonJwf29vYNlvf48WP88ccfdVpPk5KSkJubK3b+1dXVOHv2LKKjo1FRUQE9PT3Ra7V/1JiYmMDW1hZt27Zt8PjyxC7v+bZo0QLt27dHampqgzEpytLSUua+YxqR1BmGwaxZs3Do0CGcPn0abdq0aXSf5ORkAICDg4OKo5OTnh4wbRrwxRfAtm2U1GUk7Dhnbm5OSZ1otcLCQujp6YmuMefn5+Prr7+Gt7e3xO2V0fzeqlUrtGrVSuZYDQ0N4efnh8TERIwePRoA+38xMTFRrPlbktjYWNja2mLYsGFi6wcOHIibtfoWBQcHw9PTEwsWLBBL6IqQJ3Z5z7ekpAQPHjzARE38Xmc0wIwZMxgej8ecPn2aycrKEi0vX75kGIZhUlNTmcjISObKlStMWloa8/PPPzNubm7MW2+9JfUxCgsLGQBMYWGhqk6jrsxMhtHXZxiAYa5fb7rjKkllZSWTmprKpKamMmVlZaLHlZWVYtukpaWJrZNUTkpKCvP3338zZWVlTFlZGXP+/HmmrKxMtH9ZWZlYOdKUS4g2uH//PgOAiYyMZP79919m8ODBTNeuXRk7Ozvm0aNH6g6vjv379zNGRkZMXFwcc+fOHWbatGlMixYtmOzsbNE2mzdvZgYMGCB6Xl1dzTg7OzMLFiyQ6hj9+vVjZs+eLfG14uJisTwgaamqqpIr9tpxS3u+c+fOZU6fPs2kpaUx586dYwYNGsTY2Ngwubm5Up2vrBTJVxqR1AFIXGJjYxmGYZj09HTmrbfeYqytrRkjIyOmbdu2zKeffirTCaslqTMMw4wbxyb1GTOa9rhKQEmdEOWIjIxkWrZsyRgbGzNBQUFMXl4e07VrV8bT01PdoUm0efNmxtnZmTE0NGR69OjBXLhwQez1pUuXMi4uLqLnJ06cYAAwKSkpUpXfUFJfunRpvTlBuKSlpckVe+24pT3fgIAAxsHBgTE0NGQcHR2ZgIAAJjU1VapzlYci+YrDMM1j5pGioiLweDwUFhY27f3tf/4JDBwImJsDmZmACntKKltVVRUeP34MAHB0dMTTp08BAC4uLtDXZ6/cpKenY8OGDZgzZ0691/iqqqrw8OFDPHv2DH5+fgDYyyc+Pj7Q19fHkydPYG9vj+zsbLRu3Rr6+vqoqqrCkydPRM8JIaS5UCRf0cVKVXv7baB9e6CkBNi3T93RKF1MTAzu3r2LGB0b654QQrSRQkm9srISGRkZSElJQX5+vrJi0i0cDiAco37rVp2bkjUkJASenp4qva2DEEKIdGRO6sXFxdi6dSv69esHS0tLuLq6okOHDmjVqhVcXFwwdepUqQYYaFYCA9n71K9fBy5eVHc0SsXn8zF79uwmv2+UEEJIXTIl9fXr18PV1RWxsbEYNGgQDh8+jOTkZNy7dw/nz5/H0qVLUVVVhcGDB+Pdd9/F/fv3VRW3drG2BgIC2Mfbtqk3FkIIITpLph5Ily9fxtmzZ9GpUyeJr/fo0QOTJ0/Gtm3bEBsbi7/++gvt2rVTSqBab/p0YNcudua29evZRE8IIYQokUw19R9++KHehF6TkZERpk+fjsmTJ8sdmM7p2RPw8QHKywEppxEkhBBZXbt2DQYGBujfv79Kyt+6dSu6dOkiGu2sV69eOHbsmMzb1NZUc5bLE9vZs2cxYsQI8Pl8cDgcHD58WOlxKYtSer/n5eXht99+wy+//IKsrCxlFKl7OBxgxgz28ZYtQHW1euMhhOiksLAwzJs3D9evX1dJ+a1bt8ZXX32FpKQkXLlyBQMGDMCoUaNw+/ZtmbaRRN45y2WZc12e2EpLS+Ht7Y0tW7ZIdQx1UvgG4J9++gkhISFo3749KisrkZKSgi1btiA4OFgZ8emWjz4CFi4EHj4Ejh0Dhg9Xd0SEEB2yb98+WFlZITQ0FF999RUePXoEV1dXpR5jxIgRYs+//PJLbN26FRcuXBC15EqzjSTCOctVSZ7Yhg4dKjapiyaTuaZeUlIi9nz58uW4dOkSLl26hGvXriEhIQGff/650gLUKWZm7OxtALB5s3pjaSaEY8gLBAJ1h0KISpWWluKzzz7DmjVr0Lp1a/B4PNEcGbUpa8726upq7N+/H6WlpejVq5fc2wg19ZzlssSmLWRO6n5+fvj5559Fz/X19ZGbmyt6npOTA0NDQ+VEp4s++YRtij95EpBifmJ10KVEqEvnQkhDVq1ahXfffRcdOnQAAHTs2LHepD59+nQkJyc3uDR0m+rNmzdhbm4u6j916NAhdOzYUeZtahLOWX78+HFs3boVaWlp6Nu3L4qLiyWea80fIH/99RemT58u9Y8SWWPTJjI3v584cQKhoaGIi4vDli1bsGnTJgQEBKC6uhpVVVXgcrlSX9toltq0AUaOBH7+GYiOZhcNI0yE0sx9TAhRv4cPH+Lbb7/FrVu3ROu8vLzqTeqKztnu4eGB5ORkFBYW4sCBAwgMDMSZM2fEEqM029Qky5zlis65Lmts2kTmpO7q6orffvsNP/zwA/r164ewsDCkpqYiNTUV1dXV8PT0pGTQmFmz2KQeFwd8+SXA46k7IkKIFpszZw6eP3+O1q1bi9YJBIJ652NQdM52Q0ND0fzmfn5+uHz5MjZt2oRvv/1Wpm0a0tCc5fLOua6s2DSZ3B3lxo8fj6FDh2LevHno378/tm/fDh8fHyWGpsMGDAA6dgTu3GET++zZ6o6IEKKlTp48iXPnzuHatWtikx9dvnwZkydPRkFBAVq0aCG2jzLmbK9JIBCgoqJC4W1qaso5y2WNTZPJldSPHj2Kf//9F97e3vjuu+9w5swZfPTRRxg6dCgiIyNhYmKi7Dh1C4fD1tZnzGCb32fNAriaObcOh8OBpaUluBoaHyHNWWVlJcLDw/Hpp5/WqVQJZ/dKTk6uc8+6Is3vixYtwtChQ+Hs7Izi4mLs27cPp0+fxokTJ6TeJjo6GocOHUJiYqJon3nz5mHEiBFwcXFBZmYmli5dCj09PYwfP75ODCUlJWKdtvfv3w8AyM7OFq1r1aoV9PT05Iq/dnwlJSViLQZpaWlITk6GtbV1va0Z6iLzN/XcuXMRHByMy5cv4+OPP8aKFSvQr18/XL16FcbGxvD19W30Rn4CYMIEttk9NRU4flzd0dSLy+WiRYsWlNQJ0UDR0dF4/vw5Zs6cWec1JycnmJqa1ntdXV65ubmYNGkSPDw8MHDgQFy+fBknTpzAO++8I/U2eXl5ePDggVi5T548wfjx4+Hh4YH//e9/aNmyJS5cuIBWrVrVieHrr7+Gg4NDg0tGRobc8deO78qVK/D19YWvry8AICIiAr6+vliyZIn8b6SKyDyfesuWLXHy5En4+fkhPz8fb7zxBu7duyd6/c6dO/j444/x119/KT1YRahtPvWGzJ3LDhk7eDBQ41eiugnnMq85xzkA0dzqNedTl2bec3XOp07zshNCtE2TzqduZmaGtLQ0AEBGRkadTnEdO3bUuISusUJDX9/e1shIS7pKIBCgpKSEbjkjhBAlkDmpr169GpMmTQKfz0e/fv2wYsUKVcTVPLi5AaNHs483blRnJCqXk5ODnTt3IjMzU2w9JXVCCFEemZP6Rx99hIyMDPz888949OgRRo0apYq4mo+5c9l/v/8eqDGIj65JSEjAgwcPEBsbi8zMTMTExNRJ8I2hgWQIIaRhMiV14Qg9LVu2RPfu3evcJlHb06dP5Q6s2XjzTaBHD6CiAvjmG3VHozL+/v5wd3dHcHAwYmNjkZqaitjYWJnKoKROCCENkympd+/eHR9//DEuX75c7zaFhYXYsWMHvLy88NNPPykcoM7jcICICPbxN98AZWXqjaeGzMxMbNq0SeYadU3CRNyqVStMnjwZfD4fwcHBaNu2LU36QwghSiZTd+A7d+7gyy+/xDvvvANjY2P4+fmBz+fD2NgYL168wJ07d3D79m107doVa9euxXvvvaequHXL2LGAszOQng7s3QtMmaLSwwkEAhQVFTV6/3lsbCzu3r2LmJgYfPHFF3Ifq7CwEDVvsuDz+QgJCZFpcIv6ypbmPAghpLmQ6ZuwZcuWWL9+PbKyshAdHY127dohLy8P9+/fB8Beb09KSsL58+cpoctCX//1qHLr1wOy3WUoM2mbsYODg+Hp6Vln3GVNQc3xhBAiTq4bd01MTDBu3DiMGzdO2fE0XyEhwLJlwL//svesv/uuuiMCn8/H7NmzG61RC5OrtbW1zDVmgUCAsrIySsyE6LiMjAxMnDgRubm50NfXxxdffAF/f391h6VzqM1SU/B4r5vd161TbywSZGZmIioqCjk5OXVeYxgGhYWFciVmSuqENA/6+vrYuHEj7ty5g5MnTyI8PBylpaXqDkvnUFLXJGFhgJ4e8McfwNWr6o5GTExMDFJSUhAfH6/uUAghWsjBwUE0Pr29vT1sbGyQn5+v3qB0ECV1TeLqCghnTlq7tskOK8216ZCQEHh4eCAgIECu8ktLS6k2TggBACQlJaG6uhpOTk5KLffs2bMYMWIE+Hw+OBwODh8+3Og+xcXFCA8Ph4uLC0xMTPDmm29KvMPr6dOnmDBhAlq2bAkTExN07twZV65cUWr8ykBJXdPMn8/+m5AA1JrwQFWkSep8Ph9hYWGws7OTuXyGYVBaWgoZpxkghOig/Px8TJo0Cdu3b1d62aWlpfD29saWLVuk3mfKlCn4/fff8f333+PmzZsYPHgwBg0aJDbOyosXL9C7d28YGBjg2LFjuHPnDtatWwcrKyuln4Oi5Erqjx8/xq+//opLly4pOx7i4wMMGQIIBBp5bZ0QQiTx9/dHq1atxJL1xYsXYWhoiJMnTwIAKioqMHr0aCxcuBBvvvmm0mMYOnQoVq5ciTFjxki1fVlZGX766SesXbsWb731Ftq2bYtly5ahbdu22Lp1q2i7NWvWwMnJCbGxsejRowfatGmDwYMHw93dXennoCiZk/oPP/yA9u3bY9SoUejVqxe6deuGZ8+eqSK25mvBAvbf2FidHjpWEXQ7GyGaJSoqCmPHjkVkZCQAdg7yCRMmYMaMGRg8eDAYhkFQUBAGDBiAiRMn1lvOqlWrYG5u3uAiHN1UUVVVVaiurq4zMZmJiQn+/vtv0fNffvkF3bp1g7+/P2xtbeHr64sdO3YoJQZlkzmpL1++HB9++CHu3r0r+vW1cOFCpQfWrPXvD3TvDpSXA1FR6o5GI0ka1IYQoj4ODg4IDw/H06dP8fz5c4SFhcHIyAhr1qwBAJw7dw7x8fE4fPgwfHx84OPjg5s3b9YpZ/r06UhOTm5wUXTgKiELCwv06tULK1asQGZmJqqrq7Fnzx6cP38eWVlZou0ePnyIrVu3ol27djhx4gRmzJiBsLAw7Nq1SylxKJPM96k/fPgQx48fh6urK9q3b489e/bAz88PMTExqoiveeJw2Nr6uHHAli3sYwsLdUdFCCENat++PUxNTbFkyRLs3bsXly5dEtWC+/TpI1XLmrW1NaytrVUdqsj333+PyZMnw9HREXp6eujatSvGjx+PpKQk0TYCgQDdunXDqlWrAAC+vr64desWtm3bhsDAwCaLVRoy19Srqqpgamoqeu7p6QmBQIDs7GylBtbsjR4NtGsHFBQAGtrMoyzCWjc1pROi3bhcLjp37oxvvvkGK1euhLe3t8xlNGXzOwC4u7vjzJkzKCkpQUZGBi5duoTKykq4ubmJtnFwcEDHjh3F9uvQoYNS41AWuTrK7dq1C//88w9KSkoAsIMKvHz5UqmBNXt6esCnn7KP161jZ3HTUQzDoLi4WCVN6cqYlIYQIh3h/+GuXbtirnBaaQmqq6vrfa0pm99rMjMzg4ODA168eIETJ06ITSveu3dvpKSkiG1/7949uLi4KD0ORcmc1Pv27YuVK1eiT58+aNGiBdq1a4fy8nLExMTg1KlTKC4uVkWczdOkSQCfD2RmAhp47UaVlJWMY2JiRJPSEEJUa+PGjbh48SIEAkGdIaNHjhyJTz75BN27d29w2mVra2u0bdu2wUVfX/KV45KSElHiB4C0tDQkJyeLatTR0dEYOHCg2D4nTpzA8ePHkZaWht9//x1vv/02PD09xWaRnDNnDi5cuIBVq1YhNTUV+/btw/bt2xEaGirP26RSMif1M2fOoLCwECkpKdizZw/GjBmDfv36YevWrRg4cCCsrKzQoUMHVcTa/BgZvb5vffVqoLJSvfE0IWUl45CQEI2elIYQXXHz5k0sWrQIn3zyCe7cuYOqqqo6r3t4eODy5cuYoqKZKK9cuQJfX1/4+voCACIiIuDr64slS5YAAPLy8vCg1vgfhYWFCA0NhaenJyZNmoQ+ffrgxIkTMDAwEG3TvXt3HDp0CD/88AO8vLywYsUKbNy4ER999JFKzkMRHEaJbZ5paWm4cuUKrl27JupQoCmKiorA4/FQWFgIS0tLdYcjvZcv2ZHmnj1ja+uTJilcZFVVFZ48eYLWrVtDX1+/3uf29vbIzs5G69atAbDjEwCAi4uL6JdyVVWVxPU1j/Xw4UNkZWWBw+GgR48eAIDk5GT4+PigqqoKp06dQps2bWBiYiIqIz09HRs2bMCcOXPg7OxcJ+764mnsXAkhyldeXo7u3buja9eu2LhxI6ytrXHz5k14eXkBYEdt69Spk0Zeg9ZEiuQrpY4o16ZNG/j7+2tcQtdqpqaA8NrUqlVAA9eiNFVOTg7i4uKQl5cn9T7SzhBHCFG/hQsXorS0FNHR0bCysoKLiws2btwounx2+/ZtlQw2Q+rSiGFiV69eje7du8PCwgK2trYYPXp0nU4J5eXlCA0NRcuWLWFubo6xY8dKnDFMJ82YAVhZASkpwE8/qTsaMVwuFzweDxwOp95tEhIS8PDhQxw5cqTOa1lZWfjhhx+QS4PsEKKVTp48iS1btmDPnj2w+O/W28WLF+Pw4cOia843b95E586d1Rlms6ERSf3MmTMIDQ3FhQsX8Pvvv6OyshKDBw8Wm5Zvzpw5+PXXX5GQkIAzZ84gMzMT77//vhqjbkKWlsDs2ezjL78ENGjAFS6XixYtWjQ4j7q/vz/c3NwwfPjwOq/FxcXh0aNH+EnDfqwQQqQzePBgVFZWitXEp0yZgry8PBw6dAgAJfWmpBEXGY8fPy72PC4uDra2tkhKSsJbb72FwsJCxMTEYN++fRgwYAAAIDY2Fh06dMCFCxfwxhtvqCPspjVrFvD118CNG8CvvwIjR6o7IqnZ2dlh+PDhOHLkiGjqxZiYGCxatAhBQUFIS0vD2LFj1RskIURlomhkzCajETX12goLCwFANKpQUlISKisrMWjQINE2np6ecHZ2xvnz5yWWUVFRgaKiIrFFq1lbAzNnso+XLdOo2ro0jhw5ggcPHiA2NhaxsbFITU1FbGwsHBwcMH78eADsf3y6n5wQQli1c1iFFOOVKFRTT0xMRGJiInJzc+uMBrZz5065yhQIBAgPD0fv3r1FPSezs7NhaGiIFi1aiG1rZ2dX70h2q1evxvLly+WKQWPNnQtERwPXrgGHDgFadPlh+PDh+O2330T3fmZnZ4vdB/rTTz/hwYMHiImJ0b2/GyGEyKH2fPNLly7FsmXLGtxH7pr68uXLMXjwYCQmJiIvLw8vXrwQW+QVGhqKW7duYf/+/XKXAQCLFi1CYWGhaMnIyFCoPI1gYwOEh7OPly5lp2fVEjY2NggKCgKfzwefz0dISIhYz/axY8fCw8OD7icnhJD/ZGRkiOWxRYsWNbqP3DX1bdu2IS4ursEp9GQ1c+ZMHDlyBGfPnhXdiwwA9vb2ePXqFQoKCsRq6zk5ObC3t5dYlpGREYyMjJQWm8aIiAA2bwZu3QISEoCAAHVHpBS2trYICwvTyFvYBAIBioqKYGlp2WCHQEIIUSZLS8umu0/91atXSrvvkGEYzJw5E4cOHcKff/6JNm3aiL3u5+cHAwMDJCYmitalpKQgPT0dvXr1UkoMWsPKik3sAHttvQnvW8/JyWmW171p7nZCiLaQO6lPmTIF+/btU0oQoaGh2LNnD/bt2wcLCwtkZ2cjOzsbZWVlAAAej4eQkBBERETg1KlTSEpKQnBwMHr16tU8er7XFh7OJve7d4Effmiyw8bHxyMlJYXGUSeEEA0ld/N7eXk5tm/fjj/++ANdunQRGycXANavXy91WVu3bgUA9O/fX2x9bGwsgoKCAAAbNmwAl8vF2LFjUVFRgSFDhuCbb76RN3ztZmnJzuD22WfA8uXABx8ATTAEakBAAOLj4+m6NyGEaCi5M8GNGzdE9xzfunVL7LWGRheTRJrh542NjbFlyxZs2bJFprJ11qxZwIYNQGoqsHs3MHmyyg9pZ2ensde9CSGEKJDUT506pcw4iKzMzYEFC4B584AVK4AJEwBDQ3VHJTPhMLFTpkzRyLmJCSFEm1BXXm02YwZgbw88egQ0MD+xJqNhYgkhRHkUSuoFBQVYt24dpkyZgilTpmD9+vWi0eBIEzA1BYT3La5cCZSXqzceOQwdOhTFxcXo0aMH9S4nhBAFyZ3Ur1y5And3d2zYsAH5+fnIz8/Hhg0b4O7ujqtXryozRtKQadMAR0fgyRPgu+/UHY3Mjh07BnNzc5w6dUqqvhWEEELqJ3dSnzNnDkaOHIlHjx7h4MGDOHjwINLS0jB8+HCEC0c9I6pnbAx8/jn7eNUq4OVL9cYjo6CgILi6uuK9995TdyiEEKL1FKqpL1iwAPo1bqXS19fH/PnzceXKFaUER6QUEgK4uABZWezY8FrEwcEBH3zwAWxsbNQdCiGEaD25k7qlpSXS09PrrM/IyICFhYVCQREZGRoCkZHs41WrgOfP1RsPIYQQtZA7qQcEBCAkJATx8fHIyMhARkYG9u/fjylTpoim0iRN6KOPAG9voLAQ+PJLhYrKzMzEpk2bZB4OloZTJYQQ9ZI7qX/99dd4//33MWnSJLi6usLV1RVBQUEYN24c1qxZo8wYiTT09ADh+75lC5CWVu+mAoEAJSUl9SbfmJgY3L17F7Ey3iYnKalnZmYiOjoaeXl5MpVFCCFEdnIndUNDQ2zatAkvXrxAcnIykpOTRT3gdXJ2NG0weDAwaBDw6hWweHG9mzWW1ENCQuDp6Sk237m8YmJikJKSgiNHjihcFiGEkIYpPPiMqakpOnfujM6dO8PU1FQZMRF5cTjA2rXs4337ADlvLeTz+Zg9e7ZShoMNCQmBh4cHhg8frnBZ8l4WIISQ5kKmYWIjIiKwYsUKmJmZIUI4/Wc9ZJnQhSiRry87ZOyePeykL3/8wSZ7NeHz+Zg5cyaysrJE67hcLszNzWWem1x4WSAmJgZffPGFVPtwuVy0aNGC5kEnhDQLMiX1a9euobKyUvS4PrJO6EKUbOVK4McfgT//BE6cAN59V90RiZE3qYeEhKCoqEimWeKESZ0QQpoDmZJ6zUlcdu3ahdatW9f5YmYYBhkZGcqJjsjHxYWdxW3dOnbSl3feYTvSKYDL5YLH46GoqEhJQcpOmZcFCCFEF8ndJtmmTRuJPZrz8/PRpk0bhYIiSvDZZ0CLFsCNG2xTvIKoGZsQQjSf3N/Q9Y3TXVJSAmNjY7kDIkpibf16+NjFi4GyMvXGQwghROVknk9d2EGOw+FgyZIlYj3eq6urcfHiRfj4+CgtQKKAmTOBzZuB9HRg48bXM7oRQgjRSTIndWEHOYZhcPPmTRgaGopeMzQ0hLe3N+bNm6e8CIn8jI3ZYWMnTGBHmQsMBNR0PZrD4cjVOY4QQoj0ZE7qws5ywcHB2LRpEywtLZUeFFGiDz9kR5g7f56tqe/apZYwVJHUc3JyEB8fj7lz58LZ2Vlp5RJCiLaS+xs2NjaWEro24HCAqCj28e7dwMWL6o1HieLj45GSkoKYmBh1h0IIIRpB7qS+evVq7Ny5s876nTt30tjvmqZbN0A45OusWYCOTLgSEBAADw8Pme5bJ4QQXSZ3Uv/222/h6elZZ32nTp2wbds2hYIiKrBqFWBhAVy+DO6+faLV2jyzmp2dHcLCwui+dUII+Y/cST07OxsODg511rdq1UpsSFCiIeztgf+GVjX44gtwS0sBsEm9sLCw3lsUpaHomOz5+fnYvXs3cnJy5I5BlWjMeUKItpA7qTs5OeHcuXN11p87d45qTppq9mygXTtkZWcjPiJCaUmqoalahTPCNfSj4eTJk3j48CESEhKkPmZmZiaioqKa5IdAzTHnCSFEk8md1KdOnYrw8HDExsbi8ePHePz4MXbu3Ik5c+Zg6tSpyoyRKIuhIbBhA3YASLt5E7Hr1gFge5FHRUXJneQbmqpVmqQ+ePBguLm5wd/fX2x9dXU10tPTUVVVVWcf4ZSu8fHxcsUsC+H50bV7Qoimkzupf/rppwgJCcEnn3wCNzc3tGnTBrNmzUJYWBgWLlyozBiJMr33Hia/9Ra8GQbTUlIAKN6LXNEx2a2trTFp0iTY2dmJ1gknfcnOzpaY1IVTugYEBDRavqL9BmjMeUKItpA7qXM4HKxZswbPnj3DhQsXcOPGDeTn52PJkiU0S5sm43BgGxWFL/T04JSYCM6xY2rvRS7pmrpwApn67mvn8/kICwsT+yFQH23uDEgIIbJQeCSQ9PR0PH/+HA8fPsSJEyfwyy+/4JdfflFGbERFGA8PZP9Xw+XOng17S0u19iKX55o6IYSQumQeUU7o4cOHGDNmDG7evAkOhyO6ZiqspVdXVysnQqIST0JC4HD2LDhpabDauhX5/43prwzCmrG0gxMNHjwYZ86cqXNNnRBCiGzkrqnPnj0bbdq0QW5uLkxNTXHr1i2cPXsW3bp1w+nTp5UYIlEFgakpKr/+GgDQYscOGKSmKq/sWs3deXl5iIuLq7enuqRr6oQQQmQnd1I/f/48IiMjYWNjAy6XCz09PfTp0werV69GWFiYMmMkKiIYORKCYcPAqaxEq6VLASnvVZf1GvWRI0eoeZ0QQpqA3Em9uroaFhYWAAAbGxvR7VAuLi5I+a9XNdFwHA4EGzdCYGwMk4sXwakx0lxDZE3qw4cPl3jLGiGEEOWSO6l7eXnh+vXrAICePXti7dq1OHfuHCIjI+Hm5qa0AImKubriRWgoAID76afAixdKP4SNjQ2CgoIkNq9zuVyYmJjQHROEEKIEcif1xYsXi2pqkZGRSEtLQ9++fXH06FFECWcFI1qhICQEFe7u4Dx7BmbRoiY9NiV1QghRHrl7vw8ZMkT0uG3btrh79y7y8/NhZWVFX9DaxtAQucuWwWniRGD7dmDCBKB1a3C5XLRo0UKpc6DXxuVyYWZmhuLiYpUdgxBCmgu5vq0rKysxcOBA3L9/X2y9tbU1JXQtVdajB4rHjAGHYaA3fTrw6lWTJXVzc3OVHqMxNDgNIURXyPVNamBggBs3big7FqJmeYsWgbGxAef2bfA0bPpcVSZeSuqEEF0hd/VowoQJNGuVjhFYW0OwcSMAoEV0NHD7tnoDqkEZU8QSQoiukzupV1VVYevWrejWrRs+/vhjREREiC2yOHv2LEaMGAE+nw8Oh4PDhw+LvR4UFAQOhyO2vPvuu/KGThrA/O9/EAwfDk5lJfSmTQOaYGRAgUCAsrIyqikTQoiC5O4od+vWLXTt2hUAcO/ePbHXZL2uXlpaCm9vb0yePBnvv/++xG3effddsfm6jYyMZIyYSIXDgWDzZuDMGXAvXQI2bwbCwxUsktNgD3eGYVBeXk61cEIIUZDcSf3UqVNKC2Lo0KEYOnRog9sYGRnB3t5eacckDWjdGvkLF8Lm88+Bzz8HRo0C2rSRuzhN6AxHCCHNgczfsg8fPlRLjer06dOwtbWFh4cHZsyYgefPnze4fUVFBYqKisQWIr2SDz6AoF8/4OVLYNo0qYeQJYQQohy1c1hFRUWj+8ic1Nu1a4dnz56JngcEBNQ7UYeyvPvuu9i9ezcSExOxZs0anDlzBkOHDm1wJrjVq1eDx+OJFicnJ5XGqHO4XAi2bQOMjYE//gB27mx0l8zMTGzatEk0ZHBDqMc5IYQ0zMnJSSyPrV69utF9ZE7qtWvpR48eRWlpqazFyOSDDz7AyJEj0blzZ4wePRpHjhzB5cuXG5wNbtGiRSgsLBQtGRkZKo1RJ7VtC6xcyT6eMwd49KjBzWNjY3H37l2xvg/1YRgGhYWFlNQJIaQeGRkZYnlskRQjfmrlRU43NzfY2NggtYHpQo2MjGBpaSm2EDmEhwO9ewPFxeBOnQo0kISDg4Ph6emJ4OBgmQ/D4XBgYWFB190JIeQ/tXOYNB3EZf4GFd5SVntdU3ry5AmeP38OBweHJj1us6SnB8TFAaam4J4+DYvdu+vdlM/nY/bs2eDz+XIdSpmj1ynSvE+XBggh2krm3u8MwyAoKEj0i6G8vBzTp0+HmZmZ2HYHDx6UusySkhKxWndaWhqSk5NhbW0Na2trLF++HGPHjoW9vT0ePHiA+fPno23btmLjzxMVatsW+L//A0JDYbVmDSrHjgX+m3ZXUwkTs7m5uUL7UssBIUSbyJzUAwMDxZ5PmDBB4SCuXLmCt99+W/RcOHhNYGAgtm7dihs3bmDXrl0oKCgAn8/H4MGDsWLFCrpXvSlNnw7BwYPgJibCYMoUYO9edUekEIFAgKKiIrosQwjRKTIndWk6Qcmqf//+Dd4md+LECaUfsznLzMxETEwMFi1aBGdnZ+l24nIh2LED8PYG99Il8HbsAL78Uord2HvUNW2iH0Vq8oQQoqmobbEZio2NRWpqquw/0Jyc8HzpUgBAiw0bgGvXGt2l5sAzsiZ24Y8PaW6RI4QQQkm9WQoODkbbtm3l6qVe+v77qB4xgh0bftIkdnAaKck6aJHcPz6UrCmmoCWEEGWgb6lmiM/nIyQkRL5e6hwOKr/5BlW2tuDcvQvMnStzEQKBACUlJRAIBOByubCwsJBYi1fkx4cyUVInhGgL+pZqxuS+dcvGBnnr1rGPt20D59dfZT5uSUkJGIYBl8uFtbU1zMzM6iTN2j8+OBwOeDweJVdCCKkHfTs2Y/LOUZ6ZmYk1SUl4MnUqAIA7bRr0cnOl3r925zkulwsTE5NGkzXVmAkhpGFyz9IGAImJiUhMTERubm6d2t5OKcYKJ9pJOBzsjq5dsdzHB5zkZNjMmwccP15n28zMTERFRWHcuHGiGrcwqTeEBoAhhBDZyV3lWb58OQYPHozExETk5eXhxYsXYgvRPcKackhICDw9PRHy8cfAvn1gTExg8tdf0Nu0qc4+MTExSElJQUJCgkzHEt5HritzrNOPFEJIU5C7pr5t2zbExcVh4sSJyoyHaID6mrmF683NzTFr1iyYmppC0Lo1mK+/hl5oKPSXLIGRhwfw3nuifUJCQlBQUIBx48Y19WloFBqljhDSFOT+dnn16hXefPNNZcZCNIQ0165rzrLGTJ2KkhEjwKmqQqtZs4Aac93z+XyEhYXBzs5OoZikndaVasSEkOZM7qQ+ZcoU7Nu3T5mxEG3F4eD5l19C0LYt9DMzYdDIbG7yiImJwd27dxETE9PgdpTUCSHNmdzN7+Xl5di+fTv++OMPdOnSBQYGBmKvr1+/XuHgiPZgLCxQuXcvDPv2hd6xY7Dcvv31XOxKEBISgqKiIoSEhCitTEII0TVy19Rv3LgBHx8fcLlc3Lp1C9euXRMtycnJSgyRqFJOTg6io6ORk5OjcFlMly54vnQpMgFErV2LzJ9/VjzA/ygyrau0TfcNoRYAQog2kLumfurUKWXGQdQkPj5e1Dt99uzZCpdXMn48Yvftw63bt7EzOBiL791TQpSKqdl0v3z5crnKoI5uhBBtQN9OzVxAQAA8PDzg7++vnAI5HAT+8AM6tWiBacXF4AYFNXp9veakLzXl5ORg586dCk/oIroFT8qme2XU7AkhRB0USuoFBQVYt24dpkyZgilTpmD9+vUoLCxUVmykCdjZ2WHmzJkK906vid+uHT7etw8OxsbgnjgBq2++aXD7+pJ6QkICHjx4oPCELrI23UvbKY8QQjSN3En9ypUrcHd3x4YNG5Cfn4/8/Hxs2LAB7u7uuHr1qjJjJFqoskMHCKKiAABWmzbB7PRpmcvw9/eHu7u7Sid0kXStXNaaPSGEaAq5k/qcOXMwcuRIPHr0CAcPHsTBgweRlpaG4cOHIzw8XIkhEm3FBAdDMG0aOAwD+3nzYPDokcTthBO8VFVViSVYOzs7TJ48Wb7Z5KQkafx7RTrlEUKIOilUU1+wYAH09V/3tdPX18f8+fNx5coVpQRHtJ9gwwaU+flBr7gY9tOnA8XFdbepJ6kTQgiRjdxJ3dLSEunp6XXWZ2RkwMLCQqGgSNPLyclBVFSUwp3D6jRnGxoiJzoaVba2MExNBTckBKg1nntmZiZiYmKQlZUl9XG4XC6srKxgZWVFvdEJIeQ/cn8bBgQEICQkBPHx8cjIyEBGRgb279+PKVOmYPz48cqMkTSBhIQEpKSkKNw5jGGYOhOxVLdqhczNm8EYGIB76BA4a9aI7RMbG4vU1FTExcVJfRzhPOzW1taU1Akh5D9y36f+9ddfg8PhYNKkSaiqqgIAGBgYYMaMGfjqq6+UFiBpGv7+/jhw4IDKOoc9dnBAdN++CP/zTzgsWYIiExPEpadj7ty5CA4ORnZ2NoKCglRybEIIaS7kruIYGhpi06ZNePHiBZKTk5GcnCzqAW9kZKTMGEkTsLOzQ1hYmMo6hyUkJOCGgQG2+/mBwzA4/NlnuJeUhJiYGPD5fISEhMDBwUHh4wivz9N1eUJIc6Rwu6WpqSk6d+6Mzp07w9TUVBkxER3k7+8PDw8PTI6PB9O7Nz5+9Qrdb95E8KhRSu0cR0mdENKcydT8HhERgRUrVsDMzAwRERENbksTujQf0kzVKmoJcHFBdUICWvXsiciMDAjCwvBg2zaxa/BNLScnB/Hx8Zg7dy6cnZ1VcgzhKHVz5sxR2TEIIUSmmvq1a9dQWVkpelzfQhO6aBcOhwNLS0twOBy59pcmqYtp1QpZ27ej2twc3HPnYLtkSZ0e8U1JOP69KkeQo1HqCCFNQaaaes1JXHbt2oXWrVvX+SJnGAYZGRnKiY40CWFSLi8vb7JjVrZrh5yoKDhMnQreoUNwtLUFfHya7PhCAoEAw4cPB8MwCA4ORkFBgUouI9HUsYSQpiD3NfU2bdogLy+vzvr8/Hy0adNGoaBI81D21luoWrcOAOD67bfQU+JUrdISCAQwNTXFzJkzYW9vr7LBb2iUOkJIU5A7qdd3DbSkpATGxsZyB0RUr74JVFRNODDNkydPEBUVhZycHFR//DFeTJwIADCeOhWG1641aUyNoXnUCSHaROb71IUd5DgcDpYsWSLWVFldXY2LFy/CRw3NqER66krqDMOgsLAQcXFxSElJQXx8PLp27YpnCxeCuX8f1hcuwG7KFDB//w20awcOh6OSOGVJ1MJt6YcqIUQbyPxtKewMxzAMbt68KdZB7u7du/D29pZpZDDS/AQHB8PDwwMBAQHsCn19pERGotrXF3r5+dAbPhzIyZHrx4dwyNnk5OR650SXNIkLIYToAplr6sLOcsHBwdi0aRMsLS2VHhTRbFwuFzweD0VFRXLtz+fzMXPmTLx8+VJUWxaYmqL8wAEY9u8Pg4cPoTdqFDg7d4IxM5OpbOGQs5999hnMzMwQExODL774os52wtvYZs6cKdc5EEKIJpJ7mNjY2FgAwJ07d5Ceno5Xr16JvT5y5EjFIiMaS9hbvqSkpM56Ho+HFy9eNFoGwzB1BolhbG2Rs2sXHP/3P3CSkmAXFobsbdvEjtlYrV045GxgYCB++umnenubC29ji4uLw3vvvddovIQQog3kTuppaWkYPXo0bt68CQ6HI2rKFN7rXF1drZwIidYQJt7CwkK5y6hq0wbVhw9D7513YHbmDFp98QWwb5+o7MYIh5z18vJC69at6+1tHhAQgPj4eAQFBSE3N1fieXC5XOogRwjRKnL3QAoLC0ObNm2Qm5sLU1NT3L59G2fPnkW3bt1w+vRpJYZItB2HwwGPx5N+cJuePfEkOhrLOByUHDgA7vz5Sh+cRjjCnaTx5mUeTIcQQjSE3N9a58+fR2RkJGxsbMDlcsHlctGnTx+sXr0aYWFhyoyRaDlhs7wsSfK7x49xpWNH7ADA3bgRghUraIx4QghphNxJvbq6GhYWFgAAGxsbUS9jFxcXpKSkKCc6olWUeU93SEgI3AcOxNjZswEA3KVLUb1hg0Jl14xP0jV9QgjRdnIndS8vL1y/fh0A0LNnT6xduxbnzp1DZGQk3NzclBYgUZ3MzEzRIDDKIO+tYjk5Odi5cyeysrJE6/h8PsLCwmAyaxaqlywBALSMjARn1y6F4qOBZAghukzupL548WLRl2NkZCTS0tLQt29fHD16FFFRUUoLkKhObGysaBAYVZHm+nRCQgIePHhQ7/gGzOLFEAhr7NOmAQcPqiJUQgjRenIn9SFDhuD9998HALRt2xZ3795FXl4ecnNzMWDAAJnKOnv2LEaMGAE+nw8Oh4PDhw+Lvc4wDJYsWQIHBweYmJhg0KBBuH//vryhk//UGQRGBsKpRCUN7lKTNEnd398f7u7uCAoKkrwBhwPB//0fiv/3P3AEAuCDD4CTJ+tsJutIcWVlZVRrJ4ToFKV277W2tpZr+s7S0lJ4e3tjy5YtEl9fu3YtoqKisG3bNly8eBFmZmYYMmRIk84qpouETdx2dnYy76vMqUTt7OwwefJkiT3RRTgcPF+1CoJx44DKSmD0aODPP8U2EQgEuH//PmJiYsSa8iVhGIaSOiFE58id1AcMGIDly5fXWf/ixQuZa+pDhw7FypUrMWbMmDqvMQyDjRs3YvHixRg1ahS6dOmC3bt3IzMzs06NnkhHGbdshYSEwNPTs8GpRKWtOUtda9bTg2DXLmDYMKCsDBg+vE5ib6wpX17StkwQQog6yf2tfvr0aURHR2P06NEoLS0VrX/16hXOnDmjlOAAdpCb7OxsDBo0SLSOx+OhZ8+eOH/+fL37VVRUoKioSGwhLGUkdWmmEq2Z1IW3tUlqyZGp1mxoCPz0E/Dee68T+39DFwNSNOXLKTY2VmktE4QQIo3aOayioqLRfRRqfv/jjz+QnZ2NN954A48ePVKkqHplZ2cDQJ0mYjs7O9FrkqxevRo8Hk+0ODk5qSQ+Ih2lDuhiZMR2lhMm9mHDgP8GPJKqKb+G7Oxs7Ny5s9E7AIKDgxttmWgIDWhDCJGVk5OTWB5bvXp1o/so9A3j4OCAM2fOoHPnzujevbtGjSS3aNEiFBYWipaMjAx1h0QakZWVJbGJu3YzvkAgQEFZGQQJCa8T+3vvgXPmjMwd4L7//ns8ePAACQkJDW4nTctEQ9SZ1OlWPkK0U0ZGhlgeW7RoUaP7yP0NI2xGNTIywr59+zB79my8++67+Oabb+QtUiJ7e3sAqFOTysnJEb0miZGRESwtLcUWono5OTmIjo6W69733bt3S2zirn3/uyhJCZvihw4FysrAHTkSJufP15vUJY0iN3HiRLi7u8Pf319nr5tTUidEO9XOYUZGRo3uI3dSrz3AyOLFi7F3716sW7dO3iIlatOmDezt7ZGYmChaV1RUhIsXL6JXr15KPRZRnHD2M3nufZ80aZLsTdzGxmxT/NCh4Lx8Cafp02F7+bLETSUldXt7e0yePBl2dnZ03ZwQovUUmqXNxsZGbN3YsWPh4eGBpKQkmcoqKSlBamqqWNnJycmwtraGs7MzwsPDsXLlSrRr1w5t2rTBF198AT6fj9GjR8sbPlGRgIAA7N+/H1OnTpW5qdnBwUGuJm6BoSGKYmNhOW0auL/8gu6rVqGsTRvg7bdlKic4OBglJSVyXzcnhBB1k7um7uLiIvFL28vLC4GBgTKVdeXKFfj6+sLX1xcAEBERAV9fXyz5b3jQ+fPnY9asWZg2bRq6d++OkpISHD9+HMbGxvKGT1TEzs4Os2fPhqenZ53Ph6Tm7aysLGzZsgXPnz+X+5jCa+yv9uxB4bBh4FZXw3TKFJjXc528vmMqet2cEELUTaaaekREBFasWAEzMzNEREQ0uO369eulLrd///4NjhfO4XAQGRmJyMhIqcskmkc4YE1sbCwmTpwIAIiLi8O9e/dQXFwsGqFQbgYGyPzqKxRUVcHlxAnYzJ+PagMDMDNnAnh9P3xsbCzu3buH/Px8dOvWDTt37kRQUBCcnZ0VPUVCCFErmZL6tWvXUFlZKXpcH3lGlSO6LyQkBEVFRQgODkZ1dTXMzc0xefJkFBcXY+DAgcrpGa6nhxuffAK7tm1hvGUL9ObMgaCwEC3+mw64rKwMEydORFFREbp16ybW+7179+6KH58QQtRIpqR+6r9BPiorK8HlcrFt2za0a9dOJYER3SNs3ra3t8fTp09hbm4OR0dHzJw5E8+ePZOqDIFAgKKiooZ7cnM4KF+xAuX6+mixaRO4y5ahRXY2Sr76CgB77T40NBSpqamYOHEi8vLy4O/vr4QzJIQQ9ZKramRgYIAbN24oOxZCGtXY9K4cDgcmJibg6umhIDwc1VFRAIcDbNsG4wkTwK01IlPN3u+1j0O3gRFCtI3c7Z0TJkygW390AIfDAY/Ha/JBUXJychAVFaX0e8K5XC6b1P87H2bGDCAhATAygv6RI+ixeDE4L140Wg4ldUKINpL7lraqqirs3LkTf/zxB/z8/GBmZib2uiwd5Yj6CEc6a2rC+9lrdppTmbFjgVatwIwcCes7d1A9dCj065kRsLbMzExs3rwZs2bNoiFeCSEaT+5vqVu3bqFr166wsLDAvXv3cO3aNdGSnJysxBBJU2jKYUy5XC6mTp0KT09PBAcH17udUlsR3noLZSdPoqxlS+ilpMD1ww9hmZbW6G7CAWliY2Mlvk41ekKIJpG7pn6qxsxYRPs1ZY2dy+XC09MT4eHhsLe3R3p6eqMx5eTkID4+HrNnz5b7uE+trPBl375YcfMmnB48QO+FC1FibY1yB4d6k7JwQJr6fnwIk7q5uTnV5AkhakffQkQr1Gyur09mZiaio6ORl5cn8fW4uDjcz83FNyNGINXXF1+WlaF00iTY798Ppp6kTgPSEEK0idw1daE7d+4gPT0dr169Els/cuRIRYsmWiYzMxNRUVEICAiAi4uLUssOCAhAfHw8goOD6+35Hhsbi+vXr+Pvv/+Gp6cnKisrxWrgQUFBePz4MSZNm4aNFRV48PQpkJuLpd99hxdFRcB33yk1Zk0hHMlvzpw5NMAOITpO7pr6w4cP4e3tDS8vLwwbNgyjR4/G6NGjMWbMGIwZM0aZMRItERMTI9dkLrm5udizZ0+DM7vZ2dkhLCyswRpzcHAwysvLYWRkhLi4uDqTtzg4OODDDz+Eg4MD/AMCYDVwID6IiADD4cDqxx9hMGYMuEVFdcqV5rq5Jl9bF47kR3erEKL75E7qs2fPRps2bZCbmwtTU1Pcvn0bZ8+eRbdu3TRqXnXSdEJCQuDh4YGAgACZ9ktISEBaWhpiY2PrJGJZ8Pl8rFu3Dj4+PggKCqrzelZWFvbt24esrCzRurKJE3Hl888hMDGBXmIiHMaMAe7eFduPYRgUFRWhurq63sStyUk9JCRE9tnvCCFaSe6kfv78eURGRsLGxgZcLhdcLhd9+vTB6tWrEfbfkJykeeHz+QgLCxMN5CLt/OQBAQHw8PBAYGCgQkkdYGvj8+bNg6OjY53X4uLikJaWht27dyMhIQGPHj3C3r17kdOzJx7t2QPG0REGDx9C7803wfnllzr7MwyDwsLCJk/civ5goH4BhDQfcif16upqWFhYAABsbGxEX9wuLi5ISUlRTnREq0nb7CtsWndwcJCqXHkHrhk6dCiKi4vx7rvvwt/fH66urvjoo48AABUdOqDin39Q3qMHOMXF0Bs7Fi3Wrwc0oOatya0AhBDNIndS9/LywvXr1wEAPXv2xNq1a3Hu3DlERkbCzc1NaQES7aWqZl9hT3hZrxEfO3YMFhYWOH78OGxtbTF+/HjY29u/3sDWFtl79kAwaxYAoMXmzTAYN07idXZCCNFEcif1xYsXi2oOkZGRSEtLQ9++fXH06FFERUUpLUCivRRp9m1oMBxhc73wx4K0zfxBQUFo06YNJk2a9HqMeC4XHA4HlpaW7LEMDCBYvx7VsbEQGBlB79gxOL7/Pgzv35f5HAghpKnJndSHDBkimv+6bdu2uHv3LvLy8pCbm4sBAwYoLUDSPDWU1Gv3hK85T3tDavZ+rzlGfM2R60Rjxk+YgOwDB8A4OcHw0SM4+/vD4vBhpZ9nQ6jZnRAiK6UOPmNtbU1zqROlEwgEDXagEzbz1xz1jcvlSjXKW3Z2tliP+JrHeOXlhYpz5/DyzTfBLSuDw4IFaLVwIfDypRLOqnGU1AkhspJp8JmIiAipt6UJXYiyPHnyBFFRUQgKCoKenl6dAW5qz9MOSJ/Uv//+e6SlpSEuLg7Lly8XrRf9kHBzQ1ZsLHjR0Wi5ZQssDxwAc+cOO/Nbx46qO2lCCJGDTEn92rVrUm1HtXWiTLGxsXjw4AESEhLA4/FEA9z06NFD4bInTpyIp0+f1rmvXax1QE8P+TNnoqx7dzjMnQv9O3eAbt2Ab74BJkxQOAZCCFEWmZI6TeJC1CE4OBiZmZnw9/eHo6MjioqKZB7gpiYulwtLS0vk5eXB3t5edJ29pszMTMTExGD+/Pmi2n5Zz554cuQInD//HNw//gCCg8FNTARn4UKFzo8QQpRFoWvqBQUFWLduHaZMmYIpU6Zgw4YNKCwsVFZsRMfIO70rn8/H5MmTYWdnV2eAG1XFERsbi9TUVMTFxQF4XXOvtLKC4LffgC+/BLhcZO/Zg5jevZF15Ijc8QjLp+vnhBBFyZ3Ur1y5And3d2zYsAH5+fnIz8/H+vXr4e7ujqtXryozRqIiisyhLs++TTlnu6KCg4PRtm1bBAUFIScnB9HR0Xj8+DE7mQyXC3z2GXD6NLbzeLhdWIhd//sfWmzYAFRWynW8hpK6tLfsEUKI3N+uc+bMwciRI/Ho0SMcPHgQBw8eRFpaGoYPH47w8HAlhkhURRVJvebtYdqMz+cjJCQEDg4OiI+Px7179/Drr7+Kb9S3Lyb//Tc8PTwwkmGwKSoKOW+8Ady7B0B5te/Y2FiakIUQIhWFauoLFiyAvv7ry/L6+vqYP38+rly5opTgiPZRd21c1h8VNe9Xr49wsJuAgIA6nUD5np4IOXoUCcOG4bq+Pr67cQPw9QW+/RaCBiaAkUVwcDBNyEIIkYrc37yWlpZIT0+vsz4jI0M0JjwhspD2NrTGypDlRwWXy4WpqWmD29vZ2WHWrFkwNTVFdHR0nWZwhmHw3oIFcJk4ESFvvsnexz59OrjDhkHvyRO5z0WIJmQhhEhL7m/PgIAAhISEID4+HhkZGcjIyMD+/fsxZcoUjB8/XpkxEh3G5XLB4/HA4XCUktTlwTAMgMavXSckJNQ75rydnR1mfv457E+dAtatA4yMwP39dzi++y4427apdWIYdbeeEEKajtz/y7/++mu8//77mDRpElxdXeHq6oqgoCCMGzcOa9asUWaMRIepKuHI07mssVnl/P39xcacl4jLBSIigOvXwbz5JrilpdCbNQsYOBB48EDW01AKSuqENB9y/y83NDTEpk2b8OLFCyQnJyM5ORn5+fnYsGEDjIyMlBkjIVKpmbwkJWhDQ0O4uLjA0NBQ4v6NzSpXe8z5Bnl4oPrPP/F8yRIwpqbA6dNAly7Apk1AdbU8p0cIIY2SOakLBAKsWbMGvXv3Rvfu3REZGYm2bduic+fOMDU1VUWMRAfVTMDCx9nZ2QrdulWzTEkJ2tDQEJ06dao3qSv92rWeHoqDg1F99SrQvz97rT08HOjdG0hOVs4xCCGkBpmT+pdffonPPvsM5ubmcHR0xKZNmxAaGqqK2IgOk5TUlXnrlrQJukmu47u7A4mJwNatgIUFcPEi4OcHzJkDFBer7riEkGZH5m+y3bt345tvvsGJEydw+PBh/Prrr9i7dy+NhEUU1lDzt6ruf5cnqct1/zmXC0yfDty9C/zvf2zHuY0bgQ4dgJ9+Av7rrEcIIYqQ+RsyPT0d7733nuj5oEGDwOFwaLQrorCGatfCXvLKSuqK1NAVGlSGzwfi44Hjx9ka/NOnwLhxyB4yBJsjI+n/ESFEITJ/o1VVVcHY2FhsnYGBASrlHB6T6A5t6mXdWFIXTuiSk5MjuuVOqYYMAW7eBL74AjAwQMypU3gUF4edH30ElJSIYqDhYQkhspD525dhGAQFBeH9998XLeXl5Zg+fbrYOtL81E7qykrywgQrTG61m+JV8WOi5oQuKvuhYmICREYCN28ipHdvdGEYTPv7b8DDA9izBzHffYe7d+8iNjZW+ccmhOgkmb+pAgMDYWtrCx6PJ1omTJgAPp8vto4QZSVbYYKNjY0Fl8uFtbU1rK2tFUrqje0jnNAlODgYAESTuuTk5NTZtqHXpOLhAftTpzDj22/h4OYGZGYCEydi2q+/ohOPh5CQEK1pASGEqJdM86kDoFoDaXLBwcHIzs5GcHCwKBkrSlhOVVWVxNeFE7oIr+/Hx8cjJSUF8fHx6Nq1q9i2NV/r0aOHfAFxOHg5eDCqP/wQ+tHRwMqVcExOxtfJyRDo64P71VeAEs6bEKLb6Kc/0Xi1E6yqCOdMr6qqEusIx+VyMXXqVNGkLjU19FpmZiY2btyIu3fvSt+pztgYWLgQuH8fgkmT2GPs3Qu0awcsWAC8eKH4iRJCdBYlddKsSeqMJhAIUFxcDDMzM9F99J6enpg9ezbs7OxE2xQUFAAAPD09ER4ejvbt24s1kQtHtduxY4fsPeUdHCCIiUHm4cNg+vQBysuBtWsBNzfg//4PKCtT+NwJIbpHK5L6smXLwOFwxBZPT091h0V0QM3hZGv2iG+sdzzDMCgsLIRAIIBAIMDLly9haWkptn1ISIioBi/vbXCvvL1R/eefwJEjgJcXUFAAzJ8PtG8PxMbSkLOEEDFakdQBoFOnTsjKyhItf//9t7pDIjqgsfHea6pvABxhs33thM3n8xEWFgY7OzvF7m3ncIBhw9ihZePiACcn4MkTYPJkoHNn9r53GvyJEAItSur6+vqwt7cXLTY2NuoOiegAWcZ7l7WXfX3JXm56ekBgIHDvHvD114CVFfDvv8AHH7CTxSQkUHInpJnTmqR+//598Pl8uLm54aOPPkJ6enqD21dUVKCoqEhsIaQpCZM6o+whYI2NgblzgYcPgWXLAB4PuH2bHX7W2xs4cICSOyE6oHYOq6ioaHQfrUjqPXv2RFxcHI4fP46tW7ciLS0Nffv2RXEDk2GsXr1a7L55JyenJoyYKFOTTLoig8zMTERFRcl/X7qytGgBLF0KPHrE/mtpCdy6Bfj7Az4+bHKna+6EaC0nJyexPLZ69epG99GMb8lGDB06FP7+/ujSpQuGDBmCo0ePoqCgAD/++GO9+yxatAiFhYWiJSMjowkjJsqkaUk9JiZGdF+6Ksh8/b1FC7bG/ugRsGQJm9xv3mSTe8eOQEwMIMUvfJXHSQiRSUZGhlgeW7RoUaP7aMa3pIxatGiB9u3bIzU1td5tjIyMYGlpKbYQogw1e7WrgjzJUiAQoIDDgUBYc1+yhL3mfu8eMGUKeyvc118DSrwMJYyz9n39hBDlqJ3DjIyMGt1HK5N6SUkJHjx4AAcHB3WHQpqhmr3aG5OXl4fo6OhGJ2VRdEhdsR8CVlbA8uXA48fAunWAoyM79OynnwIuLsDnnwNKvHRANXZCNIdWJPV58+bhzJkzePToEf755x+MGTMGenp6GD9+vLpDIxpKkSSpzGvmv/zyC1JSUhodXlmenvWNJlILCyAigu1Qt3Mn4OnJ3ue+ahWb3ENCgBs3ZDgbQoim04qk/uTJE4wfPx4eHh743//+h5YtW+LChQto1aqVukMjGkqRpB4bGyvTNfOGrvmPHDkSHh4eeO+992SaRrWxpC1T7djQEAgOZnvIHzoE9OzJXmPfuZPtLf/228Dhw9SpjhAdoBVJff/+/cjMzERFRQWePHmC/fv3w93dXd1hESVT93zswiFj33vvvUavmefk5CAqKgqZmZkNJnUbGxvMnDkTR48eFY1cJw2GYVBUVITq6upGk7fUCZ7LBUaPBs6fB86dY2+B09MDTp8GxowB2rZlm+vlGF+emuAJ0QxakdRJ86DupC4cMvbo0aONXjMXzswmbZIODg6WeuS6mmoOR1sfgUCAwsJC6e+H53CAN99kR6JLS2MnkLG2ZjvYzZsHtG4NfPwxcPWqUuMkhKgeJXVC/iMcMjY4OLjOkLC1f3AEBATAw8MDISEhEieFAdia/86dO5GTkyPTyHVNyskJWL2aHXZ2xw52fPmXL4Ht2wE/P6B7d3Z9SYm6IyWESIGSOiH/qZl4ayfx2s/t7OwQFhYGPp8vNilMTTt37sS9e/caHE+hPrIOMVvzcoBcTEzYW99u3AD+/JMdetbAALhyBZg2DXBwAKZPB65dk698QkiToKRONJ4qm+WVUXZ9k8IEBwfDzc0N/v7+MpfJMAxKS0thYWEBDofT6PayXg6oTXRNnGHYjnM//AA8fcpO89quHVtT//ZboGtXoFs3YMsW4PlzuY5FCFEdSupE42l6Uq+vad3R0RHTpk2Dvb29TOVlZmYiOjoa+fn5EmeFE25Ts8m/5uUAeUjs6NaqFXuNPSVFvPaelATMnImc1q2x4513cH3TJkRt2KD+YXMJIZTUCdE0wmFof/nllwa3qdnkX/NygNJxOOK19w0bAB8ffFdVhbupqVi4cCEytm/Hr598wiZ8ZU9go4Gotz/RVJTUCZERl8sFj8eTqllcHiEhIXB0dMTLly+RnJwssRNezSb/Jp1gplUrIDwcuHYNwcePo33XrvjaygrelZWYdf069N94g+1st2IFO0StjqKkTjQVJXVCZCRLk33tXvTS4PP5MDc3R1ZWFj777DOJnfBqNvmreoKZeuMcOBAf//QTOqWnY+TXX2NLu3Z4amgI3LnDjj3v4QH4+gJr1rC3zukwVSV5+vFAZEVJnRAVEtbqZb1m7+/vD3d3d6xatarR+9tVPcFMo/T1EZeejiQ3N2wPD2dHqhsyhB3YJjmZvQ/ezY0dyW79ekAHZ0ykpE40BSV1QmpR9yA4AHuNfPLkyfDx8Wn0/nZZJpipT3332ktL1FEvNJQdkvb4cSA7m+0xP2AAO5rdpUvA3LmAszOb4FevBv79VyuvwSv6fmnacYjuoKROSC2qSOrq+qEgbU2vvnvtpSWxo56NDXuPe2Ii28EuOhro25fteHfpEvDZZ+x8756ewIIF7PC1aqqRylojrv1+qSr5Kvp3Ic0PJXVC5CBNks7MzERMTIxofHhZknpOTg527tyJrKwsqWOSlJikTVb13WuvNPb2QGgocPYsm+C//RYYOpSdbObePWDtWnboWj6f/SFw5AhQWqqaWCSQNanXfr9UlXxV/nchOoeSOiFykCZJx8XF4dGjR4iLi5O5/ISEBDx48EDqfTkcDiwsLFBUVCTX9dcmHcbWwYFN3EePAs+esWPQjx8PWFqy87zv2AGMGAG0bAkMHszeQpeSolHN9LXfL1UlX40dXphoLErqhPxH2U3kU6ZMgZeXF6ZMmSLzvsKOckFBQRJfb2wYW61hacnOFrdvH5vgT54EZsxg53uvqAB+/52dE97TE3B3Bz75BPj11yatxUuDki/RFPrqDoAQTSFMjMqiyBe9sKOcg4MD8vLy6rxeM1ZFe0ZrTA9rQ0PgnXfYhWGAu3eBY8fY5exZ9ra4rVvZxdAQ6N0bGDiQXbp1A/Sb7uusqX5Eae2PNaI29EkhREHqHsZW0ePLPHWrhOMrfTAeDgfo0IGtpf/+OzvO/C+/sLV4V1fg1Svg1Clg8WKgVy926tgRI9im+hs3ZO5wJ2tHN0rqRFPRJ4UQBSn7i7dmklQ0qavylijhcfX19VWWeEStCKambNL+5hvg4UO2Fr9lCzB2LJvQi4vZznUREYC3N9sxLyCArdXfutVokle0o5suJF+NabEhCtHeTyAhOorL5cLKygqWlpYyJQlJCVzaZCXPyHdNkcgkJhoOhx2t7pNPgAMH2GvxV6+yM8oNHQqYmbHrfvyR3aZzZ/b2uhEj2NHt/vmHvV5fg6Id3XQhqT958gTLly/HkydP1B0KUYD2fgIJ0WFcLhfm5uYyJQlJCVzaZKVoUlJrUuNy2eFo581je9Tn5wN//w0sXw4MHAjG1BR48YKtyS9cyF6L5/GAt94CPv8cOHYMfGPjZt/Rje6J1w3UUY4QHRESEoKioiKxBC5NZz1hDX/OnDlwdnaW69jK7mQoJFdswk50vXsDAKrLypB78iTs7t+H3j//sAn/2TPgr7/YBewXoaOrK7hvvgm88QY74p23N2BkpPRzUoRAIEBRUZHMrTjSkPT5IdqHw8jbO0bLFBUVgcfjobCwEJaWluoOhzQDVVVVePLkCVq3bg19GXtmy7OvpH2kKWfp0qW4dOkSevTogeXLl8sUp6opI7Y67wHDAPfvswn977/ZJTW17o6GhoCPD9CjB5vke/QA2rVjm//VRJHPVGNU+YOByEaRfEU1dUI0UFM2Z2tyDU0lsXE4QPv27PJfuVU5OXh+/DhsHj6E3pUrwMWLbI/7S5fYJTqa3dfSkk30vr6v/+3YETAwUF58DVBGq0p9VNXaQpoW1dQJUZGmrvnIW1NXZe1PSN73QhmxpaenY8OGDQ0mwjrxMQx7X/zFi2xSv3iR7YxXq4MdALZG7+XFJnjh0qULYG4uV7wN0eRWFaI8VFMnRANRzec1YS92WTv/KUPNDmD1JcI6fysOh50u1s2NHcIWACor2Vnlrl1jp5QV/ltYyCb8q1fF92/Thk32NZf27RW6Tq/JrSqKevLkCbZs2YLQ0FC0bt1a3eFoLUrqhOgISU32unCrlaKUlggNDNgaeJcuQGAgu05Yo792TXzJymLvp3/4kB00R0hPj03sNRN9p07sELhStETo8nC0O3bsQHJyMnbs2EGtEAqg5ndCmrmmaH6XpglcEmVcwmiK86sjNxe4fZtdbt16vRQWSt5eX59N7B4er5f27dl/W7USdc7T5c5s8n5GdJEi+YqSOiHNXFMkCmmvBasiFrUkdUkYBsjMFE/ywsT/8mX9+7VoUTfRt23LXhrQoe8yXf7BIiu6pk4IkVtTXPuXtglcFdfeNeYSBIcDODqyy5Ahr9cLBOwc8ykpr5d799h/Hz8GCgrYjnoXL9Yts2XL19f+3d0hcHXFS3t7mHp5gevszDb3awnqg6IcVFMnhKictLVlaoKtpayMvYe+drJ/+JAdQKchBgbsFLbu7mynPWfn14uTE/vjooluxSOyoZo6IUQnSNNTvVkxMWHHru/cue5rxcVsJ72HD4EHD4CHD/Hk9m1sv34dH5eUwLGykv1BIGlgHYBtOeDz6yb7ms+trdU62A6RHSV1QojKSdsErsu3bCmdhcXr3vj/2bF0KS6bmIDTrRuWT5v2OuGnpQEZGUB6+ut/X71im/2fPgXOn5d8DCMjNvELFwcH8efCxdKSkr+GoOZ3QojG0JhObVpK6ssXAgHbfJ+eLp7oay45OdIf2NRUPOk7OAC2tuxiZ/f6sa0tuy1pEPV+lwIldUI0HyV1xSi1B3lFBXu/fWYmu9R8XHMpKJCtXDMzycm+5mJjw3YCtLZulj8C6Jo6IUQnaExPdS2l1B7kRkaAqyu7NOTly7rJPyuLvVe/5pKTw/5QKC1lLwekpUkXh7Hx6wTfsqV0j62tm20nQErqhBCNQbc1aSFTU7aHvbt7w9sxDFBSIp7kayd+4frnz9mlqgooL3997V8WZmbsPf4tWgA8nvi/0qwzMdHKfgKU1AkhhKgeh8N27rOwaPwHAPD6R4AwwT9/DuTniz+XtE54OaC0lF1k/TEgZGDAJnhLy9dxW1iwE/U09Ly+bZpozABK6oQQQjRPzR8BjV0CqKm6Gnjxgk3uhYXsvzUfN7SusJBdBAJ2Ap9nzxofD0BaJiavz8fM7PViair+3MxMoR8AlNQJIYToDj09tqOdjY18+wsEbAtBzURfUsKOCyBcZHleVcWWW1bGLrm5SjtVSbQqqW/ZsgX/93//h+zsbHh7e2Pz5s3o0aOHusMihBCiK7hctsnd0pIdjEcRDMN2Dqyd6IWXBl6+fP245vLiBbBnj1yH1JqkHh8fj4iICGzbtg09e/bExo0bMWTIEKSkpMDW1lbd4RFCCCHiOBy2976xMTvbnrSKiuRO6lpz38j69esxdepUBAcHo2PHjti2bRtMTU2xc+dOdYdGCCGEaAStqKm/evUKSUlJWLRokWgdl8vFoEGDcL6e4Q0rKipQUVEhel743zzGRUVFqg2WEEIIUYAwTwnzlpCRkRGMjIwa3FcrknpeXh6qq6thZ2cntt7Ozg53796VuM/q1aslTgjhpOg1EkIIIaQJ1B7qd+nSpVi2bFmD+2hFUpfHokWLEBERIXpeUFAAFxcXpKeng8fjqTEyxXTv3h2XL1/W+mMqWqY8+8uyj7TbNrZdfa8XFRXByckJGRkZWj9sMX0m5d9f2Z9JRbahz6TmHJNhGHTt2hVJSUlioys2VksHtCSp29jYQE9PDzm1JhjIycmBvb29xH3qa6bg8Xha/YHV09Nr8vhVcUxFy5Rnf1n2kXbbxrZr7HVLS0ut/jwC9JlUZH9lfyaVsQ19JjXjmEZGRnKNrqgVHeUMDQ3h5+eHxMRE0TqBQIDExET06tVLjZE1vdDQUJ04pqJlyrO/LPtIu21j26nj79XU6DMp//7K/kwqaxttpwufSXnL05pZ2uLj4xEYGIhvv/0WPXr0wMaNG/Hjjz/i7t27da61S0KztBFNQp9HomnoM6kbtKL5HQACAgLw7NkzLFmyBNnZ2fDx8cHx48elSugA25SxdOlSqa5JEKJq9HkkmoY+k7pBa2rqhBBCCGmYVlxTJ4QQQkjjKKkTQgghOoKSOiGEEKIjKKkTQgghOoKSOoAjR47Aw8MD7dq1w3fffafucAjBmDFjYGVlhXHjxqk7FEKQkZGB/v37o2PHjujSpQsSEhLUHRKpR7Pv/V5VVYWOHTvi1KlT4PF48PPzwz///IOWLVuqOzTSjJ0+fRrFxcXYtWsXDhw4oO5wSDOXlZWFnJwc+Pj4IDs7G35+frh37x7MzMzUHRqppdnX1C9duoROnTrB0dER5ubmGDp0KE6ePKnusEgz179/f1hYWKg7DEIAAA4ODvDx8QEA2Nvbw8bGBvn5+eoNikik9Un97NmzGDFiBPh8PjgcDg4fPlxnmy1btsDV1RXGxsbo2bMnLl26JHotMzMTjo6OoueOjo54+vRpU4ROdJSin0lClE2Zn8mkpCRUV1fTjJcaSuuTemlpKby9vbFlyxaJr8fHxyMiIgJLly7F1atX4e3tjSFDhiA3N7eJIyXNBX0miaZR1mcyPz8fkyZNwvbt25sibCIPRocAYA4dOiS2rkePHkxoaKjoeXV1NcPn85nVq1czDMMw586dY0aPHi16ffbs2czevXubJF6i++T5TAqdOnWKGTt2bFOESZoReT+T5eXlTN++fZndu3c3VahEDlpfU2/Iq1evkJSUhEGDBonWcblcDBo0COfPnwcA9OjRA7du3cLTp09RUlKCY8eOYciQIeoKmeg4aT6ThDQlaT6TDMMgKCgIAwYMwMSJE9UVKpGCTif1vLw8VFdX15n0xc7ODtnZ2QAAfX19rFu3Dm+//TZ8fHwwd+5c6vlOVEaazyQADBo0CP7+/jh69Chat25NCZ+ojDSfyXPnziE+Ph6HDx+Gj48PfHx8cPPmTXWESxqhNbO0qdLIkSMxcuRIdYdBiMgff/yh7hAIEenTpw8EAoG6wyBS0Omauo2NDfT09JCTkyO2PicnB/b29mqKijRn9JkkmoY+k7pFp5O6oaEh/Pz8kJiYKFonEAiQmJiIXr16qTEy0lzRZ5JoGvpM6hatb34vKSlBamqq6HlaWhqSk5NhbW0NZ2dnREREIDAwEN26dUOPHj2wceNGlJaWIjg4WI1RE11Gn0miaegz2Yyou/u9ok6dOsUAqLMEBgaKttm8eTPj7OzMGBoaMj169GAuXLigvoCJzqPPJNE09JlsPpr92O+EEEKIrtDpa+qEEEJIc0JJnRBCCNERlNQJIYQQHUFJnRBCCNERlNQJIYQQHUFJnRBCCNERlNQJIYQQHUFJnRBCCNERlNQJIYQQHUFJnRBCCNERlNQJIYQQHUFJnRBCSL0yMjLQv39/dOzYEV26dEFCQoK6QyINoAldCCGE1CsrKws5OTnw8fFBdnY2/Pz8cO/ePZiZmak7NCIB1dSJRunfvz/Cw8PVHUaT0LVzVef5PH/+HLa2tnj06JFajq/LHBwc4OPjAwCwt7eHjY0N8vPzRa9/8MEHWLdunZqiI7VRUidyefbsGWbMmAFnZ2cYGRnB3t4eQ4YMwblz5wBoR8LicDgNLsuWLZOqHFWea1BQEEaPHi3TPup67w8ePIgVK1YoPYZ+/fph8uTJDW7z5ZdfYtSoUXB1dVXoWPK8381JUlISqqur4eTkJFq3ePFifPnllygsLFRjZERIX90BEO00duxYvHr1Crt27YKbmxtycnKQmJiI58+fqzs0qWVlZYkex8fHY8mSJUhJSRGtMzc3V0dYWsva2lrpZTIMg2vXrsHf37/ebV6+fImYmBicOHFC6cdvCq9evYKhoaFaY/Dx8UFVVVWd9SdPngSfzwcA5OfnY9KkSdixY4fYNl5eXnB3d8eePXsQGhraJPGSBjCEyOjFixcMAOb06dMSXw8MDGQAiC1paWkMwzBMeXk5M2vWLKZVq1aMkZER07t3b+bSpUuiffv168fMnj1b9PzIkSOMpaUls2fPHoZhGKa6uppZtWoV4+rqyhgbGzNdunRhEhISxI6fkJDAeHl5McbGxoy1tTUzcOBApqSkpMFzio2NZXg8nsTXGoq5oXM9duwY07t3b4bH4zHW1tbMsGHDmNTU1HrPtb73ctSoUWL7zJo1i/n0008ZKysrxs7Ojlm6dKnY9pLiaex9a6xcad5X4fnUF8OuXbsYa2trpry8XKzcUaNGMRMmTJB4/ikpKQwA5p9//qn3PUpISGBatWoltq5fv37MzJkzmdmzZzMtWrRgbG1tme3btzMlJSVMUFAQY25uzri7uzNHjx6t9/2WJ16Gke7vHhoaysyePZtp2bIl079/f4Zh2M/2mjVrGHd3d8bQ0JBxcnJiVq5cKTrHht77xv6+DZUtjfLycqZv377M7t27Jb6+fPlypk+fPlKXR1SHkjqRWWVlJWNubs6Eh4fX+cJjGIYpKChgevXqxUydOpXJyspisrKymKqqKoZhGCYsLIzh8/nM0aNHmdu3bzOBgYGMlZUV8/z5c4ZhxBPd3r17GQsLC+bXX38Vlb1y5UrG09OTOX78OPPgwQMmNjaWMTIyEv3AyMzMZPT19Zn169czaWlpzI0bN5gtW7YwxcXFDZ5TQ0m9oZgbOtcDBw4wP/30E3P//n3m2rVrzIgRI5jOnTsz1dXVdc61PpKSuqWlJbNs2TLm3r17zK5duxgOh8OcPHmywfe+sfetsXKleV+F51NfDC9fvmR4PB7z448/ivbJyclh9PX1mT///FPi+e/bt4/R09NjSktL632PwsLCmHfffVdsXb9+/RgLCwtmxYoVzL1795gVK1Ywenp6zNChQ5nt27cz9+7dY2bMmMG0bNlSrOya77c88TKMdH93c3Nz5tNPP2Xu3r3L3L17l2EYhpk/fz5jZWXFxMXFMampqcxff/3F7NixQ6r3vrG/b31lS0MgEDAffPBBnR95NR07dowxNDSU+H1AmhYldSKXAwcOMFZWVoyxsTHz5ptvMosWLWKuX78uel1SwiopKWEMDAyYvXv3ita9evWK4fP5zNq1a8X2i46OZng8nlhrQHl5OWNqalqn1hYSEsKMHz+eYRiGSUpKYgAwjx49kul86kvqssTcmGfPnjEAmJs3b0q9n6SkXrtG1L17d2bBggVi29QsV5r3rbFypXlfax63vnObMWMGM3ToUNHzdevWMW5uboxAIJBY5rx585iOHTvWe0yGYWvOkydPrhNLzfOpqqpizMzMmIkTJ4rWZWVlMQCY8+fPi9bVfr9ljVcSSX93X19fsW2KiooYIyMjiYm2sfe+sb9vQ2VL46+//mI4HA7j7e0tWm7cuCG2zfXr1+X6f0eUj66pE7mMHTsWw4YNw19//YULFy7g2LFjWLt2Lb777jsEBQVJ3OfBgweorKxE7969ResMDAzQo0cP/Pvvv6J1Bw4cQG5uLs6dO4fu3buL1qempuLly5d45513xMp99eoVfH19AQDe3t4YOHAgOnfujCFDhmDw4MEYN24crKys5DpPaWOW5P79+1iyZAkuXryIvLw8CAQCAEB6ejq8vLzkigcAunTpIvbcwcEBubm59W4vzfvWWLnKel+nTp2K7t274+nTp3B0dERcXByCgoLA4XAkbn/16lX4+fk1WGZZWRmMjY3rrK95Pnp6emjZsiU6d+4sWmdnZwcADb53ssYLSPd3r31O//77LyoqKjBw4MA65TX23jf2922obGn06dNHdA71MTExAcD2byDqRUmdyM3Y2BjvvPMO3nnnHXzxxReYMmUKli5dWm9Sl5avry+uXr2KnTt3olu3bqIv0JKSEgDAb7/9BkdHR7F9jIyMALBf3r///jv++ecfnDx5Eps3b8bnn3+Oixcvok2bNgrFJasRI0bAxcUFO3bsAJ/Ph0AggJeXF169eqVQuQYGBmLPORxOg1+60rxvjZWrrPfV19cX3t7e2L17NwYPHozbt2/jt99+q3f7q1evYsSIEQ2WaWNjgxcvXtRZL+l8aq4Tfq4aeu9kjReQ7u9e+x5vYVKUpLH3vrG/b0FBQYPxKoPwFrdWrVqp/FikYXRLG1Gajh07orS0FABgaGiI6upqsdfd3d1haGgouu0NACorK3H58mV07NhRbLtTp07h559/xqxZs8TKNzIyQnp6Otq2bSu21LzFhsPhoHfv3li+fDmuXbsGQ0NDHDp0SK5zkiZmSef6/PlzpKSkYPHixRg4cCA6dOggMfGoQu14pH3fGiPL+yrpPRGaMmUK4uLiEBsbi0GDBtUbw8OHD1FQUNBoTd3X1xd37tyR+jxkJW28gPx/93bt2sHExASJiYkSX2/ovW/s79tY2cpw69YttG7dGjY2Nio7BpEO1dSJzJ4/fw5/f39MnjwZXbp0gYWFBa5cuYK1a9di1KhRAABXV1dcvHgRjx49grm5OaytrWFmZoYZM2bg008/hbW1NZydnbF27Vq8fPkSISEhYsdo3749Tp06hf79+0NfXx8bN26EhYUF5s2bhzlz5kAgEKBPnz4oLCzEuXPnYGlpicDAQFy8eBGJiYkYPHgwbG1tcfHiRTx79gwdOnSQ61yliVnSuVpZWaFly5bYvn07HBwckJ6ejoULFyr2xktJUjyNvW+NkfV9lRQDl8vWIT788EPMmzcPO3bswO7du+s9ZlJSEjgcjmjgk/oMGTIEixYtwosXL+S+zNIQaeMFIPff3djYGAsWLMD8+fNhaGiI3r1749mzZ7h9+za8vLwafO+l+X9RX9m1/9/J66+//sLgwYOVUhZRkLov6hPtU15ezixcuJDp2rUrw+PxGFNTU8bDw4NZvHgx8/LlS4Zh2FuR3njjDcbExETsNq+ysjJm1qxZjI2NjVS3tN25c4extbVlIiIiGIZhe+Ju3LiR8fDwYAwMDJhWrVoxQ4YMYc6cOSPafsiQIaLbz9q3b89s3ry50XNqqPd7YzHXd66///4706FDB8bIyIjp0qULc/r0aQYAc+jQIYnnKomkjnK19xk1ahQTGBjYYDyNvW+NlSvN+1qzjPreE6GJEydKvF2spoULFzLt27dv8P0R6tGjB7Nt2zaJsQi5uLgwGzZsEFtX8+8hjGvs2LF1ypcmXiF5/+7V1dXMypUrGRcXF8bAwIBxdnZmVq1aJdV739jft76ylaGsrIzh8XhiHQ6J+tDY74SQJjdw4EB06tQJUVFRSinvt99+w6effopbt26JWgTk8e6776Jt27aIjo4WW6/seHXJ1q1bcejQIZw8eVLdoRBQ8zshpAm9ePECp0+fxunTp/HNN98ordxhw4bh/v37ePr0qUz9BGrGde7cOZw+fRrTp09Xeby6xMDAAJs3b1Z3GOQ/VFMnhDQZV1dXvHjxAl988QXmzZun7nBExowZg8uXLyMwMBArV64U9YzX1HgJqQ8ldUIIIURH0C1thBBCiI6gpE4IIYToCErqhBBCiI6gpE4IIYToCErqhBBCiI6gpE4IIYToCErqhBBCiI6gpE4IIYToCErqhBBCiI6gpE4IIYToCErqhBBCiI6gpE4IIYToiP8HxwU2l4CF0+gAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "PvI_plot_Step2 = HtoH.PvI(CatP, CatdP, CatI, Iscale=[1.0,400.0], Pscale=[0.0,30.0], showfit='true', weighted='false', errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "59006181",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Power law index: -0.786452214396975 ± 0.1084423950771106\n",
      "Coefficient: 1.5327408328944348e+18 ± 8.412557689506365e+18\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfUAAAF2CAYAAAB3bAFnAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABriElEQVR4nO3deVxU9frA8c8MqwugiLII7qm57+bVct9aXDOyTEHstpgbbVJduWpJWmruVlcxy5Ios0VvpubSr9TccEnFVAyUxQwBUUFhzu+Pc2cCWWeYYRae9+t1XjJnzvLMODPP+X7Pd9EoiqIghBBCCLuntXYAQgghhDAPSepCCCGEg5CkLoQQQjgISepCCCGEg5CkLoQQQjgISepCCCGEg5CkLoQQQjgISepCCCGEg5CkLoQQQjiIKpPUc3Nz+fe//01ubq61QxFCWJh830VVZRNJfdWqVbRr1w5PT088PT3p0aMH//3vfw3P5+TkMHnyZOrUqUPNmjUZPXo0aWlpRp0jNzeX2bNny5dciCpAvu+iqrKJpB4YGMjbb7/N4cOHOXToEP369WP48OH89ttvAMyYMYNvv/2W2NhY9uzZQ3JyMqNGjbJy1I5vxYoV1g7BKNaO19LnN/fxzXG8ihzD2H2t/f/r6Ozt/bV2vDb7fVdsVO3atZX//Oc/SkZGhuLi4qLExsYanjt9+rQCKPv27Sv38TIzMxVAyczMtES4Dunee++1dghGsXa8lj6/uY9vjuNV5BjG7mvM9vJ9N561vz/Gsna8tvp9dzbjhYVZ5OfnExsby40bN+jRoweHDx/mzp07DBgwwLBNy5YtadCgAfv27eO+++4r9ji5ubmFqt4yMjIAyMzMtGj8jiQ/P5+srCxrh1Fu1o7X0uc39/HNcbyKHMPYfY3ZXr+dPX1+rc3a3x9jWTteS55fURRyc3PJyMhAq/27Qt3NzQ03N7cyd7YJx48fV2rUqKE4OTkpXl5eypYtWxRFUZQNGzYorq6uRbbv2rWr8sorr5R4vMjISAWQRRZZZJFFFodYIiMjy8ylNlNSb9GiBXFxcWRmZvLFF18wYcIE9uzZY/LxIiIiCA8PNzzOzMykQYMGJCUl4enpaY6Q7cPTT8Pnn8Mjj8Ann1g7GiEqRV5eHpcvX6Z+/fo4O9vMz5wQ5ZKVlUVQUBCJiYl4eXkZ1pdZSgds5tPu6upKs2bNAOjcuTMHDx5kyZIlBAcHc/v2bTIyMqhVq5Zh+7S0NPz8/Eo8XknVFPoW9lXGrFlqUv/uO0hOhpYtrR2REBaXl5dHVlYWnp6ektSF3fLy8jI6X9lE6/fi6HQ6cnNz6dy5My4uLuzcudPwXHx8PImJifTo0cOKEdqJ1q1h+HBQFJg/39rRCCGEsCCbuISNiIhg6NChNGjQgOvXr/Ppp5+ye/dutm3bhpeXF2FhYYSHh+Pt7Y2npydTpkyhR48eJTaSE3eJiICvv1ar32fPhgYNrB2REEIIC7CJpH7lyhXGjx9PSkoKXl5etGvXjm3btjFw4EAAFi9ejFarZfTo0eTm5jJ48GBWrlxp5ajtSPfu0K8f/PgjvP02yHsnhBAOSaMoimLtICpDVlYWXl5eZGZmVq176np79kCfPuDqCufPQ2CgtSMSwmLy8vK4dOkSgYGBck9d2J2K5CubvacuzKx3b3W5fVvurQshhIOSpF6VzJql/vvhh2pLeCGEEA5FknpV0rcv9OoFubmwYIG1oxFCCGFmktSrEo0GIiPVv99/H1JSrBuPEEIIs5KkXtX07w89ekBODrzzjrWjEUIIYUaS1KuagqX11avByHnphRBC2C5J6lXRoEFq3/Vbt+Ddd60djRBCCDORpF4VaTR/t4RfuRKuXLFuPMJm6HQ6MjIy0Ol01g5FCGECSepV1dCh0KUL3LwJixZZOxphIySpC2HfJKlXVQXvrS9fDlevWjceIYQQFSZJvSp76CHo1Alu3JB+68KuSQ2DECpJ6lWZRgNz56p/L18u/daF3ZKkLoRKknpVN3So2m/91i2IirJ2NEIIISpAknpVp9HAm2+qf7//PiQmWjceIcwgOTmZJUuWkCxzHIgqRpK6UOda79tXncFNn+CFsGNr1qzhzJkzrFmzxtqhCFGpJKkLlf7e+tq1cO6cdWMRooLCwsJo2bIlYWFh1g5FiEolSV2oevZU76/n58Ps2daORogKCQgIYNq0aQQEBFg7FCEqlSR18Td9aX3DBjh1yrqxCCFswrp169BoNGg0GsO6Pn36oNFoCAkJsV5goliS1MXfOneGkSNBUf4emEYI4bD0ybm4ZfPmzQDUrVuX7t270717d+sGa0YhISFoNBr69OlTru1Leo/eeOMNywZqAmdrByBszOzZsHkzfPEFHD0KHTtaOyIhSqTT6cjKyqJ69erWDsWuubq60vGu77q3tzcADz30EA899JA1wrI5HTp0wM3NzfA4KCjIitEUT0rqorC2beHxx9W///Uv68YiRBlk0Bnz8Pf3Z//+/YWWBx54ACi++r04EyZM4J577sHDwwNXV1caNmzI1KlTycrKMmyjLyE3atSIjz76iKCgIDw9PZk+fTq3bt1i+vTpeHl50aBBA1avXl3o+MnJyUycOJGAgABcXV1p0qQJc+fOJS8vz7CNvuZh/PjxREZG4u/vT+3atRk3bhzXr18HMJwbYM+ePYbXtnv37jLfp6+++qrQe/TMM8+U6/2tTJLURVH//jc4OcGWLfDTT9aORghhB77++muuXbtG06ZNCQoKIjExkWXLlhXbAyElJYXnn38eNzc3rl+/zpIlS+jcuTPR0dF4eXmRlJTE5MmTOX36NAB//fUX9913H9HR0WRnZ3PvvfeSlJTErFmz+Oc//1nk+Bs3bmTx4sVUq1aNjIwMNmzYwNtvvw1Ax44d8fHxAcDDw8Nwa8HT07PM19ilSxeqV69O69atefvtt8nNza3IW2YRktRFUc2bg/6L+Oqr6j12IYTD+uOPP4rcLzbWnj17uHr1KnFxcZw/f57XX38dgM2bN5OTk1No29u3b/PDDz9w9uxZQxX2uXPnOHbsGKdOncLd3R2dTseePXsAWL58OUlJSfj6+nL+/HmOHTvGF198Aag1Cefu6obr7u7O6dOnOXfuHJ07dwZg586dgFra1t9O6NSpk6HU3alTp1JfX+3atQkMDMTNzY1Tp04RERHB+PHjjX6fLE2SuiheZCRUqwb79sE331g7GiGEBbm6uhpKrKY2ituxYwdt2rShWrVqaDQa3nrrLQDy8vL4888/C21bu3ZtevbsiVarpUGDBgC0adOGRo0aUbNmTerVqwdAWloaAL/++qvhcb169dBoNIwYMQIARVE4cOBAoeP369eP+vXro9VqadmyZaFjmWL//v389ddfxMXFcfnyZfr16wfA559/TlJSksnHtQRpKCeKFxAA06er48FHRKgzujnLx0UIR6S/p26qDRs28NJLLxmOFRQUxNWrV7lw4QIA+fn5hbYvWNXt/L/flYLr9DUFyl21hB4eHrRq1arI+e9uKFmrVq0ix7/7WMYoeJFTvXp1Ro4cyY8//ghAUlKSTTWYk5K6KNkrr4C3N5w+DevXWzsaIYSN0l8QeHh4kJCQwIEDBxg0aJDZjt+1a1dATdAbN240VJlv376d559/npEjRxp1PP1FwI0bN8rcdu/evXzxxReGC5OcnBy+/vprw/MNGzY06tyWJkldlKxWLXjtNfXvyEh1JjchhLhLu3btALh+/TpNmjShSZMmfP7552Y7/uTJk6lfvz7Xrl2jRYsWdOjQgaZNm1KnTh0mTJhg9PH0VfKHDh2ibdu23Hfffdwq4fftwoULjBkzBi8vL9q1a0dAQAA7duwAIDQ0lPr165v+wixAkroo3eTJEBQEly6pc64LmyHduYStCAsLIzw8HB8fH65fv06fPn2YM2eO2Y5ft25d9u/fT2hoKHXq1OG3337j1q1b3H///SxevNjo402cOJHRo0fj5eXFyZMnOXDgQJFbBHq9evXi2WefpUGDBiQkJKDT6ejcuTOrV6/mgw8+qOhLMzuNUpEbDXYkKysLLy8vMjMzy9V1QRSwbh2EhkLt2nD+vPqvsLq8vDwuXbpEYGCg4b6hLR7TkvTx+vn5kZqaaoi7rNehH7TG09MTrVbKNsK2VCRfyadZlO2pp6B1a7h2DebPt3Y04n90Oh3Z2dlSUjeB1HIIRyVJXZTNyUltBQ+wZIlaFS+sTpK6EOJuktRF+Tz8MPTqBTk5YIOTGAghhJCkLspLo4GFC9W/169XJ3sRwkhS7S2EZUlSF+XXrRuMHasOG/vSSzJ8rDCaJHUhLEuSujDOvHng5gY//ghbt1o7GiGEEAVIUhfGadQIpk1T/375ZSgw7aEQQgjrkqQujBcRAXXqqMPH/uc/1o5GmEiqwoVwPJLUhfFq1VLnXAeYNQuysqwZjTBRVUrqVem1iqrNJpJ6VFQUXbt2xcPDg3r16jFixAji4+MLbdOnT58i8/0+++yzVopY8Mwz6rzrf/4pA9IIk1RmopWkbpwVK1bQqFEj3N3d6d69u2Hq0+I0atSoyG+zRqNh8uTJhm2uX7/O9OnTadiwIdWqVeMf//gHBw8etHrsoE7Y8sgjjxAQEIBGo2Hz5s0mbWMrbCKp79mzh8mTJxtm3blz5w6DBg0qMoPO008/TUpKimFZsGCBlSIWuLiA/v1ftAgSE60bj7AZ5U2gkmhtU0xMDOHh4URGRnLkyBHat2/P4MGDuXLlSrHbHzx4sNDv8vbt2wEYM2aMYZtJkyaxfft2Pv74Y06cOMGgQYMYMGAAly9fLjOePn36sG7dOovEDupMbe3bt2fFihUV2sZmKDboypUrCqDs2bPHsK53797KtGnTTD5mZmamAiiZmZlmiFAoiqIoOp2i9O6tKKAojz9u7WiqnFu3bin79u1Tbt26ZdL+d+7cURISEpQ7d+6Uus4cxy3uOXOe69atW4WOdfexy3pcGe655x7lvvvuU27evGlYp9PplO7duyszZ86stDjK0q1bN2Xy5MmGx/n5+UpAQIASFRVVrv2nTZumNG3aVNHpdIqiKMrNmzcVJycn5bvvviu0XadOnZTXX3+9zOP17t1biY6OrpTYAeWrr76q8DYVVZF8ZRMl9btlZmYC4O3tXWj9hg0b8PHxoU2bNkRERHDz5s0Sj5Gbm0tWVlahRZiZRgPvvaf+u3Ej/PSTtSMSwmbFxMRw5MgRfv75Z8O6DRs28Mcff/CafopjYN68edSsWbPUJdFCNWO3b9/m8OHDDBgwwLBOq9UyYMAA9u3bV679P/nkEyZOnIhGowHUSXfy8/Nxd3cvtG21atX4v//7P5uJ3RbdncNyc3PL3MfmkrpOp2P69On07NmTNm3aGNY/8cQTfPLJJ+zatYuIiAg+/vhjxo0bV+JxoqKi8PLyMixBQUGVEX7V06EDPP20+vf06SBVqUIUq2PHjnTo0IEzZ84AcPPmTSIiInjzzTfx8PAwbPfss88SFxdX6hIQEGCRGK9evUp+fj6+vr6F1vv6+pKamlrm/ps3byYjI4OQkBDDOg8PD3r06MHcuXNJTk4mPz+fTz75hH379pGSklLkGHdf1Pz00088++yzZV7UVDR2WxQUFFQoj0Xp5+Aohc3NrTh58mROnjxZ5Arun//8p+Hvtm3b4u/vT//+/Tl//jxNmzYtcpyIiAjCw8MNj7OysiSxW8rcuWpJ/cgRdZrWiROtHZEQNql58+aGRsALFizAx8eH0NDQQtt4e3sXqaU01syZM5lfRgPW06dP07Jlywqd525r1qxh6NChRS46Pv74YyZOnEj9+vVxcnKiU6dOjB07lsOHDxc5xrPPPstjjz1mePzkk08yevRoRo0aZVhnqYsaW5OUlFRo6lU3N7cy97GppP7CCy/w3XffsXfvXgIDA0vdtnv37gCcO3eu2KTu5uZWrjdAmEG9ehAZCS++qPZhf/RRkDnrqyRp/Fa6Fi1asHfvXi5dusQ777zDli1bisznPm/ePObNm1fqcU6dOkWDBg1KfP7FF18sVFouTpMmTYqs8/HxwcnJibS0tELr09LS8PPzK/V4f/zxBzt27GDTpk1FnmvatCl79uzhxo0bZGVl4e/vT3BwcLEx3H1RU61aNerVq0ezZs1KPX9FYrdVnp6eRs+nbhNJXVEUpkyZwldffcXu3btp3LhxmfvExcUB4O/vb+HoRLm88AK8/z6cPQtvvvl3y3hRpeh0OkObGFFU8+bN+fDDD5k5cyaDBg2iT58+Rba5u6RanLJKqnXr1qVu3bpGx+fq6krnzp3ZuXMnI0aMANT/0507d/LCCy+Uum90dDT16tXjoYceKnGbGjVqUKNGDa5du8a2bdvM2oOpIrE7EptI6pMnT+bTTz/l66+/xsPDw3D/w8vLi2rVqnH+/Hk+/fRTHnzwQerUqcPx48eZMWMGDzzwAO3atbNy9AIAV1dYvBgeekhtPPf003DPPdaOSgib0rx5c5KSkvjiiy84efJksduYo/q9IsLDw5kwYQJdunShW7duvPfee9y4caPQbYLly5fz1VdfsXPnTkBNntHR0UyYMAFn56JpZdu2bSiKQosWLTh37hwvv/wyLVu2LHLrASA7O5vs7GzD440bNwIUui9et25dnJycjI797rj15zt37pzhcUJCAnFxcXh7extqQ8qzjc0we1t8EwDFLvpuDImJicoDDzygeHt7K25ubkqzZs2Ul19+2ajm/tKlrZIMGaJ2cXvkEWtH4vBssUvbnTt3lHPnzinnzp2r1C5t58+fV6ZPn6788ccfxb4OW+jSpiiKcuPGDUWj0SgvvvhipZ7XWMuWLVMaNGiguLq6Kt26dVP2799f6PnIyEilYcOGhsfbtm1TACU+Pr7Y48XExChNmjRRXF1dFT8/P2Xy5MlKRkZGsdtGRkaWmBP0S0JCgkmx3x23oijKrl27ij3HhAkTjNrGnCqSrzSKUjXmz8zKysLLy4vMzEyj71EII5w+De3aqRO9/Pe/MGSItSNyWDk5OcTFxdGhQ4ci3YXKIy8vj0uXLhEYGGgoXRW3zthj/vHHHwA0bNiwyDEKHh+o0LkKHm/NmjUcOnSIbt26MXv27CKvo6zHlSU9PZ06depw7NgxqWUUJapIvrK5Lm3Czt17L0ydqv49ZQqUo1+lEBUVGhpKy5YtCQsLs3YopTp27Biurq7ce++91g5FOKgKJfU7d+6QlJREfHw86enp5opJ2LvISPDzg3PnYOFCa0cjqoCAgACmTZtm812djh07RqtWrXBxcbF2KMJBGZ3Ur1+/zqpVq+jduzeenp40atSIe++9l7p169KwYUOefvppiw3UL+yEpye8+67695tvyrjwQvzP9OnTOXr0qLXDEA7MqKS+aNEiGjVqRHR0NAMGDGDz5s3ExcVx9uxZ9u3bR2RkJHl5eQwaNIghQ4bw+++/WypuYeueeAIeeABu3YICgwAJIYSwHKNaiBw8eJC9e/fSunXrYp/v1q0bEydOZPXq1URHR/PTTz9xj3Rrqpo0Gli+HDp2hC+/hO3bYeBAa0clhBAOzaik/tlnn5VrOzc3N5nrXEDbtuqgNEuWqI3mjh9X+7MLIYSwCLO0fr969Spbtmzhm2++KXaAflGFzZ4Nvr4QH68OTiOEEMJiKpzUv/zyS5o1a8bs2bOJjIykadOmREdHmyM24Qi8vP4eMnbOHPhfH2YhhBDmZ3RSLzh8H8Ds2bP59ddf+fXXXzl69CixsbG8/vrrZgtQOICnnlIbzd28+XcfdiFMIBPGCFE6o5N6586d+frrrw2PnZ2duXLliuFxWloarnLf1CGY7QdUo4FVq8DZGb75Bgp8foQwhiR1IUpndFLftm0bH3zwASNHjiQ5OZklS5YQHByMn58fPj4+zJw5k5UrV1oiVlHJzPoD2qoVvPyy+veUKXBXjY8QQoiKMzqpN2rUiC1btvDYY4/Ru3dv4uLiOHfuHNu3b2fHjh0kJiby4IMPWiJWYe/eeAMaNYKkJLUBnRBCCLMyuaHc2LFjOXjwIMeOHaNPnz7odDqTJ5YQVUT16rBihfr34sVqFzchhBBmY1JS37p1KwsXLuTQoUP85z//YcGCBTz55JO8/PLL3Lp1y9wxCkfy4IMwejTk58Ozz4LcGxUWoNVqqVWrFlqtzFklqhajP/EvvvgioaGhHDx4kGeeeYa5c+fSu3dvjhw5gru7Ox07duS///2vJWIVjuK996BmTdi3D95/39rRCAckSV1UVUZ/4tetW8fWrVvZuHEjBw8e5OOPPwbA1dWVuXPnsmnTJubNm2f2QIUDCQwE/WfklVfUe+xCCCEqzOikXqNGDRISEgBISkoqcg+9VatW/PTTT+aJTjiu55+HHj3UVvDPPQeKYu2IhIVJdzQhLM/opB4VFcX48eMJCAigd+/ezJ071xJxCUfn5AT/+Y86FvyWLVDOeQWE/dLpdGRmZqKY8QJO3602OTm52PPJRYSoaoxO6k8++SRJSUl8/fXXXLx4keHDh1siLlEVtGqldnMDmDYNrl61bjzC7qxZs4YzZ84UOzR1aUm9tIsBIeyZUUk9MTERgDp16tC1a1dq1apV6vaXL182OTBRRbz6qjqb29WrMH26taOp8uwt2YWFhdGyZUtCQ0ON2k9/MbBmzRoLRSaEdRiV1Lt27cozzzzDwYMHS9wmMzOTDz/8kDZt2vDll19WOEDh4FxdYc0a0Gphwwa1Kl5Yjb0lu4CAAKZNm0ZAQEC5tteX3kNDQ2nZsiVhYWEWjlCIymXUfOqnTp3irbfeYuDAgbi7u9O5c2cCAgJwd3fn2rVrnDp1it9++41OnTqxYMECGVlOlE/XrjBjBixcqDaaO3kSPD2tHVWVFBYWRlZWFmFhYeh0OrKysvD09LRY1zB9kvX29q6U7mf68/n5+Rl1MSCEvTDqW1SnTh0WLVpESkoKy5cv55577uHq1av8/vvvgHq//fDhw+zbt08SujDOnDnQpInavS0iwtrRVFkFS76V0dBMURQyMzPtojGbNLwT9sCokrpetWrVePTRR3n00UfNHY+wEZX+A1a9Onz4IfTvDytXwuOPw/33V865hSgH/XeiZs2aMqiNsFnyyRTFskT3ozL16wf6e5yTJkFOTuWdWzg0e2sAKISpJKkL2/Luu+DvD2fPgoyBIMyktK5vQjgSSerCttSq9fdMbvPnQ1ycNaMR2Ma95IrGYGrXNyHsjSR1YXtGjoRHH1VncgsJgdu3rR1RleYISd3Yrm9C2CuTkvoff/zBt99+y6+//mrueIRQLV8OPj5w7BjMnm3taEQZbCHxCyFMSOqfffYZzZs3Z/jw4fTo0YMuXbrw559/WiI2UZX5+sLq1erfb78N+/dbNx5RKknqQtgGo5P67NmzeeKJJzhz5gw//PADADNnzjR7YEIwejSMGwc6HYwfDzduWDsiIYSwaUYn9QsXLhAZGUnz5s3p378/n3zyCRs3brREbELAsmVQvz78/jvIxaMQQpTK6KSel5dH9erVDY9btmyJTqcjNTXVrIEJAait4fXdkJYvhx07rBqOMB+dTkd2drZU2QthRiY1lPvoo4/45ZdfyM7OBsDZ2ZmbN2+aNTAhDAYOhOefV/8ODYWMDKuGI0qWnJzM0qVLSUtLK3NbfVI3ZoAjGURGiNIZndTvv/9+3nzzTXr16kWtWrW45557yMnJYc2aNezatYvr169bIk5R1S1YAM2awaVLMHWqtaMRJYiOjiY+Pp6YmBiLHd+eZpETorIZndT37NlDZmYm8fHxfPLJJ4wcOZLevXuzatUq+vfvT+3atbn33nstEauoymrUgPXr1SlaP/4YNm2ydkSiGKGhobRo0YLg4GCLHV+mTBWiZCZN6AJwzz33cM899/D4448b1iUkJHDo0CGOHj1qluCEKKRHD3j1VYiKgmeegZ491a5vwmYEBAQw1YI1KTKIjBClM+uIco0bN2bMmDHMmzfPnIcV4m+RkdCuHVy9qib2ypxwRphVcnIya9euLdf9dyFE+djEMLFRUVF07doVDw8P6tWrx4gRI4iPjy+0TU5ODpMnT6ZOnTrUrFmT0aNHy49BVeTmpla/u7jA11//3TK+CkpOTmbNmjV222gsOjqa8+fPExsbW+FjabVaatWqJVOiiirPJr4Be/bsYfLkyezfv5/t27dz584dBg0axI0Cg43MmDGDb7/9ltjYWPbs2UNycjKjRo2yYtTCatq1+3sGt6lT1RndqqDo6GjOnTtnVzOPpaWlsXTpUpKTkwkNDaVp06aMGTOmwseVpC6Eyia+Ad9//z0hISG0bt2a9u3bs27dOhITEzl8+DAAmZmZrFmzhkWLFtGvXz86d+5MdHQ0v/zyC/tLGD40NzeXrKysQotwIC+9BH37qqPMjR0LubnWjqjSPfjgg2RlZfHggw9aO5RCCibuu8XExBAfH8+aNWsICAhg4sSJ+Eq7CCGKdXcOyy3H75xNJPW7ZWZmAuDt7Q3A4cOHuXPnDgMGDDBs07JlSxo0aMC+ffuKPUZUVBReXl6GJSgoyPKBi8rj5KRWw9epA0eOwOuvWzuiSrd161Y8PT3ZunWrtUMppGDivltwcDAtWrSQ1utClENQUFChPBYVFVXmPia3fgfYuXMnO3fu5MqVK0VGhVq7dq1Jx9TpdEyfPp2ePXvSpk0bAFJTU3F1daVWrVqFtvX19S1xJLuIiAjCw8MNj7OysiSxO5r69dV76sOGwcKFMGAADBli7agqTWhoKKmpqTY3R3hwcDAxMTHFJm5fX1+mTp1KQEAAeXl5VohOCPuRlJSEp6en4bGbm1uZ+5ic1GfPns2cOXPo0qUL/v7+aDQaUw9VyOTJkzl58iT/93//V6HjuLm5lesNEHbukUdg8mRYsQImTIDjx6tMN7eAgADCwsJsrntXwcQthDCdp6dnoaReHiYn9dWrV7Nu3TqeeuopUw9RxAsvvMB3333H3r17CQwMNKz38/Pj9u3bZGRkFCqtp6Wl4efnZ7bzCzv1zjuwZw+cPAkhIbBlizpIjSiVfsjVGTNm0KBBA2uHYxT9VK8F56EQQlTgnvrt27f5xz/+YZYgFEXhhRde4KuvvuLHH3+kcePGhZ7v3LkzLi4u7Ny507AuPj6exMREevToYZYYhB2rVg02bgR3d/j+e3jvPWtHZBfWrFljt0OuyvztQhTP5KQ+adIkPv30U7MEMXnyZD755BM+/fRTPDw8SE1NJTU1lVu3bgHg5eVFWFgY4eHh7Nq1i8OHDxMaGkqPHj247777zBKDsHOtW8PixerfM2eqjedEqcLCwmTIVSEcjMnV7zk5OXzwwQfs2LGDdu3a4eLiUuj5RYsWlftYq1atAqBPnz6F1kdHRxMSEgLA4sWL0Wq1jB49mtzcXAYPHszKlStNDV84omeegW3bYPNmtZvb4cNQs6a1o7JZMuSqEI7H5KR+/PhxOnToAMDJkycLPWdso7nyTL3o7u7OihUrWLFihVHHFlWIRgNr1sChQ+qANFOmVOkR52yNflrW4OBgGjZsaO1whHBIJif1Xbt2mTMOIczD2xs++QT69YN16+CBB9Q52IXVrVmzxjAta7du3awdjhAOSZoIC8fTuzfMmaP+/fzzajc3YXVhYWEWnZbVFMY0uNP3FrDXsfZF1VChpJ6RkcHChQuZNGkSkyZNYtGiRYbR4ISwqogIdSCanBwYMwZkmGCr00/LakvDwhqT1O25t4CoOkxO6ocOHaJp06YsXryY9PR00tPTWbx4MU2bNuWItDwW1qbVqsPIBgaq99efflqmaRUVIr0FhD0wOanPmDGDYcOGcfHiRTZt2sSmTZtISEjg4YcfZvr06WYMUQgT+fjA55+Ds7P6rzSyFBUgvQWEPahQSf3VV1/F2fnvtnbOzs688sorHDp0yCzBCVFhPXqoI84BhIfDr79aNx4hhLAgk5O6p6cniYmJRdYnJSXh4eFRoaCEMKtp02D0aLhzBx57DNLTrR2RzdLpdGRnZ8tIbULYKZOTenBwMGFhYcTExJCUlERSUhIbN25k0qRJjB071pwxCispbV5su6Lvv96sGfzxhzrxiyStYklSF8K+mZzU3333XUaNGsX48eNp1KgRjRo1IiQkhEcffZT58+ebM0ZhJaXNi213vLwgNhbc3OC778AGP6MynrnlSbc04ehMTuqurq4sWbKEa9euERcXR1xcnKEFvEx56hiCg4Np0aKF47T27dABli9X/379dXVIWRsiSd28ikvg0dHRnDlzhmgZaVA4qAoPPlO9enXatm1L27ZtZRpEB+OQ82JPmvR397bHH4fz560dkTCCVqulVq1aaMsxtW5x/cpDQ0MN3dLKexwh7IlRw8SGh4czd+5catSoQXh4eKnbGjOhixCVatkyOHEC9u+HkSNh3z6oUcPaUYly0Cf1vLy8MrcNCwsjKyurUE2TvltaYGAgzs7O5TqOEPbEqKR+9OhR7ty5Y/i7JMZO6CJEpXJzgy+/hM6d1eQ+caI6H7t8bh2K9CsXVZFRSb3gJC4fffQRgYGBRaqvFEUhKSnJPNEJYSkBAfDFF9C3rzowTefO8Mor1o5KWIgx1fZC2DOTP+GNGzfm6tWrRdanp6fTuHHjCgUlRKXo2ROWLlX/joiAH36wbjzCYiSpi6rC5E94SXOgZ2dn4+7ubnJAQlSqZ56BsDC13/rjj8OFC9aOSAghTGb0fOr6BnIajYZZs2YVavGen5/PgQMH6NChg9kCFMKiNBq1m9uJE+oQsiNHws8/Q82a1o5MlIOUwIUozOikrm8gpygKJ06cwNXV1fCcq6sr7du356WXXjJfhEJYmrs7bNqk3lc/fhzGjVMfS6KwefqkLoRQGZ3U9Y3lQkNDWbJkCZ6enmYPSohKV78+fPWV2nDu66/VwWmioqwdlVXpB2+ZMmWKlISFsBMmf1Ojo6MloQvH0qOHOkY8wNtvw/r11o3HyvSDt8joa0LYD5OTelRUFGvXri2yfu3atTL2u7BfTz6pltJBHXnul1+sG48VhYWF0bJlS0JDQy1yfJ1OR05OjgyLK4QZmZzU33//fVq2bFlkfevWrVm9enWFghLCqubMURvM3b4NI0aoM7tVQZYevOXy5cusXbuWtLQ0o/eViVmEKJ7JST01NRV/f/8i6+vWrUtKSkqFghLWl5yczNKlS036wbV7Wi18/LE6Acyff8Ijj8D169aOyqaZkmTXrVtHQkICsbGxRp+vuHHdhRAVSOpBQUH8/PPPRdb//PPPMiyjA1izZg3x8fHExMQYtZ/DzDRWowZ88w34+and3Z54AvLzrR2VzTLl/ntISAi+vr7cuHGjzIuBuy8a9LcG7p5B0JTPn8N8ZoWgAkn96aefZvr06URHR/PHH3/wxx9/sHbtWmbMmMHTTz9tzhiFFYSFhdGiRQuCg4ON2s+hfiCDgtSW8O7u6hzsM2aos7uJIky5/+7v70/16tW5fPlyiSVu/efpP//5T6GSeUm3Bkz5/CmKQmZmpmN8ZkWVZ3SXNr2XX36Zv/76i+eff57bt2+jKArVqlXj1VdfZebMmeaMUVhBQEAAU6dOtXYY1tetm9oK/rHH1NndGjdWk7sopGCSzcrKIisrq1z7DRo0iIMHDxYpcevpk3RISAjXr18vcTshhMrkkrpGo2H+/Pn8+eef7N+/n+PHj5Oens6sWbNkljbhWMaMgQUL1L9ffFGd4U0Uy9gR3urUqcPkyZPLvGVnqUZ7VbrtiHBIJpfU9RITE/nrr7+4ffs2FwqMmz1s2LCKHloI2/HSS5CQAKtWqSPOBQSo/dpFpdDpdGRnZ5u9irxg25Fu3bqZ9dhCWIPJSf3ChQuMHDmSEydOoNFoDBO86Evp+dKoSDgSjUad0S0xEbZsgWHDYN8+aNbM2pFZXHJyMsuWLWPGjBk0aNDAKjFYKqmHhYWRkZFhdNsRIWyVydXv06ZNo3Hjxly5coXq1atz8uRJ9u7dS5cuXdi9e7cZQxTCRjg7w8aN0KkTXL0KQ4eq/zq46Ohou+0+VlbDOX3bEV9f30qOTAjLMDmp79u3jzlz5uDj44NWq8XJyYlevXoRFRUlDayE46pZUy2pN2wI586pJfabN60dlUWFhoYW233MHkjLdlHVmJzU8/Pz8fDwAMDHx8fQf7Rhw4bEx8ebJzohbJGfH2zdCrVqqVXwwcGQl2ftqCzG0iPLCSHMx+Sk3qZNG44dOwZA9+7dWbBgAT///DNz5syhSZMmZgtQCJvUqhV8++3ffdj/+U/pw24BlrqXbgqZu13YA5M/nW+88YbhizZnzhwSEhK4//772bp1K0uXLjVbgELYrF69ICZGHVY2Ohpee83aETkcSepCGMfk1u+DBw82/N2sWTPOnDlDeno6tWvXln7qouoYNgw++AAmTVKna/X1henTrR2VybRaLTVr1rRI4tJqtXh5eZV7YBohhPFM+ubeuXOH/v378/vvvxda7+3tLQldVD1hYTBvnvr3jBnw2WfWjacCLJ3Uy1vSNedwww41dLEQZTDpm+vi4sLx48fNHYsQ9mvmTND3+pgwAb7/3rrx2DlzJmJpAS+qEpMvx8eNG2e2fqt79+7lkUceISAgAI1Gw+bNmws9HxISgkajKbQMGTLELOcWwiw0Gli8GMaOhTt31PnY9+yxdlQORafTSXIWogwm31PPy8tj7dq17Nixg86dO1OjRo1Czy9atKjcx7px4wbt27dn4sSJjBo1qththgwZUmhaRzc3N9MCF8JStFr46CN17vXvvoOHH4adO9VJYUQRWq2WatWqlfuWnb7EXadOHQtHJoT9Mjmpnzx5kk6dOgFw9uzZQs8Ze1996NChDB06tNRt3Nzc8PPzMy5IISqbiwvExsJDD8GPP8LgwbB7N7Rvb+3IbI5Wq8Xd3V1akwthRiYn9V27dpkzjjLt3r2bevXqUbt2bfr168ebb75Z6hV7bm4uubm5hsfS4lZUGnd3dR72wYPhl19g4ED46Sdo0cLakQkh7MjdecvNza3MWmqjL5EvXLhgmLylsgwZMoT169ezc+dO5s+fz549exg6dGipk8ZERUXh5eVlWIKCgioxYlHl6YeT7dgR/vwT+vdXZ3kTQohyCgoKKpTHoqKiytzH6KR+zz338OeffxoeBwcHW3wu4scff5xhw4bRtm1bRowYwXfffcfBgwdLnTgmIiKCzMxMw5KUlGTRGIUoolYt2LYN7r0XLl+GAQPUf21YcnIyS5YsMQz7bCt0Oh05OTnSSE5UKUlJSYXyWERERJn7GJ3U7y6lb926lRs3bhh7mApp0qQJPj4+nDt3rsRt3Nzc8PT0LLQIUenq1oUdO6BJE7hwQU3sBS6Kbc2aNWtsckY2nU7HrVu3Kr2WUAhrujuHlaeBuF22ULl06RJ//fUX/v7+1g5FiLIFBKit4AMD4cwZGDQI0tOtHVWxwsLCKjwjm60OpyqD0IiqwOhvnb6f+N3rKiI7O5u4uDji4uIASEhIIC4ujsTERLKzs3n55ZfZv38/Fy9eZOfOnQwfPpxmzZoVGqpWCJvWqJFaYq9XD+Li1BK7DSZ2c8zIVt6kri99V1aS1fdzl9K+cGRGt35XFIWQkBBDNUBOTg7PPvtskX7qmzZtKvcxDx06RN++fQ2Pw8PDAZgwYQKrVq3i+PHjfPTRR2RkZBAQEMCgQYOYO3eu9FUX9qVFC7WbW79+cPSo2nhuxw6oov2upUpdCPMzOqlPmDCh0ONx48ZVOIg+ffqU+sXetm1bhc8hhE1o3Rp27YK+ff8usf8vsesbqc2YMYMGDRpYO9JKpa8a9/T0tGi1vUajsfg5hLAmo5N6wVHdhBAmaNVKTez9+qmJ/X8l9oKN1GbPnm3tKC2ipKp5fVI392Qy+qlbPTw8Sj2/EI5CPtlCWIM+sfv6wrFj0L8/YSNHVriRmq2zdFK9u0ueoijcuHFDqvhFlSFJXQhruffevxP78eMEjR/PjKeeKncjNUtOk2rrSro4sNUueUJUlqr3ayCELbn3XnVseD8/NCdO4PfEE3DlSrl2LZjUq1p3rZKSujm65AlhzySpC2FtLVvCrl0ofn64xsfj1L8/XLpU5m4FE5sjJnVTXlPBLnlarRZPT89Su9xqtVq8vLwq3C1XCFshSV0IW9CyJfk7d5IXEIDmzBm4/344f77UXRy90ZeiKFy/ft3k++FarZabN28SHR1daCjrghcL+m2WLVtmc0PjCmEKk2dpA9i5cyc7d+7kypUrRa6m165dW6HAhKhymjcnJSaGwNBQNOfOqYl9xw61UZ0wSXR0NOfPnyc2NpauXbsCRVvaV4VeB6LqMPkSf/bs2QwaNIidO3dy9epVrl27VmgR9k+j0eDl5eWwJUFblB8YSP6uXdCmDaSkwAMPwJEjFjmXPZT0KzrqXGhoKE2bNmXMmDElbiP34YUjMbmkvnr1atatW8dTTz1lzniEDbGHH32H5OcHe/bAkCFw8KA6UM2WLdCrl1lPo///rWw6nY6srKxyfa4URalQUg8ICGDYsGHExsZSv379Ys9pjqFxhbAVJv9a3759m3/84x/mjEUIoeftrVa9P/AAZGWpk8Bs327RU1ZWY7vKbtT3zTffEB8fLwNniSrB5KQ+adIkPv30U3PGIoQoyNMT/vtftcR+6xY89BB8/rnFTueoE54MGzaMFi1aEBoaau1QhLA4k6vfc3Jy+OCDD9ixYwft2rXDxcWl0POLFi2qcHBCVHnVq8PXX8OTT8IXX8Djj0NaGkyZYu3IKiwlJYXPPvuMZ555hsaNG5vtuPqaAG9vbwB8fHx44YUXCAgIIDU11WznEcIWmVxSP378OB06dECr1XLy5EmOHj1qWPRTqAohzMDVFTZuhOefB0WBqVPh9dfVv+3Y+vXruXjxIrGxsSVuc+XKFT777DOOHz/OZ599xpVyDMyjKAqZmZkO1WdfiPIyuaS+a9cuc8YhHERVnmnMopycYPly8PeHf/0L5s2D1FR4/31wrlDPVKsZP348f/zxR6kt07/88ksuXrzI7NmzURSF2NhYevbsadL5kpOTWbZsGVOmTJHGn8JhySdbmJWMvW1BGg288QZ8+CFotbB2LYwcCTdvWjsyk9SvX5+JEyfi6+tb4jajR4+mUaNGREZG0rhx41IvAMoSHR3NmTNnpMGccGgVusTPyMhgzZo1nD59GoBWrVoRFhaGl5eXWYIT9icsLIysrCzp82tJkyZB3brq/fXvvlPnZP/qK2tHZTStVku1atVKLTXXq1ePsWPH0q5dO8aOHUu9evVMPl9oaCjZ2dnSYE44NJNL6ocOHaJp06YsXryY9PR00tPTWbx4MU2bNuWIhQbLELZP+vxWkuHD1S5vtWrBvn043X8/zhcvmnQoRxw3vjjy2RRVgclJfcaMGQwbNoyLFy+yadMmNm3aREJCAg8//DDTp083Y4hCiGL17An/93/QoAGa33/Hf9Qo9bGRHDGpF/eadDod2dnZDvU6hbhbhUrqr776Ks4FGuk4OzvzyiuvcOjQIbMEJ4QoQ+vWcOAASpcuOF27htPgwWCm8SNsYURBfev31NRU/vrrL1asWFGuiVeKawEvSV1UBSZ/Wz09PUlMTCyyPikpCQ8PjwoFJYQwgp8f+Tt3cmPwYDS3b6t92ufMqXCXN2OSuqUuAPSt3zds2MAPP/zA2bNnCzXCdMRaBiEqwuRvYHBwMGFhYcTExJCUlERSUhIbN25k0qRJjB071pwxCiHKUr06f65ciS48XH0cGQkTJkBubqWc3lJJXd/6feDAgdy8eZP69esXaoRZVlJPSUlh7dq1haZeFcKRmdz6/d1330Wj0TB+/Hjy8vIAcHFx4bnnnuPtt982W4BCiHLSatHNn4+2eXOYPBk+/hgSE2HTJnUseTukb/2+fft20tLSqF+/vlEN3datW2eYerVjx44mx6GfhMbT01P6uAubZvKn09XVlSVLlnDt2jXi4uKIi4sztIB3c3MzZ4xCCGM884w6q5uHhzrbW48e8Pvv1o6qQp588kmT+qmHhISUOfVqeUg1v7AXFb7krF69Om3btqVt27ZUr17dHDEJG2ALjaREBQweDD//DEFBcPYsdOsG27ZZOyqT+fn5ldhPvbSE6+/vX+YAN0I4EqOq38PDw5k7dy41atQgXH/vrgQyoYt9s9Zc28KM2raFX3+FUaNg3z548EGYPx9efFEdnc6GJScns3TpUnr27MnmzZtp2LAh1apVQ1NM3PqW7vJ5FcLIpH706FHu3Llj+LskxX3xhBBW4OcHu3ap99jXrIGXX4a4OHWo2WrVrB1didasWUN8fDz79u3D2dmZDRs2MGzYMLPXHEm1unA0RiX1gpO4fPTRRwQGBhb5kimKQlJSknmiE0JUnJubmsQ7dIDp02HDBjhzBjZvhsBAw2b6fty2cBstLCyMjIwMQ0n9qaeeIjMz0+zn0c8hL4SjMPmyt3Hjxly9erXI+vT0dLPOjSyEKJt+drwSB2bRaOCFF2D7dqhTBw4fhi5d1Pvu/6NP6ooNTOkaEBDA1KlTadWqFWPHjsXPz6/INlLKFqIok5N6SV/87Oxs3N3dTQ5ICHtjC8ml3LPj9e0Lhw5Bu3aQlgZ9+6L58EOrzc2u0WhM7iamL2UX/C3SarV4eXnJLUBRZRndT13fQE6j0TBr1qxCVXX5+fkcOHCADh06mC1AIWydPqnXrFnTar0FjJodr1Ej+OUXCAmBL77A6fnn8Rk9mpRZs1i7di0hISE0adLE0iED5u9loT9ecVXqOp2OW7duScleODSjv0lHjx7l6NGjKIrCiRMnDI+PHj3KmTNnaN++PevWrbNAqEKIkpQ0A1mJtQg1asDnn0NUFIpWS80vv+STgQNJOHmS2NjYSok5JSWFNWvWkJycXPbtAxNptVpq1qyJRqORpC6qBKNL6vrGcqGhoSxZsgRPT0+zByWEMI9SaxE0Gpg5E12XLvD44zxz5Qpcu0ZwcHCl1DisW7eOc+fOER0djVarNdw+mD17ttnOodVq8fDwkJHgRJVh8jCx0dHRAJw6dYrExERu375d6Plhw4ZVLDIhRKVQ+vQh5bvvCJgxg9n796v92FNTYd48cDb5J6JMISEhXLlyhdDQUJydnct/+8BI+ir5nJwcsx9bCFtj8jc2ISGBESNGcOLECTQajaGxir6BSn5+vnkiFEJYXL6fH7e2buX688/jv3EjvPMOHDgAGzdC3boWOae/vz9hYWEEBATg7Oxc7O0Dc9JqtVSrVo20tDRWr17NjBkzLHo+IazB5PqoqVOn0rhxY65cuUL16tX57bff2Lt3L126dGH37t1mDFEIUSlcXPhj2jRub9igjhu/dy906IB2+3abGzLYlB4H+qS+fv368vUUEMIOmfwt3bdvH3PmzMHHxwetVotWq6VXr15ERUUxdepUc8YohKhEulGj4OBBaNMGrlxB++CD1Jo3D+3/ZmO0BcV1ZyuvkJAQWrZsWWxVv/5iIS8vz+rdFIUwhclJPT8/Hw8PDwB8fHwMrVYbNmxIfHy8eaITQlhHixbquPHPPac+fucd6NULzp+vtBDS0tJYtWoVf/31V7m2L6v0npycbCidl1TVL0ld2DuTk3qbNm04duwYAN27d2fBggX8/PPPzJkzp9L6uAohLKhaNVi5Er78EmrVUkvvHTvCZ59VyuljYmI4e/YsP/zwQ7m2LyupR0dHc+7cuSJdbjUaDV5eXjZ1e0EIU5n8KX7jjTcMX545c+aQkJDA/fffz9atW1m6dKlRx9q7dy+PPPIIAQEBaDQaNm/eXOh5RVGYNWsW/v7+VKtWjQEDBvC7nc8PLYSpKn0Eu1Gj4Ngx6NkTrl+HJ56AiRPhxg2TDqfvk56SklLqdsHBwTRv3pxBgwaV67j6YW5Lel9CQ0Np1qwZISEhhdbLNMPCkZj8KR48eDCjRo0CoFmzZpw5c4arV69y5coV+vXrZ9Sxbty4Qfv27VmxYkWxzy9YsIClS5eyevVqDhw4QI0aNRg8eLB0URFVklWGpW3QAHbvhn/9S+3fHh0NnTurM74ZST+kbVmDVPn6+vLcc89Rp06dch23rKQeEBBAWFgY/v7+xoZsscFxhDA3s16aent7mzTm8tChQ3nzzTcZOXJkkecUReG9997jjTfeYPjw4bRr147169eTnJxcpEQvhD2z+RKjszPMmQM//gj160N8PHTrps7RbkQX1rCwMFq2bFmkxGzLyj22/v/YwnwAomoy+dejX79+xY78dO3aNaNL6qVJSEggNTWVAQMGGNZ5eXnRvXt39u3bV+J+ubm5ZGVlFVqEsGU2n9T1+vRRS+jDh8OdOzBzJvTuDRculGt3/ZC2ppSYyyMlJYWlS5eSlpZmtmPqL0TKOzhOcUldEr0w1t05LDc3t8x9TP712L17N8uXL2fEiBHcKHBv7fbt2+zZs8fUwxaRmpoKqFVxBfn6+hqeK05UVBReXl6GJSgoyGwxCVHl+fjAV1/B2rVqn/aff1ZnfjPzjG8ajYZq1aoZdaGzbt064uPjiYmJMduFUklj6xtDkrowVlBQUKE8FhUVVeY+Ffqk79ixg9TUVO677z4uXrxYkUOZXUREBJmZmYYlKSnJ2iEJYRPuTi7lbbhWhEYDoaFw/Dg88IDacO6f/4RHHlGHmTUhlrvpB4wxJimPHz+exo0bM2bMGPup/RCiGElJSYXyWERERJn7VOiT7u/vz549e2jbti1du3a1yEhyfn5+AEWq0tLS0gzPFcfNzQ1PT89CixCiaCLV3y9etmyZYdY0ozRqBLt2wbvvgqsrbNmiDlzz5Zdl7qooCllZWSaXXtPS0li+fHmh3wdfX1+GDBnC559/bnLDNp1OR1ZWFunp6VKyFlZzdw5zc3Mrcx+Tk7q+QZybmxuffvop06ZNY8iQIaxcudLUQxarcePG+Pn5sXPnTsO6rKwsDhw4QI8ePcx6LiGqIv39YsAwa5rRtFp1IpjDh6FDB/jrL3j0UbX729Wr5g24gNjYWENVe0HffPMN8fHxJg8Fq9PpuH79ulSXC7tjclK/e3jGN954gw0bNrBw4UKjj5WdnU1cXBxx/+sek5CQQFxcHImJiWg0GqZPn86bb77JN998w4kTJxg/fjwBAQGMGDHC1PCFEP+jv188ZcoUmjVrRmhoqOkHa9NGnQjmtdfURP/ZZ9CqlTp3ewXutRecF72gMWPG0KJFC4KDgwutHzZsGC1atCixYZvc3xaOyuSknpCQgI+PT6F1o0ePZv/+/axdu9aoYx06dIiOHTvSsWNHAMLDw+nYsSOzZs0C4JVXXmHKlCn885//pGvXrmRnZ/P999/j7u5uavhCiLsUnDWtQlxd4a23YP9+Ncn/+ScEB6uD2Bh73/5/9Em94L3xtLQ0YmNj6du3LzExMYXaBPj4+PDCCy8Uei36YWJTUlJQFIXMzExJ6sLhmJzUGzZsWGzjkzZt2jBhwgSjjtWnTx8URSmy6Aen0Gg0zJkzh9TUVHJyctixYwfNmzc3NXQhRGXo2lWtjv/3v9U+7ps3Q6tWaD76yCwt5GNiYoiPj+ftt98mPj6+zMFsihsmVt9I8NKlS1JyFw7BqPnUw8PDmTt3LjVq1CA8PLzUbRctWlShwIQQDsDVFSIjYeRIdWjZw4dxmjQJ3wceQLdqVZm7azQaatasibOzM7Vr1wbUUrtOpyM4OJiNGzfSt29fdu/eTUhISKmTv4SGhpKamlpo0Jvo6GjDoDJhYWFS+yfsnlFJ/ejRo9y5c8fwd0lMGVVOCOHA2rVTq+MXLkSJjKTa3r0oXbtSPySE288/X+Ju+mp3Z2dnvL29Det1Oh2+vr688MIL/zt8O3x9fUtN6gWHidW3lg8NDSU7O7ti7QjuIvfrhTUZldR37doFwJ07d9BqtaxevZp77rnHIoEJISpXwfvW+i5dnp6e5uvj7ewMr75K/kMPkRcWhvuvv9Jo5Upyd+yA999Xu8ZVMn0jQT8/v1IHszKGfq53IazBpG+ri4sLx48fN3csQggzM2YikruTusVKmy1bkrpxIzkrV3LHywu3s2dx7dePOq+/Dteumf98QlQhJl+Cjxs3zuQ+oEI4EluewcvYiUgqjUaDbsIE4jdvJmvUKDSKgsenn+LUpg18+qlZh5oVoioxqvq9oLy8PNauXcuOHTvo3LkzNWrUKPS8NJSrmqrisJwFE2dxkxxVhpLe97CwMLKysso9EUll0mq1VAsK4s8FC3B79lk0zz2H6/nz8OSTaNeuxeXVV9E1aoSXl1eV+jwJUREmJ/WTJ0/SqVMnAM6ePVvoOWkoV3Xpk0tVYguJs6T33RwTkZRHRe/BK/ffT8qWLQR9/jlO8+ah3bmToL17yQgJwePttystqScnJ7Ns2TKeL6XxnrlZpP2CqLJM/gTt2rWrxOXHH380Z4xC2LTKSpy2zCz34N3cUCIi4ORJdA8+iObOHWp/+CFOrVvDxx9DgWMnJycXmV41JSXFMLiMqfRd3Mrq825O0lpemJNcFgohbCuxNG2K7uuvSfnwQ243bIgmNRXGj4devdTBbFBvecTHxxMbG2vYbdmyZezYsYMVK1aYfOrQ0FBatmxJ165dmTJlCr/99luFX44QlanCSf3UqVN8//33fPPNN4UWIYT90M+WVnBOB2sn+pt9+5K0dSv58+ZBjRqwb586St0//0nYiBG0aNGCMWPGFNnvxo0brF27tsjMjsUNNXs3fa3LokWLSEpK4u233y5XrOV5r2y5QaVwHCYn9QsXLtC+fXvatGnDQw89xIgRIxgxYgQjR45k5MiR5oxRCLth7URoTqWNj15pr9PNDeXll+HsWXjySbVV/IcfEtSnD7Pc3fEr0I5gypQpDBgwgBo1anD+/PlCpXgoX1LXmzdvHn5+fjRu3Jjk5OQyG4CW5/2w2Z4IwqGYnNSnTZtG48aNuXLlCtWrV+e3335j7969dOnSxSLzqgthDyoj2dlCia/SL14CAuCTT+Cnn6BjRzRZWdRZsIBGQ4bg8e23oNMZJqSZPHkyTZs2LbYUX14dOnSgW7duXLt2zaQknJaWxtKlSwv9H+mnuLXFngjCcZic1Pft28ecOXPw8fFBq9Wi1Wrp1asXUVFRTJ061ZwxCiEKqNIlvl694NAh8teuJc/XF5fkZPxfeon6jz6K9uefAXW2uYkTJ+Lr61vuwxZ3oVRw+lZjL2L0k80U/D+SBpWiMpic1PPz8/Hw8ADUaQ71X4aGDRsSHx9vnuiEEEVUtMRnCyX9CtFqUZ56isQdO7g6Ywa66tVxP36c6kOG0HzmTLTnzhl9SH2r9+joaMO64qZvLa/g4OBS53MXwlJMTupt2rTh2LFjAHTv3p0FCxbw888/M2fOHJo0aWK2AIUQhVW0xGcPJX2tVouXl1epY14o1aqR/uyzJGzfTuYTT6BotXjv2YNHjx40WbQIpz//LPHYt27dYtmyZYYLG32rd3NN7OLr68vUqVOlVC4qnclJ/Y033jBURc2ZM4eEhATuv/9+tm7dytKlS80WoBDCvEoq6Ws0mlJHb7PkffS7G6IZMzJhvo8PV+fM4eaBA1z7xz/Q5OXhv2kTjQcORFvMePJarZbY2Fji4+MNJXNHrBp3pEabovxMTuqDBw9m1KhRADRr1owzZ85w9epVrly5Qr9+/cwWoBDCvEpKYOZo4W2MgqVxcwwvrLRsSfzChWR/8w1ZrVujvXUL7YIF0LgxzJsH2dmGbfUXNuYomdtq8rTVuIRlmXXwGW9vbxkiVghRLlqtlps3bxaqBjeFTqcjOzvbkLzy7r+fE6tXc3nVKpQ2bSAzE15/HacWLfCIjobcXLOWzPVTrer7+Bc32l1p7L6Ng7ApRo39Hh4eXu5tZUIXIURZzDEZjqIohZI6ABoNN/r1I3/iRJy/+AJmzUJz/jx15sxBiY5G98Yb0LevmV5FYfrR7mJiYujWrVu5trf2hEDCcRiV1I8ePVqu7aS0LoR905c2g4ODadiwocXOY/HJcLRaeOIJGDMG3Zo1MHcu2qQknJ55hvoNGkBEBJj5dmFYWBgZGRkEBweXe3trTwgkHIdRSX3Xrl2WikOIKslWZ+gytrSpr0KeMWMGDRo0KPd5KlINrm/Yd+2uhnDFcnFB++yzMGECrFyJMn8+LomJ8NxzBNavjzJzJgwcaHQMxfHz82PixIlUr169XNs7YiM9YT0V+hXJyMhg4cKFTJo0iUmTJrF48WIyMzPNFZsQDs9WGzOFhYURGBhIdnZ2ue71GtNNzhyN4vTHMXqu9WrV4MUXyf/9d9LfeAPFzw/ny5dxmTKFwL59cVq9Gk1uboXi0t/jLziOvhCVxeRv1aFDh2jatCmLFy8mPT2d9PR0Fi1aRNOmTTly5Ig5YxRCVLKAgABq1KjB5cuXy5WojRkQp7xJvTzblac/e7Fq1CArLIzcU6f4KzISxd8f55QUXGbMoP4DD+CyciXaAsm9rMZsxQ0LK4Q1mJzUZ8yYwbBhw7h48SKbNm1i06ZNJCQk8PDDDzN9+nQzhiiEsKSSEqMxo6JZogrZmKRecBv9xC3lSvTVqnE9JERN7nPnogQG4nzlCm6vvkqXMWPwXr0a/jf+e2k1EQWHhdXpdOTk5BSpfbm7lb4QllChkvqrr76Ks/Pft+WdnZ155ZVXOHTokFmCE0JYXknJs7RR0Wy5G5ZWq8XDw4PatWuXvzagenWujxtH7m+/cXXePHQNG+KSno7P4sU4NWnC08nJ3Fu/fokXOAUvgHQ6Hbdu3SpS/S5JXVQGk5O6p6cniYmJRdYnJSUZxoQXQjgmWxpq9u4pVbVaLd7e3nh7extfxe/qSvbYsdw+cYK0hQvJbdkSTXY2gf/5Dws/+oj6r79O1i+/FEnMtjgsrC1feAnLMTmpBwcHExYWRkxMDElJSSQlJbFx40YmTZrE2LFjzRmjEMLCjE0AtjSNaHFJvcIN8VxcyB4+nEvffkv+d9+h69MHTV4eTp98gmfPnvDII7B3rzq/u42ypQsvUXlM/tS/++67jBo1ivHjx9OoUSMaNWpESEgIjz76KPPnzzdnjEIICzM2AdhiN6y7k7ux+9aqVQtnZ2fDBYFGo8GrVi00Q4ei276d5M2byR81CkWrRbt1K/TuDffdhyYmBu7cscArqhhbuvASlcfkpO7q6sqSJUu4du0acXFxxMXFkZ6ezuLFi3FzczNnjELYDVur8ixvqdUREoC5k/rd793t9u25s2EDl3fuRPfMM+DuDr/+itO4cTTs04daK1dCCTPDVTQ+U9jihZewPKM/XTqdjvnz59OzZ0+6du3KnDlzaNasGW3bti33YAtCOCpbq/Isb1KXBFB+eY0aoVu+HP74A2bPVvu6p6VRZ9EinBo3xu255/C8cKHIfiUldXP12xcCTEjqb731Fq+99ho1a9akfv36LFmyhMmTJ1siNiHsjiOUeO2BTSTCevVg1izyz58n7d13yWnbFk1uLi6ffMID06bRYPx4NF99BXl5pR7GJl6LcBhGf4rWr1/PypUr2bZtG5s3b+bbb79lw4YN0k1DCKTEW1lsKhG6upI9YgSXN20ib+9eLg4Zwr81GjIPHcLpscegWTOIioJyztomREUY/Y1ITEzkwQcfNDweMGAAGo3GZu4hClHZbCrBWJCttRe4mzn/H0w6lkYDPXrwfocObGvXjvc6dUKpU0etpn/tNZwaNaLu5MlofvwRpBAkLMToT39eXh7u7u6F1rm4uHDHBlt/ClEZ7C2pmxqvrbUX0NNfbKSmppb7dZX1HhjzHt09f/r48eOp36wZY1avRvnjD1i3Dnr0QJOXR42tW3EaPBhatoSFC+HqVaNeqxBlMfpXSFEUQkJCGDVqlGHJycnh2WefLbROCFE2U0q/Fb2IMHV/W20vYMrFhrHvQWnbF5zRTr+ts7OzOvNejRrqzHC//ELeoUNkjRuH4uEBv/8OL70E9evDuHHw00823edd2A+jfxUmTJhAvXr18PLyMizjxo0jICCg0DohRNkqIyGZi622F6iMi43S3vOwsDBatGhhmD99/fr1JCUlsX79+sIbtm9P+ty55CcmwgcfQKdOcPs2bNgADzwA994L8+eDjd7eEPbBqPnUAaKjoy0RhxBVUlhYGFlZWTZX+gV1vvLyzvNuzVsQ1r7YuHv+9JCQEP744w9CQkIKbVdozvmnn4ann4ZDh2D1avjsM4iPh5kz4bXXYMgQCA1VR64zYtwPnU5HVlYWnp6e5nyJwo7Yx01AIRyUtRNSSYwZP12/vT21KzBVcbdLdDodN27cwMPDA61Wi7+/P0888QT+/v6F9i22VqZLF/jPfyAlBd0HH5DXvbvaiG7rVhgzBgICYOpUOHq01Lh0Oh0ZGRnk5eWRkZEhvZGqMLv4Bv773/9Go9EUWlq2bGntsISwCv0PuCV/uB0xSZvjfSvpdkl53q9SbxN4eqILDeXSxo3knTypltgDAiA9HZYtU6vqO3SAxYshJaVcr80R/w9F2ezmf7t169akpKQYlv/7v/+zdkhCWEVlJPWy2GPCMMf7VpH79+WulWnRQu3Xnpj4d4nd1RWOHYPwcAgMhAEDIDoaMjNLPIw9/h+JirOb/21nZ2f8/PwMi4+Pj7VDEqLKsteEUdHEXqm3S5ycYOhQ+PxztfHc8uXQo4daPb9zJ0ycCL6+MGYMms2bITfX8jEJm2c338jff/+dgIAAmjRpwpNPPlnsXO4F5ebmkpWVVWgRQlRtiqKQmZlpf/ec69SByZPhl1/g/Hl48021tXxuLnzxBU5jxtCga1ecn3sO9337ID/f2hELM7g7h+WW48LNLpJ69+7dWbduHd9//z2rVq0iISGB+++/n+vXr5e4T1RUVKEudkFBQZUYsRCiPOy1xF+alJQUPv30U1KKufdtFk2awOuvw2+/qQ3oXnoJpX59tNev47xuHX5PPIFTw4bw/PPw449ljj0vbFdQUFChPBYVFVXmPnbxTRo6dChjxoyhXbt2DB48mK1bt5KRkcHnn39e4j4RERFkZmYalqSkpEqMWAhRUEnV3vaQ1I2tsl+3bh0JCQmsW7fOsoFpNGrjuXfeIf/CBVI/+4y8iRPJ9/JCk5YGq1ZB//5qg7tnnoEdOyTB25mkpKRCeSwiIqLMfWz3m1SKWrVq0bx5c86dO1fiNm5ubnh6ehZahBCls1QjPFto3GcqY2MPCQmhcePGRfqpW5RWS85995G3YgVJBw+Sv2ULhIWBt7c6x/sHH8DAgeDnp/aP37YNZGhvm3d3DnMrx5gFdpnUs7OzOX/+fJF+oEKIijFn8rXnRG6M5ORk1qxZY+i7Xr9+fcLCwqhfv751AnJxQRk0SO3/npoKP/wA//wn+PjAX3+p64cMURvZjR8PX34JpdzKFPbFLpL6Sy+9xJ49e7h48SK//PILI0eOxMnJibFjx1o7NCEcgiWqwS2d1G3loiE6Oppz584ZRtvUarVUq1bNKrcUigyO4+KiltDff1/t375zJzz7rDoX/LVr8PHH8OijasJ/8EF1dDsZptau2UVSv3TpEmPHjqVFixY89thj1KlTh/3791O3bl1rhyZEpdPpdGRnZ5s1mdnDve27GZvU755NzVxCQ0Np1qwZoaGhgPpe1qxZ0yrvZXR0dMlzCTg7Q79+6r325GTYswdefFGd7/32bfjvf+G559RJZrp1U1vYnzghE83YGbv4Bm/cuJHk5GRyc3O5dOkSGzdupGnTptYOS4giKiM5WiKpV4ROpyM9PZ309HSbiak4d8+mVhJjLxYCAgIICwsz9F03NalXZL56/b4PPvhg+QbHcXJSJ5F59104exZOnVIHvOnRQ22Ad/Ag/Otf0K4dNG0KU6aoA+HcvGl0bKJy2UVSF8Je2GOJt6J0Oh3Xrl3j2rVrNp3U755NrSTWqtavyHz1+n23bt1q/OA4Go3a533mTLUffHIyfPihOpmMuzskJKgD3zz0kNrwbvBgeO89dQIaKcXbnKrzyyOEDTLlIuDuhllVXXnfw4CAAIKDg4mJibHJ964iQ9Dq99XfAqgQPz+YNAm++QauXoXNm9UucQ0aqIPd/PADzJgBLVuqpfjJk+Hbb+HGjYqfW1SYJHUhrMiUpH53wyxzqkgVcHlYohRszHsYExNDfHy8SaVhS/Pz8yMsLAw/Pz+j97XY8LU1asDw4WoDuosX1Wr6hQvVseddXdVS/MqVMGyYWoofOBAWLIDDh2VUOyuRpC6Enbm7YZY5lVYFbI5bC9ZusR4cHEyLFi1Mnr++uPegvPfQy3r/KtJWolLaWeir6cPDYft2tXvcN9+oI9c1bqw2ttuxA159VZ1Stl49tWX9qlXw++9SVV9JJKkLYWfubphlTqVVARub1AuW+i1RA2DKMX19fZk6dapR713B81gyqVeEVRpP1qyp3ndfsUIdj/7MGViyRC21e3io08Z++aWa9Js3h4YNITQUNmwodvpYYR6S1IUQBuasxi1Y6tf/bc5bBhVpWGaL5zG3Sq0V0WjUKWOnToWvv1YT+r59MHcu9OmjVtUnJcG6dTBunDp0bevW6v34mBhJ8mYkSV0IO2NrDeVK6v9dsNRv1oZcxRzfnO6uATD2PLbSA8KqtzqcneG+++CNN2DXLnWgm23b4JVXoFMn9SLg1Cn1fvzjj6tJvnlztYHe+vXq/XuprjeJJHUh7IwlG8qZoqT+3wVL/ZZoyGWpxmF3l8zLOo8+iQNkZGQAFJvUrd2ewKqqV4dBg2D+fLUR3Z9/qlXz06ZBx45qkv/9d1izBiZMUO/RN2yoluo/+ECt2pckXy7O1g5ACGGc0NBQUlNTLdJQzhRhYWFcunSJGzdukJycTIMGDawdUoWEhYWRlZVldMk8Ly+PjIyMEu+v65O6tUabsyl16sCoUeoCkJGh9pHfswf27oVDh9Tq+g0b1AXUoWx79Ph76dpVbZ0vCqninywh7I8lG8qZIiAggJo1a3Lp0iW7u+9cHEvVAFiyu6D+lkxaWpqhlsDS3RPNqlYtdez5+fPVe/EZGWpL+lmz1Hvy7u5qn/lvv4XXXoO+fcHLS63Kf+EFNfFfuCCleSSpCyHMoKJdxUxVUpV2ceu1Wi1eXl5oNBqLxFJWErVkgzv9LZl169YZkrq9NvAD1BJ4//4we7Z6T15fkl+4UO0mV7++2g/+6FG19f24cepAOH5+MGKEenGwZ0+VnH1OkroQosLK21XM3KVHRVHIzMwsd1K3ZAO2spJoWQ3uUlNT+fjjj0lNTTX63MWNXWCphoSmqHB7Ajc3tco9PBxiY+HSJUhMVFvOT58O3burM9JduaK2vp85Uy3he3mpfevHj4elS9ULAwcfv17uqQshDCyd+ArOIjZ79myLnKM0lnx9Zd2LL6taPzo6mgsXLhAdHW30e1PcLRmLjTJnAou0JwgKUpfHHlMf5+TAkSNq9f2+feqkNImJaiO7M2fUaWZBncymdWt1gBz90q6deuHgACSpCyEMCrbktoTQ0FCys7PNUnrUV6dnZWUZtY+lXl9Fk6ixDfQKsuZ0rzbD3R3+8Q910btyRW10p18OHoTUVDh+XF3WrlW3c3FRE3uHDn8v7dqBp6cVXkjFSFIXQlQac5Ye9Qk6OzvbDJFVjtKSbkXem+KSuq30l7eqevXUBngPPqg+VhR1Fjp9gtf/m56udrU7fLjw/k2aFE707durtQMWapdhDpLUhbAztlYqK6nEXDCpVHbfbP29+xkzZhjdxc7Sw7lWFkvXutgljUZtZFe/vjpRDaiJ/uJFNcEfOwZxceq/ly6pLeovXIBNm/4+Ru3ahZN8mzbqffvq1Sv/9RRDkroQdsYWk3pxJea7k4olEmVJybtgozVj70/fHXd5k7yUjIvS6XRkZWVR3UYS3t308Xk2bIi2cWMYM+bvJ69eLZzk4+LUUfCuXVNb5O/a9fe2Go1aqm/dWk3yrVurS8uWlX6vXpK6EHbGHpNHwQFazKmk5F2R+9N3K2+J15olY1v9TOgbyKWnp7Ns2TKTak4s6dKlSyxevLj4uHx81G51/fv/vS4nR03sBZP9yZPqjHXnz6vLN9/8vb2TEzRr9nei1/97zz3qfXwLkKQuhJ2xxWpVayWVkpK3LbX8rgy2+Jko2I3N2r0eSmJ0jY67uzrgTadOf69TFLVB3m+/qcvJk3//m5kJ8fHq8uWXf+/j7Kz2q2/ZsuhSwf9HSepCiAqzVlKxp+Rdkfv8lcFQFe3paZaLM51Ox9mzZ4mJieHxxx83W68HczJLjY5GA76+6tKv39/r9Y3y9Elen+hPnYLs7L+T/ddfFz6en59aujeRJHUhhKgEFbnPXxks0Zc8JiaG+Ph4tm7dapMXX4GBgURGRuJpia5rBRvlDR7893pFgcuX/+4/X3C5fFntcmfCAER6ktSFEJXGVu/9VgZz3ue3BEvUJAQHBxMTE2Mzkw/dzSo1TBoNBAaqy4ABhZ/LyoKzZ9VBdJ55xqTDV71vlhDCaiorqdvixUN5p3C1VszmHis+OTmZmJgYgoODba6EbqxKmzbX01Md4e7xx00+hO184oUQwkysnSBNYe2YzT1W/Jo1a4iPjycmJsYsxzOXS5cuERERwaVLl8q9T6UldTOwn0+8EELcxdqJ0JGYu9FhWFgYLVq0IDg42Kb+nz788EPi4uL48MMPrR2KRVj/HRZCCBPZUrKwd+Z+L/UXCc2bN8fZ2dlm/p9MqZGwp7nprf8OCyFEFWDrFyCWiM8WX7O+xXtgYGC59ymtvYGtVc1L63chhKgElmppbe7+5eZiiwkdTPt/KK3ngkWmla0A60cghBDCZJcuXWL27NlGNfyqDLaa1E1RWnsDW6uat/93WwghqjBzd0UTxrG191+q34UQwo7Z+qA2jqC0Wgdbe/81iqIo1g6iMmRlZeHl5UVmZqZlhgQUQggryMvL49KlSwQGBuLsLOW0ymaJ978i+Uqq34UQwo450r1rUXFyWSeEEHbMFqddrUps7aJKkroQQghhIlu7qLKNSwshhBBCVJgkdSGEEMJBSFIXQgghHIRdJfUVK1bQqFEj3N3d6d69O7/++qu1QxJCCCFsht0k9ZiYGMLDw4mMjOTIkSO0b9+ewYMHc+XKFWuHJoQQQtgEuxl8pnv37nTt2pXly5cD6iD6QUFBTJkyhZkzZxbZPjc3l9zcXMPjzMxMGjRoQFJSkgw+I4QQwmZlZWURFBREYmIiXl5ehvVubm64ubmVvrNiB3JzcxUnJyflq6++KrR+/PjxyrBhw4rdJzIyUgFkkUUWWWSRxSGWyMjIMvOlXfRTv3r1Kvn5+fj6+hZa7+vry5kzZ4rdJyIigvDwcMPjjIwMGjZsWOTKR5Ssa9euHDx40NphlJu147X0+c19fHMcryLHMHZfY7bXl3SkZq78rP39MZa147Xk+RVFoVOnThw+fLjQoDZlltJx4MFnSqqm8PLyki95OTk5OdnVe2XteC19fnMf3xzHq8gxjN3XlHN5enra1WfYmqz9/TGWteO19Pnd3NxMGtTGLhrK+fj44OTkRFpaWqH1aWlp+Pn5WSkqxzd58mRrh2AUa8dr6fOb+/jmOF5FjmHsvtb+/3V09vb+WjteW/2+21VDuW7durFs2TJAbSjXoEEDXnjhhWIbyt1NZmkTouqQ77uoquym+j08PJwJEybQpUsXunXrxnvvvceNGzcIDQ0t1/5ubm5ERkaW656EEMK+yfddVFV2U1IHWL58Oe+88w6pqal06NCBpUuX0r17d2uHJYQQQtgEu0rqQgghhCiZXTSUE0IIIUTZJKkLIYQQDkKSuhBCCOEgJKkLIYQQDqJKJvWkpCT69OlDq1ataNeuHbGxsYbnRo4cSe3atXn00UetGKEQwlxK+r6X9jsghL2qkq3fU1JSSEtLo0OHDqSmptK5c2fOnj1LjRo12L17N9evX+ejjz7iiy++sHaoQogKKun7npWVVeLvgBD2qkqW1P39/enQoQMAfn5++Pj4kJ6eDkCfPn3w8PCwYnRCCHMq6fte2u+AEPbKLpP63r17eeSRRwgICECj0bB58+Yi26xYsYJGjRrh7u5O9+7d+fXXX4s91uHDh8nPzycoKMjCUQshTFEZ33f5HRCOwi6T+o0bN2jfvj0rVqwo9vmYmBjCw8OJjIzkyJEjtG/fnsGDB3PlypVC26WnpzN+/Hg++OCDyghbCGECS3/f5XdAOJQyZ1y3cYDy1VdfFVrXrVs3ZfLkyYbH+fn5SkBAgBIVFWVYl5OTo9x///3K+vXrixxz165dyujRoy0WsxDCNOb+vpf2OyCEPbLLknppbt++zeHDhxkwYIBhnVarZcCAAezbtw9QJ6APCQmhX79+PPXUU9YKVQhRQRX5vsvvgHBEDpfUr169Sn5+Pr6+voXW+/r6kpqaCsDPP/9MTEwMmzdvpkOHDnTo0IETJ04AMGDAAMaMGcPWrVsJDAw0/DAIIWxPRb7vpf0OCGGv7GbqVXPq1asXOp2u2Od27NhRydEIISyptO97SeuFsFcOV1L38fHBycmJtLS0QuvT0tLw8/OzUlRCCEuQ77sQhTlcUnd1daVz587s3LnTsE6n07Fz50569OhhxciEEOYm33chCrPL6vfs7GzOnTtneJyQkEBcXBze3t40aNCA8PBwJkyYQJcuXejWrRvvvfceN27cIDQ01IpRCyFMId93IYxg5db3Jtm1a5cCFFkmTJhg2GbZsmVKgwYNFFdXV6Vbt27K/v37rRewEMJk8n0Xovyq5NjvQgghhCNyuHvqQgghRFUlSV0IIYRwEJLUhRBCCAchSV0IIYRwEJLUhRBCCAchSV0IIYRwEJLUhRBCCAchSV0IIYRwEJLUhRBCCAchSV0IIYRwEJLUhRBCCAchSV0IIWxAUlISffr0oVWrVrRr147Y2FhrhyTskEzoIoQQNiAlJYW0tDQ6dOhAamoqnTt35uzZs9SoUcPaoQk7IiV1IYrRp08fpk+fbu0wbJ6jvE9//fUX9erV4+LFi1aLwd/fnw4dOgDg5+eHj48P6enphucff/xxFi5caKXohL2QpC4cTmpqKlOmTKFJkya4ubkRFBTEI488ws6dO60dmlmEhISg0WjQaDS4uLjg6+vLwIEDWbt2LTqdrlJj2bRpE3PnzjU8NmeS7927NxqNhs8++6zQ+mXLlhEQEGCWc+i99dZbDB8+nEaNGpn1uKY6fPgw+fn5BAUFGda98cYbvPXWW2RmZloxMmHrJKkLh3Lx4kU6d+7Mjz/+yDvvvMOJEyf4/vvv6du3L5MnT7Z2eGYzZMgQUlJSuHjxIv/973/p27cv06ZN4+GHHyYvL6/S4vD29sbDw8Psx1UUhaNHj+Lv78+XX35Z6LnDhw/TqVMns53r5s2brFmzhrCwMLMdsyQdOnSgTZs2RZbk5GTDNunp6YwfP54PPvig0L5t2rShadOmfPLJJxaPU9gxRQgHMnToUKV+/fpKdnZ2keeuXbumKIqi5OTkKFOmTFHq1q2ruLm5KT179lR+/fXXQtv27t1bmTZtmuFxw4YNlcWLFxfapn379kpkZKRh+xdeeEGZNm2aUqtWLaVevXrKBx98oGRnZyshISFKzZo1laZNmypbt24tcp4pU6YoL7/8slK7dm3F19fXcMySTJgwQRk+fHiR9Tt37lQA5cMPPzSsy8/PV+bNm6c0atRIcXd3V9q1a6fExsYadf7Y2FilTZs2iru7u+Lt7a3079/f8P4WfJ8mTJigAIWWhIQE5aOPPlK8vb2VnJycQscdPny4Mm7cuGJfY3x8vAIoK1asUKpXr67cuHHD8Fzbtm2VWbNmlfoeGSM2NlapW7dukfX5+fnK/PnzlaZNmyqurq5KUFCQ8uabbyqKYvr/d1lycnKU+++/X1m/fn2xz8+ePVvp1auX8S9SVBlSUhcOIz09ne+//57JkycX27ioVq1aALzyyit8+eWXfPTRRxw5coRmzZoxePDgQvcvTfHRRx/h4+PDr7/+ypQpU3juuecYM2YM//jHPzhy5AiDBg3iqaee4ubNm0X2q1GjBgcOHGDBggXMmTOH7du3G33+fv360b59ezZt2mRYFxUVxfr161m9ejW//fYbM2bMYNy4cezZs6dc509JSWHs2LFMnDiR06dPs3v3bkaNGoVSTPvaJUuW0KNHD55++mlSUlJISUkhKCiIMWPGkJ+fzzfffGPY9sqVK2zZsoWJEycW+1oOHz6Mu7s7kyZNwtPTk//+978A5OTkcPr0abOW1H/66Sc6d+5cZH1ERARvv/02//rXvzh16hSffvopvr6+hudN/f8uiaIohISE0K9fP5566qlit+nWrRu//vorubm5pr1Y4fisfVUhhLkcOHBAAZRNmzaVuE12drbi4uKibNiwwbDu9u3bSkBAgLJgwQLDOlNK6gVLUHl5eUqNGjWUp556yrAuJSVFAZR9+/YVOs/dJa+uXbsqr776aomvoaSSuqIoSnBwsHLvvfcqiqKW+qpXr6788ssvhbYJCwtTxo4dW67zHz58WAGUixcvFnu+u9+nux/rPffcc8rQoUMNjxcuXKg0adJE0el0xR73pZdeUrp162bY9/HHH1cURVH279+vAEpiYmKx+5li+PDhysSJEwuty8rKUtzc3ArVehRk6v93aX766SdFo9Eo7du3NyzHjx8vtM2xY8dK/f8QwtmaFxRCmJNSjt6Z58+f586dO/Ts2dOwzsXFhW7dunH69OkKnb9du3aGv52cnKhTpw5t27Y1rNOX8q5cuVLifqC2gr57m/JSFAWNRgPAuXPnuHnzJgMHDiy0ze3bt+nYsWO5zt++fXv69+9P27ZtGTx4MIMGDeLRRx+ldu3aRsX19NNP07VrVy5fvkz9+vVZt26docFfcY4cOWIojY8aNYpRo0aRm5vLkSNHqFu3bqEGZIsWLSI9PZ0333wTgLy8POrUqcNff/2Fs3PZP3G3bt3C3d290LrTp0+Tm5tL//79S9zP1P/vkvTq1avMho7VqlUDKHfpX1Q9ktSFw7jnnnvQaDScOXPG7MfWarVFLhru3LlT6LGLi0uhx/rW6QUfA0V+uIvbz9RW7KdPn6Zx48YAZGdnA7Blyxbq169faDs3N7dynd/JyYnt27fzyy+/8MMPP7Bs2TJef/11Dhw4YDhPeXTs2JH27duzfv16Bg0axG+//caWLVtK3P7IkSOMHTsWUFvUu7i4sG3btmIbyZ08eZIHH3yw0HvQqFGjciV0AB8fH65du1ZonT55lsbU/++K0N8iqlu3rtmOKRyL3FMXDsPb25vBgwezYsUKbty4UeT5jIwMmjZtiqurKz///LNh/Z07dzh48CCtWrUq8dh169YlJSXF8DgrK4uEhATzvoAK+vHHHzlx4gSjR48GoFWrVri5uZGYmEizZs0KLQVLumXRaDT07NmT2bNnc/ToUVxdXfnqq6+K3dbV1ZX8/Pxin5s0aRLr1q0jOjqaAQMGlBjDhQsXyMjIMCRvZ2dnhg0bxpdffsmRI0eK3P8+efIk7du3NzyOi4srVGIuS8eOHTl16lShdffccw/VqlWzuW6QJ0+eJDAwEB8fH2uHImyUlNSFQ1mxYgU9e/akW7duzJkzh3bt2pGXl8f27dtZtWoVp0+f5rnnnuPll1/G29ubBg0asGDBAm7evFlql6Z+/fqxbt06HnnkEWrVqsWsWbNwcnKqxFdWWG5uLqmpqeTn55OWlsb3339PVFQUDz/8MOPHjwfAw8ODl156iRkzZqDT6ejVqxeZmZn8/PPPeHp6MmHChDLPc+DAAXbu3MmgQYOoV68eBw4c4M8//+Tee+8tdvtGjRpx4MABLl68SM2aNfH29karVcsOTzzxBC+99BIffvgh69evL/Gchw8fxtXVlTZt2hjWjR492tDo7PXXXzesVxSFM2fOFLrFcO3aNSIiIsp8bXqDBw8mIiKCa9euGW4ruLu78+qrr/LKK6/g6upKz549+fPPP/ntt98qpetbSX766ScGDRpktfML2ydJXTiUJk2acOTIEd566y1efPFFUlJSqFu3Lp07d2bVqlUAvP322+h0Op566imuX79Oly5d2LZtW6n3iSMiIkhISODhhx/Gy8uLuXPnWrWk/v333+Pv74+zszO1a9emffv2LF26lAkTJhiSKMDcuXOpW7cuUVFRXLhwgVq1atGpUydee+21cp3H09OTvXv38t5775GVlUXDhg1ZuHAhQ4cOLXb7l156iQkTJtCqVStu3bpFQkKCYUAXLy8vRo8ezZYtWxgxYkSJ5zxy5Aht2rTB1dXVsG7gwIHk5+dz+/btQtXvCQkJNG/enEOHDhnWDR06lLZt2xIbG8v27du5du0akZGRhS4SCmrbti2dOnXi888/55lnnjGs/9e//oWzszOzZs0iOTkZf39/nn322fK8bRaRk5PD5s2b+f77760Wg7B9Mva7EKLS9O/fn9atW7N06VKzHO+bb74hNjaWjz/+2LCuYcOG7N27l4YNGwJqdfwPP/zAK6+8UuJxtmzZwssvv8zJkycLXRTZklWrVvHVV1/xww8/WDsUYcOkpC6EsLhr166xe/dudu/ezcqVK8123JMnTxZqC3Hjxg0yMzMNCV2n07Fy5UpmzZpV6nEeeughfv/9dy5fvmxUe4PK5OLiwrJly6wdhrBxktSFEBbXsWNHrl27xvz582nRooXZjnvy5Ekee+wxw+MzZ84YkrxOp2PGjBk899xzBAYGlnksW5+YZtKkSdYOQdgBqX4XQjikBQsWsHnzZtq1a8fAgQMNvQKEcGSS1IUQQggHYZstQoQQQghhNEnqQgghhIOQpC6EEEI4CEnqQgghhIOQpC6EEEI4CEnqQgghhIOQpC6EEEI4CEnqQgghhIOQpC6EEEI4CEnqQgghhIOQpC6EEEI4CEnqQgghhIP4f5zOIe5FzM8TAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Polarization fraction as a function of the column density\n",
    "PvN_plot_Step2 = HtoH.PvN(CatP, CatdP, CatN, Nscale=[1e+21,5e+22], Pscale=[0.0,30.0], showfit='true', weighted='false', errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "c49996f1",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfkAAAF2CAYAAABtUOHpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABSpklEQVR4nO3deVxU9f4/8NeMCAgyo+Aeq2ZqgoBGRmmadl1um0vfvOaCW8sNReW6gJletRwM09RMv5VbXcqya9Ztsa6k2KLmMqiouSTImKgpCoKKwJzfH/1mvgzrzJkzc87MvJ6Pxzx0zhzOvM8MfN7ns5zPRyUIggAiIiJyO2q5AyAiIiLHYJInIiJyU0zyREREbopJnoiIyE0xyRMREbkpJnkiIiI3xSRPRETkppjkiYiI3BSTPBERkZtikneQsrIy/POf/0RZWZncoTiUJ5wnz9E98BzdA8/RNiolTGu7Zs0arFmzBnl5eQCArl27Yt68eRg8eDAA4Pbt2/jHP/6BzZs3o6ysDAMHDsTbb7+N1q1byxh1/YqLi6HValFUVASNRiN3OA7jCefJc3QPPEf3wHO0jSJq8sHBwUhLS8PBgwdx4MAB9OvXD0899RSOHTsGAJg+fTr+85//YMuWLcjKysKFCxcwbNgwmaMmIiJSNi+5AwCAJ554wuL5a6+9hjVr1mDv3r0IDg7GunXr8OGHH6Jfv34AgA0bNqBLly7Yu3cvHnjgATlCJiIiUjxFJPmqKisrsWXLFpSWliI+Ph4HDx5EeXk5Hn30UfM+nTt3RmhoKPbs2VNnki8rK7PozzAajSgsLERQUBBUKpXDz6O4uNjiX3flCefJc3QPPEf34AnnWFRUBAC4fv26xXYfHx/4+PjYdjBBIY4cOSL4+/sLjRo1ErRarfDVV18JgiAIGRkZgre3d4394+LihFmzZtV5vPnz5wsA+OCDDz744MMtHvPnz7c5tyqmJt+pUydkZ2ejqKgIn376KRISEpCVlSX6eKmpqUhOTjY/LyoqQmhoKAwGg9sO1lCSiooK/P7777jrrrvg5aWYXzMiIsUrLi5GSEgI8vPzodVqzdttrsVDQc313t7euPvuuwEAPXr0wP79+7FixQqMGDECd+7cwfXr19GsWTPz/pcuXUKbNm3qPF5dzRoajYZJ3gkqKipQXFwMjUbDJE9EJIJWq3WP0fW1MRqNKCsrQ48ePdC4cWNkZmaaXzt58iTy8/MRHx8vY4RERETKpogqVmpqKgYPHozQ0FDcuHEDH374IXbt2oVvv/0WWq0WEydORHJyMgIDA6HRaDBlyhTEx8dzZD0REVE9FJHkL1++jLFjx6KgoABarRbdunXDt99+i7/85S8AgOXLl0OtVmP48OEWk+EQERFR3RQx450zeMIsSUpSUVGB8+fPIzg4mH3yREQ2cLsZ74iIiEh6TPJERERuikmeiIjITTHJExERuSkmeSIiIjfFJE9EROSmmOSJiIjcFJM8kQcwGo24fv06jEaj3KEQkRMxyRN5ACZ5Is/EJE9EROSmmORJNqxdEhE5FpM8yYZJnojIsZjkiYiI3BSTPBERkZtikieH0Ov1GDlyJPR6vdyhEBF5LCZ5coiZM2ciNzcXM2fOlDsUIiKPxSRPDpGeno6IiAikp6fLHQoRkcdikieHiI2NxUcffYTY2Fi5QyEiO2zcuBEqlQoqlcq8rW/fvlCpVBg3bpx8gZFVmOSJiDyUKVnX9ti2bRsAoGXLlujZsyd69uwpb7ASGjduHFQqFfr27WvV/nV9RnPnznVsoBLwkjsAIiKSl7e3d41Wt8DAQADAY489hscee0yOsBQnJiYGPj4+5uchISEyRmMd1uTJKTjxDZFytW3bFnv37rV4PPzwwwBqb66vTUJCAjp27IiAgAB4e3sjLCwMSUlJKC4uNu9jqkGHh4dj06ZNCAkJgUajwbRp03Dr1i1MmzYNWq0WoaGhWLt2rcXxL1y4gAkTJqBdu3bw9vZG+/btsWjRIlRUVJj3MbVMjB07FvPnz0fbtm3RvHlzjB49Gjdu3AAA83sDQFZWlvncdu3a1eDn9Nlnn1l8Ri+88IJVn6+cWJMnpzAl+aZNm0Kt5rUlkbv5/PPP4eXlhQ4dOuDGjRs4e/YsVq1ahYKCAmzZssVi34KCArz00kto27Ytbty4gRUrVuC7777D77//Dq1WC4PBgMTERPTp0wddunTB1atX8cADD8BgMCAgIABdunTB8ePHMW/ePOTm5mL9+vUWx9+8eTN8fX3RokULXLx4ERkZGQgLC8Nrr72G2NhYlJaW4sqVKwgICMC9994LANBoNA2e43333YebN28iIiICY8aMwfTp0y1q9krE0paIyMOdO3euRn+zrbKysnDlyhVkZ2fjt99+w8svvwwA2LZtG27fvm2x7507d/Ddd9/h1KlT5ibvM2fO4PDhwzh+/Dh8fX1hNBqRlZUFAHjrrbdgMBjQunVr/Pbbbzh8+DA+/fRTAH+2NJw5c8bi+L6+vjhx4gTOnDmDHj16AAAyMzMB/FkbN3U/dO/e3Vwr7969e73n17x5cwQHB8PHxwfHjx9Hamoqxo4da/Pn5GxM8kREHs7b29s8uE7sILsdO3YgMjISTZo0gUqlwmuvvQYAqKiowB9//GGxb/PmzfHQQw9BrVYjNDQUABAZGYnw8HA0bdoUrVq1AgBcunQJAPDLL7+Yn7dq1QoqlQpDhgwBAAiCgH379lkcv1+/frjrrrugVqvRuXNni2OJsXfvXly9ehXZ2dn4/fff0a9fPwDAJ598AoPBIPq4zsDmeiIiD2fqkxcrIyMDM2bMMB8rJCQEV65cwdmzZwEAlZWVFvtXbRr38vKqsc3UkiAIgsXPVW1er8rPz8/iebNmzWocv/qxbFH1osfPzw9Dhw7F999/DwAwGAyKHoDHmjwREdnFdIEQEBCA3Nxc7Nu3DwMGDJDs+HFxcQD+TNibN282N7H/97//xUsvvYShQ4fadDzTRUFpaWmD++7evRuffvqp+ULl9u3b+Pzzz82vh4WF2fTezsaaPBER2aVbt24AgBs3bqB9+/bw8fFBUVGRZMdPTEzEe++9h99//x2dOnVCly5dcOPGDRgMBpSXl9vcN25qwj9w4ACioqLg7++PnTt3okmTJjX2PXv2LMaPHw9/f3+0b98e58+fx7Vr1wAA48ePx1133WX/CToQa/JERGSXiRMnIjk5GS1atMCNGzfQt29fLFy4ULLjt2zZEnv37sX48eMRFBSEY8eO4datW+jduzeWL19u8/EmTJiA4cOHQ6vVIicnB/v27avRpWDSq1cvvPjiiwgNDUVubi6MRiN69OiBtWvX4p133rH31BxOJdjTUeFCiouLodVqUVRUZNWtEmSfiooKnD9/HsHBwfDy8qrxvLZ9HM1oNKK4uBgajcbjbuNz9mdNROJJma88q6Qjybji5DauGDMRkT2Y5EkUV0yYRqMRJSUlLhUzEZE9mOTJYzDJE5GnYZInIiJyU0zyREREbopJnmTlin37RESugkmeZCUIAoqKipjkiYgcgEmeiIjITTHJE1EN7EYhcg9M8uRwTBiuh98ZkXtQRJLX6XSIi4tDQEAAWrVqhSFDhuDkyZMW+/Tt2xcqlcri8eKLL8oUMdnCaDSiqKjIpqUemWSIiOyniCSflZWFxMRE89KB5eXlGDBgQI1lAJ977jkUFBSYH6+//rpMEZOjWZPkeSFARFQ/RaxUsX37dovnGzduRKtWrXDw4EE8/PDD5u1+fn5o06aNs8MjhTIl+aZNm3rcgjNERNZQZMloWoc4MDDQYntGRgZatGiByMhIpKam4ubNm3Ueo6ysDMXFxRYPIiIiV1E9h5WVldl8DMUleaPRiGnTpuGhhx5CZGSkefuzzz6Lf/3rX9i5cydSU1PxwQcfYPTo0XUeR6fTQavVmh8hISHOCJ+IiEgSISEhFnlMp9PZfAxFNNdXlZiYiJycHPz4448W259//nnz/6OiotC2bVv0798fv/32Gzp06FDjOKmpqUhOTjY/Ly4uZqInIiKXYTAYLNaT9/HxsfkYikrykydPxpdffondu3cjODi43n179uwJADhz5kytSd7Hx0fUB0JkL6PRiOLiYmg0Go4VICLRNBqNRZIXQxElkCAImDx5Mj777DN8//33iIiIaPBnsrOzAQBt27Z1cHREtuGofyJSCkXU5BMTE/Hhhx/i888/R0BAAC5evAgA0Gq1aNKkCX777Td8+OGH+Otf/4qgoCAcOXIE06dPx8MPP4xu3brJHD0REZEyKSLJr1mzBsCfE95UtWHDBowbNw7e3t7YsWMH3nzzTZSWliIkJATDhw/H3LlzZYiWiIjINSgiyTc0E1pISAiysrKcFA05i16vx5QpU5CamoqwsDC5wyEicjuK6JMnzzRz5kzk5uaKui2EiIgaxiRPsklPT0dERARSU1PlDoWIyC0pormePFNsbCwyMjLkDoOIyG2xJk9EROSm7KrJl5eX4+LFi7h58yZatmxZY655IiIiko/NNfkbN25gzZo16NOnDzQaDcLDw9GlSxe0bNkSYWFheO6557B//35HxEpEREQ2sCnJL1u2DOHh4diwYQMeffRRbNu2DdnZ2Th16hT27NmD+fPno6KiAgMGDMCgQYNw+vRpR8VNREREDbCpuX7//v3YvXs3unbtWuvr999/PyZMmIC1a9diw4YN+OGHH9CxY0dJAiUiIiLb2JTkP/roI6v28/HxwYsvvigqICIiIpKGJLfQXblyBfv27UNlZSXi4uK4aAwpltFoRFFREby9vblCHBG5PbuT/L///W9MnDgR99xzD8rLy3Hy5EmsXr0a48ePlyI+IkkJgoCioiIEBQUxyZPL4PLFJJbNvy0lJSUWzxcsWIBffvkFv/zyC/R6PbZs2YKXX35ZsgCJiDwdly8msWxO8j169MDnn39ufu7l5YXLly+bn1+6dAne3t7SREdERESi2dxc/+233yIxMREbN27E6tWrsWLFCowYMQKVlZWoqKiAWq3Gxo0bHRAqERER2cLmJB8eHo6vvvoKH330Efr06YOkpCScOXMGZ86cQWVlJTp37gxfX19HxEouRK1Wo1mzZlCr1WxiJCKSiegRHCNHjsT+/ftx+PBh9O3bF0ajETExMUzwBMAyyRMRkTxEja7/+uuvceLECURHR+O9995DVlYWRo0ahcGDB2PhwoVo0qSJ1HESERGRjWyuZv3jH//A+PHjsX//frzwwgtYtGgR+vTpg0OHDsHX1xexsbH45ptvHBErkUNxBDMRuRubk/zGjRvx9ddfY/Pmzdi/fz8++OADAIC3tzcWLVqErVu3YvHixZIHSq4tJycHo0aNgl6vlzuUOjHJE5G7sTnJ+/v7Izc3FwBgMBhq9MHfe++9+OGHH6SJjhRLr9dj5MiRViVto9GIV199Fbm5uZg5c6YToiMiIkBEktfpdBg7dizatWuHPn36YNGiRY6IixRu5syZVidto9GIl156CeHh4UhPT3dCdEREBIhI8qNGjYLBYMDnn3+OvLw8PPXUU46IixQuPT0dERERFkm7vubue+65Bx988AFiY2OdGSb9f7a0vBCR+7Apyefn5wMAgoKCEBcXh2bNmtW7/++//y46MFK22NhYfPTRRxZJ29F92rUdn/3o1rGl5YWI3IdNST4uLg4vvPAC9u/fX+c+RUVFePfddxEZGYl///vfdgdIZFJbQjctOMMkX7/aWl6IyP3ZdJ/88ePH8dprr+Evf/kLfH190aNHD7Rr1w6+vr64du0ajh8/jmPHjqF79+54/fXX8de//tVRcZMTcQUs12dqeQkODpY7FCJyIptK7KCgICxbtgwFBQV466230LFjR1y5cgWnT58G8Gd//cGDB7Fnzx4meDdib5O4Xq/HlClTcOzYMYkjI3Id7FoiOYia8a5JkyZ4+umn8fTTT0sdD7mhlJQUGAwGpKWl4YknnpA7HCJZmJJ806ZN2SJGTsPfNHK4tLQ0hISEICUlpd79WNNxDn7ORJ6DSZ4cLjY2FqtWrULXrl3r3Y/Jxzn4ORN5DiZ5IiIiN8UkTy6HNVEiIuuIGnh37tw5HDlyBK1bt8b9998vdUxE9TIl+errJhARkSWbk/xHH32EcePGoby8HCqVyry0bMuWLR0RHxEREYlkc3P9ggUL8Oyzz+LXX3/Fd999BwANjpomsgeb54mIxLE5yZ89exbz58/HPffcg/79++Nf//oXNm/e7IjYiABIm+RPnTqFMWPGcKEWIvIINif5iooK+Pn5mZ937twZRqMRFy9elDQwIkdYtWoV8vLyHLpQC1d8IyKlEDW6ftOmTfj5559RUlICAPDy8sLNmzclDYzIEaZMmeLwde254hsRKYXNSb5379549dVX0atXLzRr1gwdO3bE7du3sW7dOuzcuRM3btxwRJxEknDGuvZc8Y2IlMLm0fVZWVkAgNOnT+PgwYM4dOgQDh06hDVr1kCn00GtVqNjx444ceKE5MESuQKu+EZESiHqPnkA6NixIzp27Ii//e1v5m25ubk4cOAA+yLJbnq9HklJSVi5cqVDa91ERO5M0hnvIiIi8D//8z9YvHixTT+n0+kQFxeHgIAAtGrVCkOGDMHJkyct9rl9+zYSExMRFBSEpk2bYvjw4bh06ZKU4ZOTqdVqaLVaqFSqGq85ol+bS94SkadRxLS2WVlZSExMxN69e/Hf//4X5eXlGDBgAEpLS837TJ8+Hf/5z3+wZcsWZGVl4cKFCxg2bJiMUZO91Go1cnNzMWrUqBqtP47o16665C0RkSdQRJLfvn07xo0bh65duyI6OhobN25Efn4+Dh48CAAoKirCunXrsGzZMvTr1w89evTAhg0b8PPPP2Pv3r21HrOsrAzFxcUWD1Keumrspn7t+prq9Xo9Ro0ahZycHKvea9y4cbh9+zaeeeYZu2ImInKG6jmsrKzM5mMoIslXV1RUBAAIDAwEABw8eBDl5eV49NFHzft07twZoaGh2LNnT63H0Ol00Gq15kdISIjjAyeb2VNjT0lJQW5uLnQ6nVX7b9y4Eb6+vvjkk09sfi8iImcLCQmxyGPWlnVVKS7JG41GTJs2DQ899BAiIyMBABcvXoS3tzeaNWtmsW/r1q3rnIQnNTUVRUVF5ofBYHB06CSCNTX2uqSlpSEiIgKpqalW7x8SEsJpmInIJRgMBos8Zm1ZV5Xo0fUAkJmZiczMTFy+fLnGlKPr168XdczExETk5OTgxx9/tCc0+Pj4wMfHx65jkLLFxsYiIyPDpv1XrVrFxZSIyCVoNBpoNBq7jiG6Jr9gwQIMGDAAmZmZuHLlCq5du2bxEGPy5Mn48ssvsXPnTot7jNu0aYM7d+7g+vXrFvtfunQJbdq0EXsKRBaMRiNKSkq4EA44NS+RuxBdk1+7di02btyIMWPG2B2EIAiYMmUKPvvsM+zatQsREREWr/fo0QONGzdGZmYmhg8fDgA4efIk8vPzER8fb/f7EwFM8lVVHRC5a9cuucMhIpFEJ/k7d+7gwQcflCSIxMREfPjhh/j8888REBBg7mfXarVo0qQJtFotJk6ciOTkZAQGBkKj0WDKlCmIj4/HAw88IEkMRPR/0tPTkZSUxKl5iVyc6Ob6SZMm4cMPP5QkiDVr1qCoqAh9+/ZF27ZtzY+PP/7YvM/y5cvx+OOPY/jw4Xj44YfRpk0bbN26VZL3JyJL9gyIJCLlEF2Tv337Nt555x3s2LED3bp1Q+PGjS1eX7ZsmdXHEgShwX18fX2xevVqrF692uZYiYiIPJHoJH/kyBHExMQAQI3JSGqbppRI6fR6PSZPnoy33nqLYz2IyC2ITvI7d+6UMg4i2aWkpCA/Px8pKSnm1RaJiFyZ4ibDIeWx93Yqo9GIW7duKX7UelpaGkJDQzm3PRG5Dbsmw7l+/TrWrVtnXjv+3nvvxcSJE6HVaiUJjpTB3tupTEnemrEXcoqNjcXq1avN3VBERK5OdE3+wIED6NChA5YvX47CwkIUFhZi+fLl6NChAw4dOiRljCQzR6wIR0REjie6Jj99+nQ8+eSTePfdd+Hl9edhKioqMGnSJEybNg27d++WLEiSl+l2qqqzEBIRkfKJTvIHDhywSPAA4OXlhVmzZuG+++6TJDgiIiIST3RzvUajQX5+fo3tBoMBAQEBdgVF5ExGoxHXr19X/MBA8lxcS4DEEp3kR4wYgYkTJ+Ljjz+GwWCAwWDA5s2bMWnSJIwcOVLKGIkcSs4kzwsMskbVwa9EthCd5JcuXYphw4Zh7NixCA8PR3h4OMaNG4enn34aS5YskTJGchGsbdiOSZ6swcGvJJboJO/t7Y0VK1bg2rVryM7ORnZ2tnmEPddx90x11Tays7MxZ84cHD9+XKbIiFwb1xIgseyeDMfPzw9RUVGIioqCn5+fFDGRC1Cr1WjWrBnU6v/7FaqrtjFnzhxcuHBBES08KpUKWq3WIm4iIndl0+j65ORkLFq0CP7+/khOTq53X1sWqCHXY0ryVdV1q93ixYvx97//HbNnz3ZihLVTq9VM8kTkMWxK8nq9HuXl5eb/14UL1FBV3bt3x/Lly3mfPRGRk9mU5KsuSrNp0yYEBwfXqBEJggCDwSBNdOQW1Go1mjRpwtozEZGTiS51IyIicOXKlRrbCwsLERERYVdQREREZD/RSb6uxUZKSkrg6+srOiAiIiKShs3T2poG3KlUKsybN89iRH1lZSX27dvHVbyoBpVKBY1GwyZ7IiInsjnJmwbcCYKAo0ePwtvb2/yat7c3oqOjMWPGDOkiJLdQ2y13RETkWDYnedPgu/Hjx2PFihXQaDSSB0VUG71ej+TkZCxbtgxt27aVOxwiIsUTXa3asGEDEzw5VUpKCnJzc5GSkiJ3KERELkF0ktfpdFi/fn2N7evXr1fEzGbkftLS0hAREYG0tLQG9+Wc8EREdiT5//3f/0Xnzp1rbO/atSvWrl1rV1BEVZkWvgFg9fzdYpJ8dnY2EhMTucAOEbkN0Un+4sWLtfaLtmzZEgUFBXYFRVSVaeEbUzO9o1a7mzNnDvLz89kdQA7BVRpJDqKTfEhICH766aca23/66Se0a9fOrqDIPdham66rEDQtfGNqpq/aN5+Tk4NRo0ZJUnAuXrwYoaGhVnUHkGeRovuHa8KTHEQn+eeeew7Tpk3Dhg0bcO7cOZw7dw7r16/H9OnT8dxzz0kZI7koWwvG6oWg6eejo6Mtmumr9s3rdDrJCs6YmBisXr3apuU82ffvGaT4nrkmPMnB5lvoTGbOnImrV6/ipZdewp07dyAIApo0aYLZs2ezuZNESU9PR1JSkrkQFAQBRUVFaNq0qcV+ptXu2rRpg9TUVOh0OtkKTlPh37RpU/McAJwTgGpT1yqNRI4kuhRSqVRYsmQJ/vjjD+zduxdHjhxBYWEh5s2bx1XoyIK1A9pMhaAtNenIyEhkZGTY9DOOxiRPREohuiZvkp+fj6tXr+LOnTs4e/asefuTTz5p76HJTVQd0JaVlVXnfkajESUlJTAajUyQREQSEJ3kz549i6FDh+Lo0aNQqVTmBWtMtfjKykppIiSXt3jxYkydOrXBAW1M8s6h1+uRlJSElStXIi4uTu5wiMiBRJekU6dORUREBC5fvgw/Pz/k5ORg9+7duO+++7Br1y4JQyRXVHWgkpgBbeQ4HOVN5DlEJ/k9e/Zg4cKFaNGiBdRqNRo1aoRevXpBp9MhKSlJyhjJBRmNRhQVFdW5JDHJh6O8iTyH6CRfWVmJgIAAAECLFi1w4cIFAEBYWBhOnjwpTXREJDkxAxyJyDWJ7pOPjIzE4cOHERERgZ49e+L111+Ht7c33nnnHbRv317KGElmrjBa3NQ9wDXriYj+j+gkP3fuXJSWlgIAFi5ciMcffxy9e/dGUFAQPv74Y8kCJPmZkryS1Xa/OhGRpxOd5AcOHGj+/913341ff/0VhYWFaN68Oe+TJ5ekVqvd9iKhamsMZ+cj8hyiSrPy8nL0798fp0+fttgeGBjIBE8uS8okr7Tpbl2hy4WIpCfqL75x48Y4cuSI1LGQh1JaQpSCO54TEbke0Zf1o0ePxrp166SMhTyUaY56JkQiImmJ7pOvqKjA+vXrsWPHDvTo0QP+/v4Wry9btszqY+3evRvp6ek4ePAgCgoK8Nlnn2HIkCHm18eNG4dNmzZZ/MzAgQOxfft2seGTC6s6Mx4REdVNdJLPyclB9+7dAQCnTp2yeM3WfvnS0lJER0djwoQJGDZsWK37DBo0CBs2bDA/9/HxsTFiciVqtRparRbFxcU1XmOSJyKyjugkv3PnTsmCGDx4MAYPHlzvPj4+PmjTpo1k70nKZhooVlJSIncoREQuy+Y++bNnz8oyVemuXbvQqlUrdOrUCX//+99x9erVevcvKytDcXGxxYOIiMhVVM9hZWVlNh/D5iTfsWNH/PHHH+bnI0aMwKVLl2x+Y1sMGjQI77//PjIzM7FkyRJkZWVh8ODB9a50p9PpoNVqzY+QkBCHxkj2Mc11L0cTvF6vx8iRIxtc756IyJlCQkIs8phOp7P5GDYn+eq1+K+//to8852j/O1vf8OTTz6JqKgoDBkyBF9++SX2799f72p3qampKCoqMj8MBoNDYyT7CIKA4uJiWVqJTKuypaSkOP29iYjqYjAYLPJYamqqzcdwyZkx2rdvjxYtWuDMmTN17uPj4wONRmPxIPdiGpxn7wRMplXZGlrvnojImarnMDEDzm1O8iqVqkah6uxZ7s6fP4+rV6+ibdu2Tn1fkk9tzflSzeLGVdmIyF3ZPLpeEASMGzfOfEVx+/ZtvPjiizXuk9+6davVxywpKbGolefm5iI7OxuBgYEIDAzEggULMHz4cLRp0wa//fYbZs2ahbvvvtti/nxyb6YJcxo1aiR3KDYzGo0oLi7mCnlE5HQ2J/mEhASL56NHj7Y7iAMHDuCRRx4xP09OTja/15o1a3DkyBFs2rQJ169fR7t27TBgwAAsWrSI98qTS+AKeUQkF5uTfNUJaaTSt2/fegdcffvtt5K/J0mPK53VzZToAwMDLRK9Xq9HUlISVq5cibi4OBkjJCJ3xGoFSUZsH7lKpZJlVL0z1TU/v2lk/8yZM2WKjIjcGZM8yU4QhFoHb3rC/eumkf3p6elyh0JEbohJnhRL6vvXG1ovXo411+0d2c8lbYmoPkzypFhS37+uxCRvLyZ5IqqP65Rm5LK8vLzQrl07eHnZNs6zrlqu0WjErVu3bE5srpjEiYjsIXoVOgDIzMxEZmYmLl++XKPAXb9+vV2BkWvT6/WYMGECgD9/F6QcOW5vkici8hSiqzQLFizAgAEDkJmZiStXruDatWsWD/JsM2fOxMmTJ3Hy5MlaR46zmZmIyPFE1+TXrl2LjRs3YsyYMVLGQ24iPT3dXJOvbeS4Kcn7+vo6OzSHqTqznatj14ay8PsgsUT/xty5cwcPPviglLGQG4mNjcW2bduwbdu2BkeO5+TkYMyYMTh+/LjVx5dzadq6uFPrBJOKsvD7ILFE/8ZMmjQJH374oZSxkIfS6XTIy8vDkiVLrP4Z09K0SkqoRqMRJSUlioqJiDyb6Ob627dv45133sGOHTvQrVs3NG7c2OL1ZcuW2R0ceYbU1FQsXrwYs2bNkjsUu1RN8qxxEZESiE7yR44cQUxMDIA/m1urcvbSs+TaIiMj8cEHH6CyshKFhYVyhwOAK8cRkXsQneR37twpZRxEOH78OF577TWsXr0a8fHxdh3L3oVfrF05rur7REVF2RMyEZHkWEUhxViyZAkMBoMk09g6a+EXLjBDREpmV5K/fv063njjDUyaNAmTJk3CsmXLUFRUJFVs5GFmz56NkJCQGtPYqtVqnDt3DmPGjKnRNVQXZy38wgVmPIMnLJZE7kl0kj9w4AA6dOiA5cuXo7CwEIWFhVi+fDk6dOiAQ4cOSRkjeYh7770Xq1atqnHLnVqtxsKFC5GXlwedTmfVsexd+MVaznofkhdbbMhViU7y06dPx5NPPom8vDxs3boVW7duRW5uLh5//HFMmzZNwhCJ/q/GnJqaKncoZu50XzzVjy025KpED7w7cOAA3n33XYtFR7y8vDBr1izcd999kgRHnsU0J31FRQWuX78OPz8/82uxsbHIyMiQMbqaTBPykPsztdgEBwfLHQqRTUTX5DUaDfLz82tsNxgMCAgIsCso8kyCIFgk+bpqyIcPH8acOXMsZshjnykRUU2ik/yIESMwceJEfPzxxzAYDDAYDNi8eTMmTZqEkSNHShkjEYxGI0pLSxEQEIBXXnkFBQUFeP31182vp6SksM+UiKga0Ul+6dKlGDZsGMaOHYvw8HCEh4dj3LhxePrpp22anpQ8m16vx6hRo3Ds2LF69zMajbh58yY0Gg2WLFmC0NBQi1vt0tLS2GdKRFSN6D55b29vrFixAjqdDr/99hsAoEOHDhb9qEQNMdXAX3nlFfj5+WHkyJHo3bt3vT8TGxuLVatWoWXLlhbb2GdKRGTJ7slw/Pz8EBUVhaioKCZ4smDNgi2mGjgA5OXlYf369c4Kj4jI7dlUk09OTsaiRYvg7++P5OTkevflAjUkCEKDSd40av7IkSPQ6XRWj+ew5thERJ7OpiSv1+tRXl5u/n9duEAN2apr167YsGEDCgoKrFoQRolJXq/XY/LkyXjrrbfg5eVlntOeE+UQwDXhSR42Jfmqi9Js2rQJwcHBNX5hBUGAwWCQJjryKGq1Gk2aNHHZQjAlJQX5+flISUmBSqUyj/bfsWOH3KGRApiSPJEziS5NIyIicOXKlRrbCwsLzX2sRJ40K1xaWhpCQ0ORlpZm1wxpnvSZEZFjiR5dLwhCrdtLSkrg6+srOiByL4IgoLi4GC1atHDZGrq1YmNjsXr1asTExMDLy0v0aH9rl7klImqIzUneNOBOpVJh3rx5FiPqKysrsW/fPsTExEgWIJGJKflpNBq5Q1GMquvZx8XF2XUso9GI4uJiaDQaXlwQuQmbk7xpwJ0gCDh69Ci8vb3Nr3l7eyM6OhozZsyQLkLySLUlHFOrQNOmTc37efpgpqqro+3atcuuY7EFgcj92JzkTYPvxo8fjxUrVrBWRQ4hCAJu3LiBsLCwehOOKclXVFQ4MbraVR1db2+t2lrp6elISkriTH9EVCvRffIbNmwAABw/fhz5+fm4c+eOxetPPvmkfZGRS1Or1dBqtbh06ZLcoTjN7Nmzce7cOcyePRvff/+9ebter8eECRMAAOvXr5f0AoAz/RFRfUQn+dzcXAwZMgRHjx6FSqUyD8Qz3SNfWVkpTYTkkkxJ3ppmX5VKBY1Gg+vXrzs+sFpI1eS/ePFivPTSS1i8eLHF9pkzZ+LkyZPm/9vbrE5EZC3RpVpSUhIiIiJw+fJl+Pn54dixY9i9ezfuu+8+FmJkE7n71aV6/5iYGKSlpdUYeJqeno5OnTqhU6dObFYnIqcSXart2bMHCxcuNN8apVar0atXL+h0OiQlJUkZI7khW9d/V6vVaNq0qUvOphgbG4tt27Zh27Ztipj9Tu6LKiJyHtF/5ZWVlQgICAAAtGjRAhcuXAAAhIWFmZsmyfNYO5GLaVR41eVi6ztmcXEx/Pz8mJgkwCRP5DlE/5VHRkbi8OHDAICePXvi9ddfx08//YSFCxeiffv2kgVIrsXaJJ+eno7w8HC8/PLLDe5rNBpRVFRU5wRMRERUO9FJfu7cuebCeeHChcjNzUXv3r3x9ddfY+XKlZIFSK7JlJjrSuCm1efCw8PNyfvYsWOYMmUKsrOz6z12dnY2kpKScOrUKanDJiJyK6KT/MCBAzFs2DAAwN13341ff/0VV65cweXLl9GvXz/JAiTXZGpi9/f3t7ofPS0tDQaDAXPmzKl3vzlz5uD8+fNYvXq1FKESEbktSTvlAgMDRQ2M2r17N5544gm0a9cOKpUK27Zts3hdEATMmzcPbdu2RZMmTfDoo4/i9OnTEkVNjmIaLGdt329KSgpCQkJq3IJW3eLFixEYGIhbt25ZPXCPpMHFc4hci+gk369fPyxYsKDG9mvXrtlcky8tLUV0dHSdNbPXX38dK1euxNq1a7Fv3z74+/tj4MCBuH37tqjYyTlycnIwZswY5OTkWLV/165dsWrVqgbXPoiJiYGfnx+uXr1a68A9JiLH4WdL5FpEJ/ldu3bhrbfewpAhQ1BaWmrefufOHWRlZdl0rMGDB+PVV1/F0KFDa7wmCALefPNNzJ07F0899RS6deuG999/HxcuXKhR4/d0SiuAdTod8vLyoNPpJD/25MmTERwcjLS0tBqvudJAPY50JyJHsqtk2bFjBy5evIgHHngAeXl5EoVkKTc3FxcvXsSjjz5q3qbVatGzZ0/s2bOnzp8rKytDcXGxxcPdmZJ8RUWFIpJ9amoqwsPDkZqaKvmx77nnHqxYscJh953beh+/WEzyRFSX6jmsrKzM5mPYVbK0bdsWWVlZiIqKQlxcnENmurt48SIAoHXr1hbbW7dubX6tNjqdDlqt1vwICQmRPDalUkqNPjIyEh988AEiIyPr3EelUlk9/a0zVV3djYhIDiEhIRZ5TEyrqOiS1TTAzsfHBx9++CGmTp2KQYMG4e233xZ7SEmlpqaiqKjI/DAYDHKHRLVQak02PT0dERER9U5Dq9frMWrUKPOYA7VajSZNmlicS0MXXEq5ICMi5TEYDBZ5TEyrqOiStXp/59y5c5GRkYE33nhD7CFr1aZNGwCosZrZpUuXzK/VxsfHBxqNxuJB8qiayKxJatnZ2RbJE7B+EJ9USdO0ult93QGm2r7p6rquJF/f+IDa4nVWVwERKVv1HObj42PzMUQn+dzcXLRo0cJi2/Dhw7F3716sX79e7GFriIiIQJs2bZCZmWneVlxcjH379iE+Pl6y9yHHEQTBPDGONUl4zpw5FskTsH4QnzMH3Zlq+1KPOWBXARFJRXSSDwsLq7WJNTIyEgkJCTYdq6SkBNnZ2eaZznJzc5GdnY38/HyoVCpMmzYNr776Kr744gscPXoUY8eORbt27TBkyBCx4ZOCLV68uEbydOQgPrFMs/bVN+ZADGu6CoiIrGHTevLJyclYtGgR/P39kZycXO++y5Yts/q4Bw4cwCOPPGLxPgCQkJCAjRs3YtasWSgtLcXzzz+P69evo1evXti+fTt8fX1tCZ8UTqVSoWnTpujQoQP+93//F35+fubXIiMjkZKSAp1Oh6CgIBmjdDxTV0FwcLDcoRCRi7Mpyev1epSXl5v/XxdbZ73r27dvvc2rKpUKCxcuxMKFC206Lsnv2LFjSEtLw5w5cxASEoLi4uI6m+pNM+QBwKFDh7B27Vq89dZb5n7xtLQ05OXlYc6cOUhISMBbb72Ft99+m902RER1sCnJ79y5EwBQXl4OtVqNtWvXomPHjg4JjNyDKTHrdDoMHDgQRUVFVv3cqlWr8Pvvv2PmzJnYsWMHgD+nvV2yZAnS0tKQlJSE8+fPIyUlxebJl4iIPIWoPvnGjRvjyJEjUsdCbiglJcWqvnS1Wg2tVmtuBZoyZQrCw8Mt+qW7du2KjIwMdO/eHS+//DJCQkJqnfHOmZR6nz8REWDHwLvRo0dj3bp1UsZCTuSs+7O7du3a4IQ4QM375e+55x588MEHdd7CZprn3lEz3llLqff5ExEBNjbXV1VRUYH169djx44d6NGjB/z9/S1et2XgHTmfKcnbskoc1c/WVffqYjQaUVJSwglyiMhuopN8Tk4OunfvDgA4deqUxWtilpslkpIczeiHDx9GUlISVq5caVcLA5M8EUlFdJI3DcIjUoLqzeZyNKNXncTGNFjQlej1evNFSlxcnNzhKAq7ZchV8TeW3IISCuGGJrFxxDgIKc+bM+3VTQm/X0Ri2P0be/z4cWzfvh1ffPGFxYOUT+7FUXJycjBq1KgG52g3LQRz7NgxJ0UmTkPz3Ved3lcqUiYfzrRH5H5Elwxnz55FdHQ0IiMj8dhjj2HIkCEYMmQIhg4diqFDh0oZIzmII5KOLYur6HQ6q2qOphqmtbfLWXvxQJasWZSHiFyL6CQ/depURERE4PLly/Dz88OxY8ewe/du3HfffQ5ZV56sJ+cqZqaEnJiYaK59Hzt2rNakm5qaalXN0VTDTElJsSoGay8enKG+Vgi5W1KIyP2JTvJ79uzBwoUL0aJFC6jVaqjVavTq1Qs6nQ5JSUlSxkg2SklJkS3JLVmyBKGhoQD+XGho6dKlSEtLQ25ubo0kHRkZiYyMDHPNsa6mZ9NCMF27drUqBmsvHpzB9F3U1gphNBpx7do13Lhxg4meiBxCdJKvrKxEQEAAAKBFixa4cOECgD9Xpzt58qQ00ZEoaWlpsiW56OhorFy5EitXrkT79u2xdOlSzJkzBxEREUhLS6u3KV2q/uXqFw/OVL12bvou6mqFEARBEbfLWduqwLXuiVyL6NI0MjIShw8fBgD07NkTr7/+On766ScsXLgQ7du3lyxAsp6pAAYge99qTEyMOYaqSdfUlG5t03tVrnD/ePVxDra2QsjFaDSiqKio3oWiAI7AJ3I1opP83LlzzQXZwoULkZubi969e+Prr7/GypUrJQuQrGcqgMUkUGfQ6/UoLS1Fq1atRM05r5RaryfjCHwi1yI6yQ8cOBDDhg0DANx999349ddfceXKFVy+fBn9+vWTLECynqkAlnvRlrqkpKTg8uXL8Pf35wju/+/YsWOYMmWKyzR/cwQ+kWuRdGaHwMBATmkrI6UXwKb+6dTU1BqrzlWl1+sxZcoUxd8XX5/q4wtM0+xWP9+0tDQYDAbFtr4QkWuzaVrb5ORkq/flAjX2MRqNKC4uhkajcZtZtkz908D/JcGSkpIa+6WkpMBgMCAtLQ1PPPGEs8OUhOn8KioqzM+1Wi1u3bplsV9KSgoWLVpk0foi1UI37sQd/x6InMGmJG9tkyJr8/azd5U4088HBga6XKGYlpaGxMREt6ndqlSqOpNTbUvmMsnXxFUTicSxKclzURrXYRrlXdctaXq9HhMmTAAArF+/XlELksTGxmLVqlVo2bJlg/uaCn9TEjU19aempiIsLMwJ0TbMVKsnInI20avQAcD169exbt06nDhxAsCftZIJEyZAq9VKEhw5zsyZM83zGcycOdNlZymsXsMz3WGg0+nw2GOPyR0eEZGsRLd7HThwAB06dMDy5ctRWFiIwsJCLFu2DB06dMChQ4ekjJEcID09HZ06dUKnTp0UeTtUfU3c9THdYZCamuqgyByPE84QkVREJ/np06fjySefRF5eHrZu3YqtW7ciNzcXjz/+OKZNmyZhiOQIsbGx2LZtG7Zt2+a00fi2TGZT2+x3KpWqwT5Z0+C+yMhISWKWkrXJ2x0nnOE8/UTysKsmP3v2bHh5/V+Lv5eXF2bNmoUDBw5IEhy5F1OSb2hWtbpUHZBmTcJXGmuTtztOOMMkTyQP0SWkRqNBfn5+je0Gg8E8pz15Hr1ej8TERGRnZ0t+bNNqdocPH0bz5s0REBDgUkm+tuRd22Q4Sp/vgIhch+gScsSIEZg4cSI+/vhjGAwGGAwGbN68GZMmTTLPn06ew1RTmz17NvLz8zFnzhxJj69SqbB06VLk5uZi9uzZkixkYy9bF9SpLXlzMhwiciTRpeTSpUsxbNgwjB07FuHh4QgPD8e4cePw9NNPY8mSJVLGSC7AlOQXL16M0NBQvPrqq+aJbpo3b47mzZtbnQxru09crVZj6dKlaN++vVOasa1J4LXtU98qe7VJSUlBSEiIYqciJiLXJjrJe3t7Y8WKFbh27Rqys7ORnZ2NwsJCLF++HD4+PlLGSC4kJiYGq1evRrdu3cxJPjAwEIGBgQCAW7duNdgvW9dkMFVrwlItS1tfDGKOb1plz9pBc7VNhkNEJBWbS0ij0YglS5bgoYceQlxcHBYuXIi7774bUVFR8PPzc0SMHskdb6MyGo24deuW6IF3VTk6yYuVmppq0e9eX5x6vR5jxozBqVOnnB0mEXkIm0vI1157DXPmzEHTpk1x1113YcWKFUhMTHREbB7NHW+jAv7sW3f0gDnTYjByXABERkYiIyPDXDOvL8nPnDkTeXl5WLVqlbPDJCIPYXMp+P777+Ptt9/Gt99+i23btuE///kPMjIyeGuMxOy5jaq+Fd7kpFar4e/vb55P31G1cblq+bZ+7unp6QgPD8eUKVMcHFnD9Ho9Ro0ahZycHLlDISIJ2VwK5ufn469//av5+aOPPgqVSoULFy5IGpins+c2KqU2ZavVajRp0sQcl1LjFEvMaPsPPvgA99xzj4Mja1jV6YCJyH3YXLpWVFTA19fXYlvjxo1RXl4uWVDkHHI2a9vC3S4GlMgdpgMmoppsXqBGEASMGzfOYgT97du38eKLL8Lf39+8bevWrdJESA7jqORpNBpRVFQkWRdO9bXZSXqm6YCJyL3YnOQTEhJqbBs9erQkwZDrM/W7l5aWwt/f3+ICQq/XY/LkyXjrrbcQHx8vY5RUG7GLAhGRctmc5Dds2OCIOMhNmO5xr6yshI+Pj0XCSElJQX5+PlJSUpCVlSXJ++n1eiQlJWHlypWIi4uT5JieSsndIvyeicRR3l8ziVa1kFbiffZpaWkIDQ1tcHa3upJNbdtd/VZDV1toR64LAVf/nonk4holC1mlagHcUKHo6MLatJhM9YVXVq9e3eAdA7YkeVdfsa2u2f2USuzvjb0Xna7+PRPJxTVKFrJZQ4Wio5N8Wloazp49ixkzZjjk+CauvmKbq9zhYC97a+Ku/j0TycW9SxYXVXXtbbHrcDu7UKxeU5s1axaCg4O58EoVtV1YKbkfXEqsiRPJw71LFhdlNBpRXFwMjUYDAKKSfPXj2XuMhphqaikpKdBqtbUuvGJaa97aJltnxO1M1iZ0d0z8rIkTycMlSpF//vOfUKlUFo/OnTvLHZZDSVnQOyNZmmpqS5curTPuqqPrreFuSd5aSkzynvpdELk65ZQiDejatSsKCgrMjx9//FHukKgKa2pq1o6ud2VKTNBSYJInck023ycvFy8vL7Rp08bq/cvKylBWVmZ+Xlxc7IiwyAam0fUxMTEW26t2T7h6cjQleSIie1XPWz4+PhazzVrDZUrU06dPo127dmjfvj1GjRqF/Pz8evfX6XTQarXmR0hIiJMiJVvXSWctkYioppCQEIs8JmYBKZdI8j179sTGjRuxfft2rFmzBrm5uejduzdu3LhR58+kpqaiqKjI/DAYDE6M2LPJuU66uzaXE5HnMRgMFnlMzAJSLtFcP3jwYPP/u3Xrhp49eyIsLAyffPIJJk6cWOvPiGnWIGmkp6dj8uTJePHFF+06jtFoRElJiU01fDaXk1ju1G1E7kGj0ZjvshLLJZJ8dc2aNcM999yDM2fOyB0K1cK0Tvoff/xhsd3W2d3EJHkSx5M/a1Ny9/Pzw/Xr111qBkKihrjkb3JJSQl+++03tG3bVu5QqA7Hjh3DlClTLO6Jd7UpXD2Jpyd5jgkhd+USpe2MGTOQlZWFvLw8/Pzzzxg6dCgaNWqEkSNHyh2aQyhxcRlbpaWlwWAwWH1PvFI5MgFw/AAROZpLlC7nz5/HyJEj0alTJzzzzDMICgrC3r170bJlS7lDcwh3WHErNTUV4eHheP3110Ufw9YZ8hyBSd4xWHsmcg6XKF02b96MCxcuoKysDOfPn8fmzZvRoUMHucNymKrzfEtRq5ejZSAqKgqbN29Gjx49rP6Z6gW/rTPkkSVbEqmzu1LcNcm763mR63KJJO9pqs4eJ0Wt3lVaBqoXkJ4wQ54jKTnJuysmeVIa/kUrnBSrdyl9BbC6WhqsWX/ek5u8ncne1iB+T0Ty4F+cwkmxeld0dDTeffddREdHSxiZdKquYGcrJg9pNJTE7W0N4vdEJA/+xbmpqoWqvbdHWdMEaU0hXtc+ppYGd2qWd7Vm24aSuNytQbxIIBKHfzFuytlL1dqT5N1xrXFXS/INJXG5vyNnJHl3uHWVqDomeVIsJdxC5ykFv9xJXE6m7zgxMdElBqgS2YJJniSnVquh1WqhUqls+rnqCVXuW+iMRiOmT5+Os2fPumTB7ykXKPYydVUAUPQAVSIxmORJcmq1GoGBgQgMDLSpeTUlJcWiJiX3LXRGoxEzZ8502YJfybdOSn0BUlf3iDXdJqauitWrV3tsawa5LyZ5cggxfahpaWkWCdWaW+gcLTIyEhkZGS5Z8Ms9WK4+Ul+A2JPkpeyqYOsJKQ2TPDXIWQWX0vqF9Xo9Ro0ahZycHLlDEUVpn2dVSr4AsYeSW0/IMzHJK5DSbheSq+CSexY203nrdDrzNlcbNa9UUl+A1HUh6uzV9dz14oVclzKyCFkkD6mT/OHDh5GUlITDhw+L+nm5Ci65k7zpvFNTU83brE3ybLZ1rrouRA8ePIiEhAQcPHjQKXEoufWEPBOTvEI4soY4e/Zs5OfnY/bs2aJ+3tEFl9JaLkxiY2ORkZGByMhIm3/W1ZptlfodWKuuC1G579Agkptr/kV7ECkKX6U3Ibp6gqmN0j/z6lz9O6jrQlTuOzSI5OYldwBUP1Phaw9TARgcHCxNUB5CrVajefPm5v/bgp+5c9U3m+Lq1asRExMjT2BEMmOSJ8Vw9iCphpju93dVUtbOlV7Tl+JimMgdKfMvlhSDo8mVQcz34ElJvi7WDDo1nZuXl5dLniNRffjbrBBKHY1tNBpRXFwMjUbj8MKv+mh6JSYWub4nd7vYctZ3a82gUyZ5cmf8bVYIJY/GlivZKjHJy/U9KfUiUCxnfbfOHgCpxN9Z8mz8TVQIVxuN7ShKr7HK9T0p+SJQyZx93zqTPCkNfxMVQqmTaDizBqnX6/Hss8/ip59+UmySl+t74kVg7ZR+UUgkNyZ5qpcza5C1TSNLf+rRowe++eYb9OjRQ+5QFIVJnqh+TPISUHpBY08TojNrkOnp6WjVqhVKS0tdvu9Z6mZbNgPXzt3GKhBJjSWGBNw5yTuzeTo2Nhb+/v64fPmyy/c9Myk7R0MtTfweyNPxN58UJTU1lX3PZLWGWpqcneQPHDiA3r1748CBA055P6KGMMlLQIomQ9Y4/hQZGYmMjAzFDUA04fekLEobsDpjxgycPn0aM2bMkDsUIgBM8pKor8nQ2qZ8pSYPpcYlF34eyqK074N3QZDSKOMvw8XV94et9P76hjizEFWr1dBqtVCpVA5/L3IPSkjyVf/GeRcEKQ0XqJEAVxyThmlBGC8vL8XUzIgaYkryTZs2NU+NS6QULEkdjLf42EYJNTMiW/BvnJSMJamDcTpSIvfGv3FSMiZ5CdRX++RAHCL3xr9xUjKVIAiC3EE4Q3FxMbRaLYqKiqDRaJz2vhUVFTh//jyCg4Ph5cUhEETuhn/jJDUp8xVr8kREduA4ElIyXnY6GAsAIvdm+hsnUiImeQdjAUBERHJh9ZKIiMhNMckTERG5KSZ5IiIiN+VSSX716tUIDw+Hr68vevbsiV9++UXukIiIiBTLZZL8xx9/jOTkZMyfPx+HDh1CdHQ0Bg4ciMuXL8sdGhERkSK5TJJftmwZnnvuOYwfPx733nsv1q5dCz8/P6xfv17u0IiIiBTJJW6hu3PnDg4ePIjU1FTzNrVajUcffRR79uyp9WfKyspQVlZmfl5UVATgz5mEiIiIlMqUp0x5y8THxwc+Pj42HcslkvyVK1dQWVmJ1q1bW2xv3bo1fv3111p/RqfTYcGCBTW2h4SEOCRGIiIiKYWGhlo8nz9/Pv75z3/adAyXSPJipKamIjk52fzcaDSisLAQQUFBUKlUDn//4uJihISEwGAwOHWufGfzhPPkOboHnqN78IRzLCoqQmhoKM6dO2cxmZqttXjARZJ8ixYt0KhRI1y6dMli+6VLl9CmTZtaf6a2Zg05Zp7TaDRu+4tYlSecJ8/RPfAc3YMnnGOzZs08Y4Eab29v9OjRA5mZmeZtRqMRmZmZiI+PlzEyIiIi5XKJmjwAJCcnIyEhAffddx/uv/9+vPnmmygtLcX48ePlDo2IiEiRXCbJjxgxAn/88QfmzZuHixcvIiYmBtu3b68xGE8pfHx8MH/+fFF9KK7EE86T5+geeI7ugedoG5UgCIIEMREREZHCuESfPBEREdmOSZ6IiMhNMckTERG5KSZ5IiIiN8Ukb6fdu3fjiSeeQLt27aBSqbBt27Ya+5w4cQJPPvkktFot/P39ERcXh/z8fOcHK1JD51hSUoLJkycjODgYTZo0MS8g5Ep0Oh3i4uIQEBCAVq1aYciQITh58qTFPrdv30ZiYiKCgoLQtGlTDB8+vMYETUrW0DkWFhZiypQp6NSpE5o0aYLQ0FAkJSXVmD9byaz5Hk0EQcDgwYPr/LtVKmvPcc+ePejXrx/8/f2h0Wjw8MMP49atWzJEbDtrzvHixYsYM2YM2rRpA39/f3Tv3h3//ve/ZYrYdmvWrEG3bt3Mk/rEx8fjm2++Mb8uVXnDJG+n0tJSREdHY/Xq1bW+/ttvv6FXr17o3Lkzdu3ahSNHjuCVV16Br6+vkyMVr6FzTE5Oxvbt2/Gvf/0LJ06cwLRp0zB58mR88cUXTo5UvKysLCQmJmLv3r3473//i/LycgwYMAClpaXmfaZPn47//Oc/2LJlC7KysnDhwgUMGzZMxqht09A5XrhwARcuXMDSpUuRk5ODjRs3Yvv27Zg4caLMkVvPmu/R5M0333TKFNdSs+Yc9+zZg0GDBmHAgAH45ZdfsH//fkyePBlqtWsU+dac49ixY3Hy5El88cUXOHr0KIYNG4ZnnnkGer1exsitFxwcjLS0NBw8eBAHDhxAv3798NRTT+HYsWMAJCxvBJIMAOGzzz6z2DZixAhh9OjR8gTkALWdY9euXYWFCxdabOvevbvw8ssvOzEyaV2+fFkAIGRlZQmCIAjXr18XGjduLGzZssW8z4kTJwQAwp49e+QK0y7Vz7E2n3zyieDt7S2Ul5c7MTLp1HWOer1euOuuu4SCgoJaf6ddSW3n2LNnT2Hu3LkyRiWt2s7R399feP/99y32CwwMFN59911nhyeZ5s2bC++9956k5Y1rXNa5KKPRiK+++gr33HMPBg4ciFatWqFnz54u1TRojQcffBBffPEFfv/9dwiCgJ07d+LUqVMYMGCA3KGJZmqiDgwMBAAcPHgQ5eXlePTRR837dO7cGaGhoXUud6x01c+xrn00Gg28vFxm3iwLtZ3jzZs38eyzz2L16tV1rn3hSqqf4+XLl7Fv3z60atUKDz74IFq3bo0+ffrgxx9/lDNMu9T2PT744IP4+OOPUVhYCKPRiM2bN+P27dvo27evTFGKV1lZic2bN6O0tBTx8fHSljcSX4h4NFSrEZhqCX5+fsKyZcsEvV4v6HQ6QaVSCbt27ZIvUDtUP0dBEITbt28LY8eOFQAIXl5egre3t7Bp0yZ5ApRAZWWl8NhjjwkPPfSQeVtGRobg7e1dY9+4uDhh1qxZzgxPErWdY3V//PGHEBoaKsyZM8eJkUmnrnN8/vnnhYkTJ5qf1/Y77SpqO8c9e/YIAITAwEBh/fr1wqFDh4Rp06YJ3t7ewqlTp2SMVpy6vsdr164JAwYMMJc7Go1G+Pbbb2WKUpwjR44I/v7+QqNGjQStVit89dVXgiBIW9645uW5izAajQCAp556CtOnTwcAxMTE4Oeff8batWvRp08fOcOTzKpVq7B371588cUXCAsLw+7du5GYmIh27dpZXIm6isTEROTk5Lh0zachDZ1jcXExHnvsMdx77702r1+tFLWd4xdffIHvv//eZfptG1LbOZrKnRdeeMG8tkdsbCwyMzOxfv166HQ6WWIVq67f1VdeeQXXr1/Hjh070KJFC2zbtg3PPPMMfvjhB0RFRckUrW06deqE7OxsFBUV4dNPP0VCQgKysrKkfRNJLkdIEISaNYKysjLBy8tLWLRokcV+s2bNEh588EEnRyeN6ud48+ZNoXHjxsKXX35psd/EiROFgQMHOjk6+yUmJgrBwcHC2bNnLbZnZmYKAIRr165ZbA8NDRWWLVvmxAjtV9c5mhQXFwvx8fFC//79hVu3bjk5OmnUdY5Tp04VVCqV0KhRI/MDgKBWq4U+ffrIE6xIdZ3j2bNnBQDCBx98YLH9mWeeEZ599llnhmi3us7xzJkzAgAhJyfHYnv//v2FF154wZkhSqp///7C888/L2l5wz55B/L29kZcXFyNWz9OnTqFsLAwmaKSVnl5OcrLy2uM2m3UqJG5RuEKBEHA5MmT8dlnn+H7779HRESExes9evRA48aNLZY7PnnyJPLz811mueOGzhH4swY/YMAAeHt744svvnCpu0CAhs8xJSUFR44cQXZ2tvkBAMuXL8eGDRtkiNh2DZ1jeHg42rVr59LlTkPnePPmTQBw+XKnOqPRiLKyMmnLG2muPzzXjRs3BL1eL+j1egGAue/93LlzgiAIwtatW4XGjRsL77zzjnD69Glh1apVQqNGjYQffvhB5sit19A59unTR+jatauwc+dO4ezZs8KGDRsEX19f4e2335Y5cuv9/e9/F7RarbBr1y6hoKDA/Lh586Z5nxdffFEIDQ0Vvv/+e+HAgQNCfHy8EB8fL2PUtmnoHIuKioSePXsKUVFRwpkzZyz2qaiokDl661jzPVYHF+uTt+Ycly9fLmg0GmHLli3C6dOnhblz5wq+vr7CmTNnZIzceg2d4507d4S7775b6N27t7Bv3z7hzJkzwtKlSwWVSmXu11a6lJQUISsrS8jNzRWOHDkipKSkCCqVSvjuu+8EQZCuvGGSt9POnTsFADUeCQkJ5n3WrVsn3H333YKvr68QHR0tbNu2Tb6ARWjoHAsKCoRx48YJ7dq1E3x9fYVOnToJb7zxhmA0GuUN3Aa1nR8AYcOGDeZ9bt26Jbz00ktC8+bNBT8/P2Ho0KFCQUGBfEHbqKFzrOt7BiDk5ubKGru1rPkea/sZV0ry1p6jTqcTgoODBT8/PyE+Pt6lKhbWnOOpU6eEYcOGCa1atRL8/PyEbt261bilTskmTJgghIWFCd7e3kLLli2F/v37mxO8IEhX3nCpWSIiIjfFPnkiIiI3xSRPRETkppjkiYiI3BSTPBERkZtikiciInJTTPJERERuikmeiIjITTHJExERuSkmeSIiIjfFJE9EROSmmOSJiCRy9epVtGrVCnl5eTb93N/+9je88cYbjgmKPBqTPHmccePGQaVSQaVSoXHjxmjdujX+8pe/YP369ZIuU9m3b19Mmzat3n1McdT1+Oc//ylZPEpjzefjSN9++22Dn/93331n0zFfe+01PPXUUwgPDzdvGzduHIYMGWKx36effgpfX19zYp87dy5ee+01FBUV2XtaRBaY5MkjDRo0CAUFBcjLy8M333yDRx55BFOnTsXjjz+OiooKp8VRUFBgfrz55pvQaDQW22bMmOG0WKRy584dl3i/hx9+2OKzDgoKwiuvvGKxrX///lYf7+bNm1i3bh0mTpxY737vvfceRo0ahTVr1uAf//gHACAyMhIdOnTAv/71L1HnQlQn6RbOI3INCQkJwlNPPVVje2ZmpgBAePfddwVBEISwsDBh+fLlFvtER0cL8+fPNz/fsmWLEBkZKfj6+gqBgYFC//79hZKSEiEhIcHm5Vo3bNggaLXaWl+rrKwUFi9eLISHhwu+vr5Ct27dhC1btphf79OnjzB58mRh6tSpQrNmzYRWrVoJ77zzjlBSUiKMGzdOaNq0qdChQwfh66+/tjhunz59hMTERCExMVHQaDRCUFCQMHfuXPMywQ29b9VjTJ06VQgKChL69u0rCIIgfPPNN8JDDz0kaLVaITAwUHjsscfM65nX9/k09LnX9X7WxFqX8+fPCwAslvq01ZYtW4SWLVvW2F71923JkiWCr6+vsHXr1hr7LViwQOjVq5fo9yeqDWvyRP9fv379EB0dja1bt1q1f0FBAUaOHIkJEybgxIkT2LVrF4YNGwZBELBixQrEx8fjueeeM9cKQ0JCRMem0+nw/vvvY+3atTh27BimT5+O0aNHIysry7zPpk2b0KJFC/zyyy+YMmUK/v73v+N//ud/8OCDD+LQoUMYMGAAxowZg5s3b1oce9OmTfDy8sIvv/yCFStWYNmyZXjvvfesfl/TMby9vfHTTz9h7dq1AIDS0lIkJyfjwIEDyMzMhFqtxtChQ2E0Gu3+fGp7P2tjrY1erwcAdO/e3eoYqvvhhx/Qo0ePOl+fPXs2Fi1ahC+//BJDhw6t8fr999+PX375BWVlZaJjIKpB7qsMImerqyYvCIIwYsQIoUuXLoIgNFyjPHjwoABAyMvLq/VYffr0EaZOnWp1XHXV5G/fvi34+fkJP//8s8X2iRMnCiNHjjS/V9VaYEVFheDv7y+MGTPGvK2goEAAIOzZs8cixi5duphr7oIgCLNnzxa6dOli1fuajhEbG9vg+f3xxx8CAOHo0aPmn6vt87GmJl/9/ayNtS4LFiwQQkJCGtyvf//+wqlTp2p97amnnhImTJhQY3tCQoLg7e0tABAyMzPrPPbhw4fr/X0iEsNL1isMIoURBAEqlcqqfaOjo9G/f39ERUVh4MCBGDBgAJ5++mk0b95c0pjOnDmDmzdv4i9/+YvF9jt37iA2Ntb8vFu3bub/N2rUCEFBQYiKijJva926NQDg8uXLFsd54IEHLM45Pj4eb7zxhtXvC6DWGuzp06cxb9487Nu3D1euXDEPaszPz0dkZKRV516X6u9nS6y1OXTokFW1+NOnT6NDhw61vnbr1i34+vrW+lq3bt1w5coVzJ8/H/fffz+aNm1aY58mTZoAQI2WFiJ7MMkTVXHixAlEREQAANRqNQRBsHi9vLzc/P9GjRrhv//9L37++Wd89913WLVqFV5++WXs27fPfAwplJSUAAC++uor3HXXXRav+fj4mP/fuHFji9dMdw9UfQ7A6jsIrH1fAPD396/x80888QTCwsLw7rvvol27djAajYiMjGxwoFxDn3tt72dLrLU5dOgQJk2aVGP7sWPHMGHCBNy6dQsjRoxAmzZtoFbX3svZokULXLt2rdbX7rrrLnz66ad45JFHMGjQIHzzzTcICAiw2KewsBAA0LJlywbjJbIWkzzR//f999/j6NGjmD59OoA/C9uCggLz68XFxcjNzbX4GZVKhYceeggPPfQQ5s2bh7CwMHz22WdITk6Gt7c3Kisr7Y7r3nvvhY+PD/Lz89GnTx+7j1fdvn37LJ7v3bsXHTt2tOt9r169ipMnT+Ldd99F7969AQA//vijxT51fT7WfO7V2RPrlStXYDAYatTky8rKMGLECHz00UeIiorCkCFDLFpLqouNja13dHxYWBiysrLMiX779u0WiT4nJwfBwcFo0aKFTfET1YdJnjxSWVkZLl68iMrKSly6dAnbt2+HTqfD448/jrFjxwL4cyDexo0b8cQTT6BZs2aYN28eGjVqZD7Gvn37kJmZiQEDBqBVq1bYt28f/vjjD3Tp0gUAEB4ejn379iEvLw9NmzZFYGBgnbXA+gQEBGDGjBmYPn06jEYjevXqhaKiIvz000/QaDRISEiw67PIz89HcnIyXnjhBRw6dAirVq3CG2+8Ydf7Nm/eHEFBQXjnnXfQtm1b5OfnIyUlxWKfuj6fhj53qT+jQ4cOAag56G7btm3o06ePucujS5cuaNeuXZ3HGThwIFJTU3Ht2rU6u2xCQkKwa9cuPPLIIxg4cCC2b98OjUYD4M+BewMGDKj3PIlsxSRPHmn79u1o27YtvLy80Lx5c0RHR2PlypVISEgwJ+LU1FTk5ubi8ccfh1arxaJFiyxqlBqNBrt378abb76J4uJihIWF4Y033sDgwYMBADNmzEBCQgLuvfde3Lp1C7m5uRaTpNhi0aJFaNmyJXQ6Hc6ePYtmzZqhe/fumDNnjt2fxdixY3Hr1i3cf//9aNSoEaZOnYrnn3/ervdVq9XYvHkzkpKSEBkZiU6dOmHlypXo27eveZ+6Pp+GPve6iI1Vr9ejdevWNRL40aNHERMTY35+8OBBDBo0qM7jREVFoXv37vjkk0/wwgsv1LlfcHCwRaL/9ttv4e3tjW3btmH79u0NnieRLVRC9c4vIvIYffv2RUxMDN588025Q1GcZcuW4dy5c1ixYgW+++47DBo0CFevXq13YOVXX32FmTNnIicnx6ZWmzVr1uCzzz6zeYY9ooawJk9EVIvRo0dj8ODBiI2NRWRkJMLCwhq8c+Kxxx7D6dOn8fvvv9t033/jxo2xatUqe0MmqoE1eSIPxpo8kXtjkiciInJTnNaWiIjITTHJExERuSkmeSIiIjfFJE9EROSmmOSJiIjcFJM8ERGRm2KSJyIiclNM8kRERG6KSZ6IiMhNMckTERG5qf8H4BWYtutQ9DkAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 512x384 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Polarization fraction as a function of the dust column density\n",
    "PvT_plot_Step2 = HtoH.PvT(CatP, CatT, Tscale=[15,30], Pscale=[0.0,30.0], CatdP=CatdP, errorbars='true', Source=FullName)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "491d0620",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Saving figures\n",
    "# PNG\n",
    "PvI_plot_Step2.savefig('Plots/'+SourceName+'_PvI_plot_Step2.png',dpi=300)\n",
    "PvN_plot_Step2.savefig('Plots/'+SourceName+'_PvN_plot_Step2.png',dpi=300)\n",
    "PvT_plot_Step2.savefig('Plots/'+SourceName+'_PvT_plot_Step2.png',dpi=300)\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "PoLiteWIP",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
